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Intergalactic space is filled with a pervasive medium of ionized gas, the Intergalactic Medium 
(IGM). A residual neutral fraction is detected in the spectra of Quasi-Stellar Objects at both 
low and high redshifts, revealing a highly fluctuating medium with temperatures characteristic 
of photoionized gas. The statistics of the fluctuations are well-reproduced by numerical gravity- 
hydrodynamics simulations within the context of standard cosmological structure formation sce- 
narios. As such, the study of the IGM offers an opportunity to probe the nature of the primordial 
density fluctuations on scales unavailable to other methods. The simulations also suggest the IGM 
is the dominant reservoir of baryons produced by the Big Bang, and so the principal source of 
the matter from which galaxies formed. The detection of metal systems within the IGM shows 
that it was enriched by evolved stars early in its history, demonstrating an intimate connec- 
tion between galaxy formation and the IGM. The author presents a comprehensive review of the 
current understanding of the structure and physical properties of the IGM and its relation to 
galaxies, concluding with comments on prospects for furthering the study of the IGM using future 
ground-based facilities and space-based experiments. 



Contents 

I. INTRODUCTION Q 

II. Observations 

A. Resonance absorption lines 

B. Absorption line properties ^ 

1. Physical properties of absorption line systems ^ 

2. Evolution in the number density of the Lyo forest 13 

3. Characteristic sizes and spatial correlations 1^ 

4. Deuterium absorption systems 3 

5. Helium absorption systems 3 

6. Metal absorption systems 1^ 

C. Direct flux statistics 

1. Mean transmitted flux , 

2. Statistics based on pixel fluxes 1^ 

III. Ionization Equilibrium 

A. Ionization 

1. Hydrogen and helium 

2. Metals 

B. Thermal equilibrium 

IV. The Metagalactic UV Background 

A. Mean energy density of the UV background 

1. Origin of the UV background 

2. Sources 

B. Fluctuations in the UV background 

C. Observational constraints on the UV background 

1. The proximity effect 

2. Constraints from helium ionization 

3. Constraints from metal ionization 



VI. Numerical Simulations 

A. Cosmological structure formation 

1. Dynamical evolution of the dark matter and 
baryons 

2. Methods of numerical simulations 

B. Physical properties of the intergalactic medium 

C. Statistical properties of the absorption systems 

1. Mean metagalactic ionization rates 

2. Spectral properties of the IGM 

3. Intergalactic helium absorption 

4. Metal absorption systems 

5. The flux power spectrum 

VII. Reionization 

A. The epoch of reionization 

B. The growth of cosmological H II regions 

C. Sources of reionization 

1. Galaxies 

2. QSOs 

3. Other sources 

D. 21cm signature of reionization 

VIII. The Absorber— Galaxy Connection 

A. Galaxy-associated IGM absorption 

1. Association of DLAs with galaxies 

2. Association of metal absorbers with galaxies 

B. Galactic winds and the IGM 

IX. Prospects for the Future 
Acknowledgments 
References 



V. Early Models 

A. Pressure-confined clouds 

B. Dark matter minihalos 

C. Caustics and sheets 

D. Galactic models 



'Submitted to the Reviews of Modern Physics. 
t Scottish Universities Physics Alliance 
t Electronic address: lA.Meiksin@ed.ac.ukl 



I. INTRODUCTION 

According to the Big Bang theory, the primordial hy- 
drogen and hehum created in the Universe first materi- 
ahzed in the form of an extremely hot ionized gas. By 
the time the Universe was three hundred thousand years 
old, adiabatic expansion cooling brought the tempera- 
ture of the primordial plasma down until the hydrogen 
and helium recombined. The radiation last scattered at 
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this time appears today as the Cosmic Microwave Back- 
ground (CMB). The search for the IGM began w ell be- 
fore the discovery of the CMB with an attempt by iFieldl 
( 1959al ) to detect the hyperfine 21cm absorption signa- 
ture from hydrogen along the line of sight to the extra- 
galactic radio galaxy Cygnus A. Although no detection 
was made, combining the optical depth with the mea- 
sured temperature of the CMB discovered in 1965 would 
have been sufficient to exclude the possibility that the 
Universe was closed by baryons, with an upper limit on 
the baryon density of only 20% of the closure density 
required for an Einstein-deSitter (flat) universe. ^ 

Nearly coincident with the discovery of the CMB, how- 
ever, a considerably improved measurement of the den- 
sity of intergalactic neutral hydrogen was made. Soon af- 
ter the dis covery of th e first Q uasi-Stell ar Object (QSO) 
Hazard et a l. 1963: Schmidtl . ir963 ). G unn and Petersonl 



19651 ) reported a small decrement in a QSO spectrum 
shortward of its Lya emission line. Attributing the 
decrement to the Lya resonance scattering of radiation 
from the QSO by intergalactic neutral hydrogen, Gunn 
and Peterson demonstrated that the cosmic mass den- 
sity of neutral hydrogen was exceedingly smaller than the 
Einstein-deSitter density. In fact, it was far smaller than 
the spatially averaged hydrogen of all the stars in the 
Universe. If the Big Bang theory was correct, it meant 
either that galaxy formation was an extraordinarily effi- 
cient process, sweeping up all but a tiny residue of the 
primordial hydrogen, or that the gas was reionizcd. 

These two themes, the detection of intergalactic gas 
through the 21cm signature in the radio or through Ly- 
man resonance scattering lines in the optical or ultravio- 
let (UV) , continue to dominate studies of the Intergalac- 
tic Medium (IGM). To date, almost all that is known 
about the structure of the IGM has been derived from op- 
tical and UV data. This situation is expected to change 
dramatically in the near future with the development of 
low frequency radio arrays like the LOw Frequency AR- 
ray (LOFAR)^ the Murchison Widefield Array (MWA)^, 
the Primeval Structure Telescope (PaST/21CMA) ^, the 
Precision Array to Probe the Epoch of Reionization (PA- 
PER) ^, and a possible Square Kilometre Array (SKA). ^ 
A primary science driver of these instruments is the direct 
imaging of the I GM prior to the com pletion of the Epoch 
of Reionization ( Madau et aLl . ll997h . Most of this review 
focuses on the current understanding of the state of the 
IGM as determined from optical and UV measurements. 



^ This conclusion requires making the (at the time) reasonable as- 
sumptions that the hydrogen was all neutral and that the hyper- 
fine structure levels of the hydrogen were in thermal equilibrium 
with the CMB. 

^ http://www.lofar.org 

3 http:/ /www. haystack. mit.edu/ast/arrays/mwa/ 

* http://web.phys.cmu.edu/~past/ 

^ http:/ /astro. berkeley.edu/~dbacker/eor/ 

^ http://www.skatelescope.org 



The observations have relied almost exclusively on the 
spectra of QSOs, although IGM absorption features have 
also been detected in the spe ctra of Gamma-Ray Bursts 
(GRBs) (Tota ni et oil . \2Q0&), and indeed played a key 
role in establishing the extr agalactic character of some 
bursts (|Metzger eii aLl . ll997l ). 

Almost immediately after Gunn and Peterson pub- 
licized their finding, it was recognized that individ- 
ual Lya absorption features should appear from neu- 
tral hydrogen concent r ated i nto cosmo l ogical structures 
()Bahcall arid Salpeteil . 1 19651 : IWagoneil . Il967t) . Absorp- 
tion features had in fact been detected in higher resolu- 
tion QSO spectra, but these were identified with inter- 
vening ionized metal absorption systems fiBahca ll et aZj, 
Il968l ). as was expected if galaxies had hot halos of ion- 
ized gas: the lines of sight to background QSOs were 
expected to pass through such hot halos an d intercept 
any i onized gas clouds within them (Bahcall and Spitzed . 
I1969D . The features, however, were uncomfortably com- 
mon, hinting at a class of unknown structures not as- 
sociated with galaxies. The Lya resonanc e line features 
continued to prove elusive until 1971, when lLvndsl ( 197ll ) 
recognized several absorption features shortward of the 
Lya emission line of a QSO as Lya lines^ arising in a 
population of discrete absorption systems also showing 
metal features. The Lya lines form a plethora of distinct 
absorption features in the spectra of high redshift QSOs; 
they are collectively known today as the Lya forest. 

The properties of the Lya forest came under increas- 
ing scrutiny, with th e first systema, t ic sur vey conducted 
several years later bv lSargent et all i 1980f) . convincingly 
demonstrating through the homogeneity of the observed 
properties of the absorbers their intergalactic origin, as 
opposed to clouds ejected by the QSOs observed. Al- 
though limited by the resolution of the spectrograph, 
the measured velocity widths of the Lya features cor- 
responded to gas temperatures of a few to several 10'' K, 
characteristic photoionization temperatures for gas of a 
primordial composition, for which there is no significant 
cooling by metals. The number of features per comoving 
length was shown to increase with redshift, demonstrat- 
ing that the systems were evolving (Young et al. . ,1982) . 

The past decade has witnessed a dramatic improve- 
ment in precision studies of the Lya forest with the ad- 
vent of 8-10 m class telescopes, particularly the Keck 
10-m and the 8.2m Very Large Telescope (VLT). For 
the first time, the individual absorption features in the 
Lya forest were spectroscopically resolved over their full 
range. Velocity widths of ~ 25 km s^'^ are typical. The 
neutral hydrogen column densities of the absorbers range 
from roughly 10'^ — 10^^ cm~^. The highest column den- 
sity systems, the Damped Lya Absorbers (DLAs), are 



These are not true absorption features involving the net destruc- 
tion of a photon, but the scattering out of the line of sight of 
resonance line photons. 
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of particular interest for galaxy formation, as they are 
suspected of containing the neutral gas that formed the 
bulk of the stars in present day galaxies. 

As the number per length of absorption systems in- 
creases along a line of sight with increasing redshift, so 
does the mean flux decrement in a background QSO spec- 
trum due to Lya scattering. The QSO spectra measured 
as part of the Sloan Digital Sky Survey (SDSS) * show 
a rapid rise in the flux decrement at 5.5, suggesting 
that the epoch of H i reionization may lie not far beyond 
z ~ 6. Many of the hydrogen features also show absorp- 
tion lines from metals, including carbon, silicon, nitrogen, 
oxygen, magnesium, iron and others. The abundances of 
the metals, however, are at most about 10% of solar at 
low redshifts, and as low as 1% at high redshifts, indi- 
cating that the absorption systems are comprised largely 
of primordial material. The primordial nature of the gas 
received furt her important confirma tion in 1994 with the 
discovery bv IJakobsen et ail (1993) of intergalactic he- 
lium using the Hubble Space Telescope (HST) through 
the detection of He ii Lya absorption. 

Because the baryons in the IGM are detected only 
though their absorption signatures, the physical struc- 
tures that give rise to the features must be modeled. 
Early models characterized the systems as discrete clouds 
of gas, with most of the fo cus on either clouds pressure- 
confined by a hot medium feargent et all 119801) . or grav- 
itationally c onfined in a dark matter halo ( Ikeuchil . Il986l : 
iReesl . Il986l ). At the time it was believed that the ab- 
sorption systems accounted for only a few percent of 
the baryons produced in the Big Bang, much like galax- 
ies, their visible counterparts. A paradigm shift in the 
models occurred in the mid 1990s. Measurements of co- 
incident absorption features along neighboring lines of 
sight suggested si zes of tens to hundreds of kiloparsecs 
for the absorbers (Bechtold et all . Il994t iDinshaw et all . 
119971: ISmette et alllldd^ much larger than expected for 
clouds confined by pressure or dark matter halos. A 
radical transformation in the understanding of the na- 
ture and structure of the IGM was initiated by numeri- 
cal simulations of cosmological structure formation. To- 
day essentially all of the baryons produced in the Big 
Bang are believed to have undergone the same gravita- 
tional instability process initiated by primordial density 
fluctuations that was responsible for the formation of 
galaxies. The computation of the structure of the IGM 
has been converted into an initial value problem simi- 
lar to that of the CMB fluctuations. Fluctuations in the 
CMB are solved for by following the gravitational insta- 
bility growth of a postulated spectrum of primordial dark 
matter density fluctuations. The growth of structure in 
the IGM is now treated as the nonlinear extension of 
these computations. The result is a networ k of filamen- 
tary structures, the so-called "cosmic web" ( Bond e^"aLl . 



Il996l ). The Lya forest is believed to be a signature 
of the cosmic web. Early simulations broadly repro- 
duced the statist ics of the Lya forest spectacularly well 
n et al, 1994: Hcrnquist et aLl. 11996: Petitiean all . 



|199: 



1995: Theuns et al . 1998; Zhang et aL .1995. . ,19971 ). An 



* http://www.sdss.org 



immediate conclusion was that at z ^ 1.5, some 90% of 
the baryons produced in the Big Bang are contained 
within the IGM, with only 10% in galaxies, galaxy clus- 
ters or possibly locked up in an early generation of com- 
pact stars. 

Soon after the discovery of intervening absorption fea- 
tures, it was recognized that they provided potentially 
powerful tests of fundamental properties of the Universe. 
The split in the fine structure lines of the metals was used 
to set co nstraints on the vari ability of the fine structure 
constant ( Bahcall et aZ.I . IT967l ). A bunching of absorption 
features near z ~ 1.9 (now known to be fortuitous), gave 
rise to the (re) introduction of a cosmological constant 
to account for the numerous features as multiple images 
due to lensing (jPetrosian et aLl . [19671 ). The expected pri- 
mordial composition of the IGM offered the potential 
of placing constraints on the photon-to-baryon ratio of 
the Universe through measurements of the intergalactic 
deuterium abundance. More recently, the success of the 
models has inspired attempts to exploit the Lya forest 
as a new means to constrain cosmological structure for- 
mation models and obtain stringent constraints on the 
cosmological parameters. 

The description of the IGM by the simulations, how- 
ever, is far from complete. There remain many unsolved 
problems. The current simulations do not reproduce all 
the observed properties of the IGM. The absorption lines 
are predicted to be substantially narrower than mea- 
sured. This likely stems from the principal outstanding 
missing piece of physics, the reionization of the IGM. Not 
only must hydrogen be ionized, but helium as well. The 
ionization heats the gas through the photoelectric effect. 
Detailed radiative transfer computations are required to 
recover the temperatures, for which there is still limited 
success. The sources of the reionization and the epochs of 
reionization, both of hydrogen and helium, are still not 
firmly established. The origin of the metal absorption 
systems in the diffuse IGM is still unknown, although it 
is widely expected they were deposited by winds from 
galaxies, possibly driven by intense episodes of star for- 
mation. As such, the metal absorption lines in principle 
offer an important means of studying the history of cos- 
mic star formation. Most fundamentally, the relation of 
the IGM to the galaxies that form from it is still mostly 
unknown, but offers perhaps the most exciting prospects 
for new research directions. 

The purpose of this review is to describe the progress 
made in the understanding of the origin of structure 
within the IGM, with a view to presenting the underlying 
physics that determines the structure. An understanding 
of the physics is necessary for future progress. The past 
decade has revealed the IGM to be a complex dynam- 
ical arena involving interactions between the IGM gas. 
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FIG. 1 Number of refereed papers on the intergalactic 
medium and QSO intervening absorption systems published 
from 1965-2007, as provided by the ADS abstract service 
(adsabs.harvard.edu). The sharp rises follow key develop- 
ments in astronomical instrumentation and observations. The 
lower histogram shows the rise of papers on the reionization 
and subsequent ionization state of the IGM. 



galaxies and QSOs. It is becoming increasingly appar- 
ent that the separation of these systems into distinct and 
isolated entities is an artificial construct. Galaxies and 
QSOs originated from the IGM, and their radiation and 
outflows impacted upon it. Any complete understanding 
of the origin of these systems requires a unified treatment. 
In this way, the IGM resembles the interstellar medium 
of disk galaxies in which the gaseous component is inti- 
mately linked to the stars and their evolution and im- 
pact. Interpreting the increasingly refined observations 
requires detailed modeling, which relies on large-scale nu- 
merical computations involving gas, radiative processes, 
and gravity. The physics involved is intermediate in com- 
plexity between that of the CMB and galaxy formation, 
rendering the IGM a bridge between these extreme scales 
of cosmological structure formation. Unraveling the pro- 
cesses that led to the formation and structure of the IGM 
may thus serve as a crucial step in the solution of the 
much more involved problem of galaxy formation. 

A search of the literature for papers on the IGM since 
the Gunn & Peterson (1965) measurement produces close 
to 6000 references.^ While this review does not have the 
space to describe the observational methods used to mea- 
sure the IGM, it should be recognized that progress in the 
understanding of the properties of the IGM is indebted to 



^ Based on a boolean search of the Astrophysics Data System ab- 
stract service (adsabs.harvard.edu) for all refereed papers with 
abstracts or keywords containing "(intergalactic and medium) or 
(quasar and absorption and line)." Reionization papers are se- 
lected as the subset discussing the reionization of the IGM or the 
subsequent ionization structure of the IGM. 



advances in observational techniques. This is well illus- 
trated by a plot of the number of refereed publications in 
the field against time. Periods of relatively steady output 
are punctuated by four leaps. The first occurs at the end 
of the 1960s and beginning of the 1970s with the introduc- 
tion of image tube spectrographs coupled with integrat- 
ing television systems for photon counting ( Boksenberd . 
|l972; Morton, 1972), which greatly facilitated the taking 
of spectra. The next occurs in the mid-1970s with the de- 
velopment of x-ray astronomy following the launch of the 
UHURU satellite in 1970 and the recognition that galaxy 
clusters contain a n extended and pervasive medium o f 
hot, radiating gas ( Kellogg et aLl . [19731 : iLea et al. 
Another sharp rise occurs in the early 1990s following the 
launch of the Hubble Space Telescope in 1990, the instal- 
lation of EMMI and its echell e spectroscope on the New 
Technology Telescope (NTT) (lD'Odoricol.ll990l). and the 
delivery of the HIRES spectrometer (|Vogt et all Il994l ) 
to the Keck telescope. The fourth occurred in 2000 with 
the introduction of the UV and Visible Echelle Spec- 
trograph (UVES) to th e Very Large Telescope (VLT) 
(|D'Odorico et ali |2000[ ). the launch of the Far Ultrav i- 
olet Spectroscopic Explorer {FUSE) (jMoos et ali l200Clf) . 
and the beginning of oper ations of the Sloan Digital Sky 
Survey (|York et aLl . boOOD . The latter in particular trig- 
gered a surge of activity in reionization studies following 
the discovery of high redshift QSOs (z > 6) with spec- 
tra indicating a rapid rise in the effective optical depth 
of the IGM to Lya photons, hinting that the Epoch of 
Reionization wa s being approached (Beck er et al\ . l200ll : 
iFan et all l2002f) . This is indicated by a rise in IGM 
reionization papers in Fig. [Jl a trend which continues 
today, fueled in part by the growing interest in the influ- 
ence of reionization on the CMB fluctuations measured 
by the Wilkinson M icrowave Anisotropy Probe ( WMAP) 
(|Kogut et aZ.l . l2003l ). The next major leap may well come 
from the anticipated radio measurements. 

The rapidly rising tide of IGM studies has brought 
along with it several reviews. This is fortunate, as it 
is impossible in a single review to cover all areas thor- 
oughly today. Early reviews of the QSO absorption 
line literature, now largely histor ical, were provided by 
IStrittmatter and Williamsl (1 19761) an d Wevmann et al\ 
(19811). iRauchI (^1998^ and iBechtoldl (I2003D ])rovide re- 
views of the Lya forest that are still largely up-to- 
date in the sense that most of the topics currently en- 
gaging the community are treated, with developments 
since mostly improvements in accuracy and in the de- 
tails of the numerical models. Reviews of the low red- 
shift IG M are provide d bv I Shull (2003. ) and lStocke et al\ 
(|2006bl ). IWolfe et al\ ([20051) have reviewed the current 
understanding of Damped Lya Absorbers, an absorber 
class of special concern as it represents the closest link 
to the gaseous component forming present day galaxies. 
A new series of reviews followed the recent explosion of 
activity on the reionization of the IGM in anticipation of 
the detection of the Epoch of Reionization (EoR) through 
its high redshift 21cm signature. An early review in this 
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area is by iLoeb and Barkanal ( 200lt ). Since then, cur- 
rent observational and theoretical as pects of rei o nizatio n 
have been exhaustivel y cov ered by Fan et al. ( 2006cJ) . 
iFurlanetto et al\ 12006), and lBarkana and Loebl (|2007l) . 

Rather than repeating the wide range of IGM phe- 
nomenology already covered by previous reviewers, I con- 
centrate here on the physics underlying the structure of 
the IGM. One aim is to describe the physical underpin- 
nings of current numerical simulations as required for fu- 
ture simulations to further progress. As observations are 
crucial for constraining any models, I begin by giving 
a broad overview of the current observational situation. 
The next section describes the physics of ionization equi- 
librium, followed by a discussion of the UV metagalactic 
background that maintains the ionization of the IGM. A 
brief review of early models of the Lya forest absorbers 
is then presented, followed by a discussion of current nu- 
merical models. The reionization of the IGM is then 
summarized along with means for its detection. This is 
followed by a discussion of the connection between galax- 
ies and the IGM before concluding with observational 
and theoretical prospects for the future. 

Unless stated otherwise, the cosmological parameters 
— 0.3 and 0„ = 0.7 are assumed for, respec- 
tively, the ratios of the present day matter density and 
the present day vacuum energy density to the current 
Einstein-deSitter density Pciit(O) = 3Hq/8ttG, a Hub- 
ble Constant of Hq ^ lOO/ikms^^ Mpc^i with h = 0.7 
(Mpc = Megaparsec ~ 3.0856 x 10^^ m), and a baryon 
density of ilf,/i^ = 0.022. These values are consistent 
within the errors with the current best estimates for a 
flat Universe based on CMB measurements using Year-5 
WMAP data of f7,„ 0.279 ± 0.014, = 0.721 ± 0.015, 
h = 0.701 ± 0.013, and flbh'^ = 0.02265 ± 0.00059 



(1) 



( Komatsu et all |2008| ) , or intergalactic D/H measure- 
ments, giving f^^/i^ ^ o.021 ± 0.001 (|0'Meara et al\ . 
120061 ). 



II. OBSERVATIONS 

A. Resonance absorption lines 

The IGM is detected through the absorption features it 
produces in the spectrum of a background bright source 
of light (typically a QSO). The production of the ab- 
sorption features is governed by the equation of radiative 
transfer through the IGM, conventionally expressed in 
terms of the specific intensity of a background radiation 
source. 

The specific intensity /i/(r, t, n) is defined as the rate at 
which energy crosses a unit area per unit solid angle per 
unit time as carried by photons of energy hpv traveling 
in the direction n relative to some fiducial direction. The 
equation of radiative transfer for I^{r, t, h) in a medium 
with attenuation coefficient a^{r,t,n) is 



ldl^ir,t,n) 
c dt 



Because of its use below, a radiative source term has 
been included through the emission coefficient j^{r, t, n), 
which describes the local specific luminosity per solid an- 
gle per unit volume emitted by sources. Normally the 
random orientation of atoms will ensure ji, is isotropic, 
but not always, as for instance if the source includes a 
scattering term, so that a dependence on n is included 
to account for the possibility of anisotropic sources. A 
central source like a QSO or galaxy will in fact generally 
radiate anisotropically. 

In general, the attenuation of the radiation field is due 
both to the absorption of photons and their scattering 
out of the beam. The attenuation coefhcient is then 

a^{r, t, n) = p{r, t)K^{r, t, n) + n(r, t)a^{r, t, h), (2) 

where p(r, t) is the mass density of the medium, 
Ki/(r, t, n) is the opacity, which expresses the absorption 
cross section per unit mass, n{r,t) = p(r,t)/m is the 
number density of scattering particles of mean mass ifi, 
and iTi/(r, t, n) is the scattering cross section. In general, 
the scattering could arise from more than one type of par- 
ticle, in which case nai, would be replaced by a sum over 
particle species i of density m and cross section a^ii)- 
In a static medium, a^{r,t,n) is generally isotropic, but 
not always. As an example, the alignment of atoms in 
a strong magnetic field would absorb anisotropically. In 
a moving medium, an anisotropic contribution to the at- 
tenuation will be produced by the dependence of Doppler 
scattering on the relative motion of the radiation and the 
fluid. 

The value of Ii, at any given time t and position s along 
the direction n will be given by any incoming intensity 
I™'^ at position sq at time tret = t — {s — so)/c, attenuated 
by intervening material at positions s' at the retarded 
times tj.gj = t — {s — s')/c, along with contributions from 
sources at positions s" that emitted at the retarded times 
t"(,j — t~ {s~ s")/c, followed by attenuation. The formal 
solution to Eq. ^ is then 



exp 



+ / ds" {j,^)s",t'j^^ exp - f 

<y Sq \_ J s' 



ds' {a^)s',t'^^^ 
ds' (a^)s',t' 



(3) 



as may be verified by direct substitution. In a uni- 
formly expanding (or contracting) medium described by 
the scale factor a{t), Eq. ([T]) must be modified to 



ldli,{r,t,h) ^ Id 
c 



dt 



c a 



31„{r, t,h) — v 



dl^{r,t,h) 
dv 



+n • V/,. (r , t , n) = -a,, (r , i , n)/^ (r , i , n) + > (r , i , n) , (4) 

where the frequency redshifts according to cx a{t)~^ . 
The solution becomes 



n • V/,.(r, i,n) 



"a(io)' 


3 




[ a{t) \ 


exp 


[-/■ 

L J So 



ds' (a^')s',t' 
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a{t) 



exp 



ds' {a^')s',t' 



,(5) 



where vq = i^o,{t) / o.{to), v' = va{t) / a{t')^ and v" = 
va{t)/a{t"). 

In the case of no intermediate sources, ju = and the 
received intensity depends only on how the incident in- 
tensity is modified by the intervening gas. This is the 
situation for a single background source such as a QSO. 
The Lyman absorption features arise from the scattering 
of resonance line photons received from the background 
QSO through a medium of scatterers with number den- 
sity n{Y,t). The frequency dependent scattering cross 
section for a scatterer at rest is given by the Lorentz pro- 
fileio 



(-) 


1 


\mecj 


47reo. 



flu 



(r„i/47r2) 



(6) 



where vi^ is the resonance line frequency, correspond- 
ing to the wavelength A;„, of the transition be- 
tween an upper level broadened by radiation damp- 
ing to a sharp lower level (the ground state), V^i = 
(5i/iu/ffn)[l/47reo]87r2e^/(mecA^„) is the damping width 
of the upper level, where gi and gu are the respective sta- 
tistical weights of the lower and upper levels, //„ is the 
upward oscillator strength, e is the electric charge of an 
electron, and mg the electron mass.^"'^ The cross section 
is normalized according to 



dvon — 



Tie 

rUeC 



1 



47reo 



flu ■ 



(7) 



For Xiu in m, T^i = 6.670 x 10'^ {gi fin /9u)XTu- For the 
Lyman series. 



2n-4 



51 /i 



2^n^{n- 1) 
3(n+l)2"+4 ' 



(8) 



where n is the q uantum number of the upper state 
(|Menzel and Pekeris . 1935.). For the Lya transition, 
gi = 2, gu — 6 and fiu = 0.4162 • • •. For hydrogen Lya, 
= 1215.67A and F^; = 6.262 x 10^ s'K For He ii 



^ SI units for electrostatic quantities throughout the text may 
be accommodated by including the square-bracketed quantity 
[l/47reo]i where eq is the vacuum permittivity. Without the fac- 
tor, the expressions correspond to the forms appropriate for cgs 
units. 

The Lorentz profile neglects the frequency dependence of the 
scattering cross-section far from line center. It does not, for 
example, recover the low frequency Rayleigh limit (ctj/ oc u"^). 
The profile, however, is fully adequate for all practical applica- 
tions to the IGM. Small deviations may become detectable for 
broad absorption featur es such as Damped Lya Absorbers at 
high spectral resolution l|Led . [20o3 ). Such deviations, however, 
would prove difficult to distinguish from the effects of density 
inhomogeneities in the absorbing gas and the associated gas ve- 
locities, except possibly in a statistical sense over many systems. 



Lya, Xiu is smaller by a factor of 4, and r,j larger by a 
factor of 16. 

In general, the atoms will not be at rest. At the very 
least they will generally undergo thermal motions de- 
scribed by a Maxwellian velocity distribution correspond- 
ing to their temperature T. They may additionally take 
part in an ordered flow of velocity v, but this may be 
accounted for by shifting the line center frequency v to 
1^(1 ~ ^^||/c), to first order accuracy in v/c, where is 
the line-of-sight velocity. In addition, there may be a 
disordered turbulent, or so-called microturbulent, com- 
ponent in some collapsed or shocked regions. Ignoring 
non-thermal velocities, the profile including thermal mo- 
tions is found by convolving Eq. © with a Maxwellian. 
The result, which neglects additional frequency depen- 
dences far from line center, is the Voigt function 



i?(a, x) = — 
n 



dy- 



yf 



(9) 



where a — r„i/(47rAz/u) is the ratio of the damping 
width to the Doppler width, and x = {y — viu)II^V£, is the 
frequency shift from line center in units of the Doppler 
width Ai/D = viuhj c with Doppler parameter 



2fcBT 



1/2 



(10) 



where ma is the mass of the scattering particle and /cb 
is Boltzmann's constant. A kinematic component such 
as microturbulence is often accounted for by adding in 
quadrature the characteristic kinematic velocity to the 
thermal component of the Doppler parameter 



1/2 



(11) 



where 6th is the thermal contribution from Eq. (jlOp . An 
expansion a pproximation for the Voigt function in a is 
provided by iHarris (1948'). Fast m ethods f or arbitrary a 
are described by Humli'cck (1979'). iHumlf cck (1982) and 
IShipponv and ReadI (Il993l) . Various evaluation methods 
are compared by Schreieil (|1992( ). 

The resonance line scattering cross section is 



7re 



47reo 



jiu'^v{a,v) ^ a<^v{a,v), (12) 



where 



(^y (a, v) = -^—<i)v{a, x) = i?(a, x) (13) 



rOO 



is the Voigt profile, normalized to dv tfv{o-:^) = 1- 
For hydrogen Lya , a ~ 0.0472r"-^/2. For He ii Lya, a 
is larger by a factor of 8. 

The intensity of the background QSO is attenuated by 
the factor e~'^", where Eq. ([5]) gives 



= ds' n{s' , t')ay 

J Sn 



(14) 
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For a homogeneous and isotropic background Universe 
expanding with scale factor a{t), radiation emitted by 
the source at time to and rest frame frequency vq > viu, 
where i^i^ is the resonance hne frequency of a (rest frame) 
Lyman transition, will be scattered by the medium at 
epoch t' given by i^' ~ i'Qa{to)/a{t') — viu- Thus all 
of the received QSO intensity will be attenuated at the 
time t at all observed frequencies v in the range > v > 
viua{to) / o,{t) , or wavelengths A in the range \iu < X < 
Xiua{i) / aitfj) , where A;^ is the (rest frame) wavelength 
of the particular Lyman transition, producing a trough 
in the spectrum of the QSO. This is the Gunn-P eterson 
effect ( Gunn and Petersonl Il965t IScheueil . Il965l ). Inter- 
vening metal absorption, assuming solar abundances dis- 
tributed uniformly throughout intcrgalactic space , would 
produce similar troughs (Shklovskij 1 1964 Il965l ). The 
corresponding optical depth is given by noting that ds' = 
cdt' = cH-^{t')da{t')/a{t'), where Hit') = a{t')la{t') is 
the Hubble parameter at time t' . Changing the integra- 
tion variable to a' = a{t') and setting z = viu/v — 1, it 
follows from Eq. ([T^ that, taking (n) to denote the mean 
number density of particles in the lower state. 



/do,' 
-^H{a')-^a,,a'-,Ma')) 

9u 1 Af„r„i(n(z)) 



91 Stt H{z) 



1.191 X 10 



[l + zf 



[n„^{l + Zf + fix(l + Zf + r!„]l/2 

^ f gu/gi\ ( (n 



506. OA 



(15) 



after noting from Eq. ([7]) that a — {gu/ gi)TuiXfy^/8'K. 
The Hubble parameter H{z) at redshift z is related to 
the value of the Hubble constant today Hq by H{z) = 
HnE(z) where E(z) = \nrr,.(l + z f + nK(l + z)"^ + 
(jPeeblesl . Il993l : ISargent et d\ . Il980[ ). Here, f7,„ is the 
mass density of the Universe expressed as a fraction of 
the Einstein-de Sitter critical density, and Qk and Qy arc 
similar fractions arising from the curvature and vacuum 
energy contributions, respectively. They satisfy the iden- 
tity Ujn + ^K + ^v = 1. In the last line of Eq. (fT5|) . the re- 
sult was scaled to the present day mean cosmic hydrogen 
density nH(0) = pH(0)/mH = (1 - F)nf,Pcrit(0)/mH ^ 
0.186 m~3^ where mn is the mass of a hydrogen atom 
and assuming a baryon density of fl^h^ = 0.022 (see 
§ III. B. 41 below) and cosmic helium ma ss fraction Y ~ 
0.2477 ± 0.0029 (|Peimbert et all I2007D . and to a Hub- 
ble Constant of Hq = 70kms-i Mpc"! IILB.2I below). 
Eq. (|15p is the inverse of the Sobolev parameter for a ho- 
mogeneously and isotropically expanding medium, and 
was first derived in the context of cosmological hydrogen 
Lya photon scattering bv . Field (1959b ). 

Observations show that in fact the IGM is not homo- 
geneous, but that the baryons have coalesced, creating a 
fluctuating density field n{s,t). The discreteness of the 
absorption lines measured in the spectra of high redshift 
QSOs suggests they originate in distinct localized regions. 



The resulting collection of absorption features, shown in 
Fig. [21 is known as the Lya forest. Denoting the centers 
of the absorbing regions by the positions Si, the optical 
depth becomes Ti, — J2i '''v{i), where 



i+As/2 

Tu(i) = I ds' n(s\ti)ay 

'si-As/2 



(16) 



Here, the integral has been localized to a region of width 
As around Si, and i'' — va{t) / a{ti), where ti corresponds 
to the epoch of coordinate Si in an expanding universe. 
The absorption features thus probe discrete spatial struc- 
tures in the IGM. Expressing the optical depth as a 
sum of discrete absorption systems is, however, an ap- 
proximation, as the intcrgalactic gas forms an evolving 
spatial continuum. Gas distributed over a wide spatial 
range, and even at different epochs, can therefore affect 
the same observed frequency u. However, for z^3.5, 
observations show the discrete absorber approximation 
provides a reasonably good description of the hydrogen 
Lya forest. Eq. pB)) may be further simplified to 



.^1/2 



where 



NaX 



TOgC 



47ren 



(17) 



(18) 



is the optical depth at line centre. Here, N is the col- 
umn density J ds' n{s' ,t'), and {(pv{a,x)) is averaged 
over the line of sight, weighted by the density. This may 
differ from 0y(a, x) when the temperature or large-scale 
macroscopic velocity field varies along the line of sight. 
Tabulations of A/;„ and related atomic data for resonance 
lines for a variety o f elements are available in the litera- 
ture (|Mortonl . [T99ll . 2000, 200l). For hydrogen Lya , 



-HI 



0.38 



101 



20 km s- 



(19) 



At redshifts z > 3.5, the blending of absorption lines 
makes it increasingly difficult to attribute the absorp- 
tion to distinguishable systems. By z > 5.5, the indi- 
vidual lines have essentially all merged, forming in effect 
a Gunn-Peterson trough (though still not corresponding 
to a completely neutral IGM). Absorption features at in- 
termediate redshifts not easily deblended into subcompo- 
nents may be characterized by an equivalent width wx- 
This is defined to be the width, expressed in units of 
wavelength, a square-well absorption feature with zero 
flux at its bottom must have to match the integrated 
area of the detected feature under the continuum: 



Wx 



dX{l 



\2 roo 
C Jo 



(20) 

For a Voigt profile, the equivalent width is related to 
the column density through the line-center optical depth 
according to 



^ = 2^na,ro), 
A c 



(21) 
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FIG. 2 The Lya forest measured in the spectrum of the z = 3.22 quasar HS0741 + 4741. The spectrum was obtained using the 
HIRES spectrograph on the Keck 10 m telescope. The integration time was 4.0 hours with a spectral resolution of -R = 37, 000 
covering 3900A to 66OOA with a signal-to-noise ratio of 145 per resolution element. The sharp peak at A = 5130A is the 
Lya emission line of th e quasar. The Lyg fore st lies towards the blue side of the feature. A Damped Lya Absorber at z = 3.02 
appears at A ~ 4890A (iProchaska et all\200l\ ). Intervening C iv systems give rise to rarer absorption features longward of the 
emission line l|Songailal . l2006l ). (Figure courtesy of A. Songaila.) 



where 



F(a,ro) 



dx 



|l-exp -TT'^/^To<j>v{a,x) |. (22) 



The relationship between the equivalent width and the 
column density is known as the "curve of growth." For 
To < 10^, the absorption profiles are well- approximated 
as Doppler in shape, corresponding to (/)y(a, x) ~ 
7r~-'^/^e^^ , and F{a,To) may be denoted more simply 
by F{to). T he funct i on F (to) may be expressed as the 
power series (|Munchl . [l968h 



, '-JKJ 



(23) 



For very optically thin lines (tq << 1), F{to) ~ 
(A/7r/2)ro, and the equivalent width is related linearly 
to the column density, 



X 



NXfi^ ~ 8.85 X lO-'^NXfi^, (24) 



where iVA has the dimension cm ^. Such features are 
said to lie on the "linear" part of the curve of growth. 

For 6 < To < 300, bette r than 1% acc uracy is provided 
by the asymptotic series ( Miinchl . Il968l ) 



i^(To) - (In To) 



1/2 



0.2886 0.1335 



In To (hiTo)^ 



(25) 



For H I Lya at a temperature T — K, F(a,To) de- 
viates from the logarithmic approximation Eq. (j25p by 
less than 10% for 2 < To < 5000. Such features are said 
to lie on the "logarithmic" part of the curve of growth. 
They are also referred to as "saturated" lines, since an 
increase in the column density has only a small change 
on the shape of the absorption feature. The measure- 
ment of accurate column densities is notoriously difficult 
for these features, necessitating the search for higher or- 
der resonance lines on a more linear part of the curve 
of growth if accurate column densities are desired. On 
the other hand, the equivalent width is nearly directly 
proportional to the Doppler width, allowing an accurate 
determination of the Doppler parameter. 
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FIG. 3 The function F{a,ro) describing the equivalent width 
of an absorption feature. Curves are shown for hydrogen 
Lya at three temperatures. Except for large values of tq, 
the curves are independent of a. 



FIG. 4 Distribution of Lya forest column densities with fits 
for different ranges, for absorption systems in the redshift 
range 0.5 < 2 < 1.9. For 12.90 < log^Q iVm < 15.70, /3 = 
1.60 ± 0.03. For 12.90 < l ogmA^m < 17.20, /3 = 1.59 ± 
0.02. From IJanknecht et al\ |2006\ (Fi gure reproduced by 
permission of Astronomy & Astrophysics.) 



For very large values of tq, the Lorentzian radia- 
tion damping wings of the Voigt profile dominate the 
line profile, which in this limit is well approximated by 
0y(a,x) ~ a/iirx^). In this case, the equivalent width is 
given approximately by 



A c V 47r 



1/2 



(26) 



Because of the square-root dependence on column den- 
sity, this limit is referred to as the "square-root" part 
of the curve of growth. While the independence of the 
equivalent width on the Doppler width of the feature 
leaves b poorly determined, the stronger dependence on 
the column density permits a more accurate determina- 
tion of the column density from the equivalent width, 
or from line-profile fitting, than is possible for systems 
on the logarithmic part. The square- root approximation 
Eq. I|26p . however, is accurate only for very large values 
of Tq. It is more than 25% too low for tq < lO**. Better 
than 1% accuracy is achieved for tq > 2 x 10^. 



B. Absorption line properties 

1. Physical properties of absorption line systems 



Systems (LLSs) . The classification is based primarily on 
the physical origin of the features. The Lya forest sys- 
tems, by far the most common, are generally well-fit by 
Doppler line profiles. The much rarer DLAs have suffi- 
ciently high hydrogen column densities that they show 
the radiation damping wings of the Lorentz profile, and 
require the Voigt line profile for accurate fitting. The in- 
termediate LLSs have a sufficiently large column density 
to absorb photons with energies above the photoelectric 
edge, or Lyman limit. The classifications are not strictly 
exclusive. Damped absorbers of course produce Lyman 
Limit Systems. Lower column density Lyman Limit Sys- 
tems will produce a Lya forest feature. For convenience, 
the features, however, are treated as distinct, correspond- 
ing typically to the column density ranges shown in Ta- 
ble |T1 Recently, a subclass of super Lyman Limit Sys- 
tems (sLLSs) has been introduced, correspo nding to sys- 
tems with lO^^ < iVni < 2 x 10^" cm^^ (jWolfe et all 
l2005l ).^^ These systems are convenient for statistical 
studies since their column densities are more easily de- 
termined than their lower column density counterparts 
because the damping wings begin to appear. They are 
sometimes also referred to as sub-damped absorbers, but 
there is an important physical distinction between them 



The absorption features comprising the Lya forest are 
broadly classified into three main types: Lya forest sys- 
tems, Damped Lya Absorbers (DLAs), and Lyman Limit 



^'^ Although the SI system is used throughout this review, column 
densities are expressed with dimensions of cm~^ to be consistent 
with the convention in the literature. 
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TABLE I Summary of absorption line system properties 





Line paramet 


ers 




Physical characteristics 




dJV' _ 
dz 


No{l + z)"' 


Absorber class 


Nm (cm"^) 6" 


(kms*^) 


(m-3) 


T*" (K) 


Size (kpc) 




No 


1 


LyQ forest 




15-60 


0.01 - 1000 


5000 - 50000 


15-1000(?) 


-3.5 - -2 


6.1 


2.47 


LLS 


W^' - lO^^ 


~ 15 


~ 10^ - 10* 


~ 30000 




-3 --2 


0.3 


1.50 


Super LLS 


10*^ - 2 X 10^" 


~ 15 


~ 10* 


~ 10000 




-1 - +0.6 


0.03 


1.50 


DLA 


> 2 X 10^" 


~ 15 


~ 10^; ~ 10* 


~ 100; ~ 10000 


~ 10-20(?) 


-1.5 - -0.8 


~ 0.03 


~ 1.5 



"Approximate ranges. Not well determined for most Lyman Limit Systems and super Lyman Limit Systems. 
''Values not well constrained by direct observations. 

'^Approximate metallicity range, expressed as a logarithmic fraction of solar: [M/H] — logjQ(A//J/) — logjo(M/_ff)0. 
"^For the following H i column density and redshift ranges. For the Lya forest: 13.64 < logjQ A'^ei < 17 and 1.5 < z < 4; 
for Lyman Limit Systems: logj^Q A^hi > 17.2 and 0.32 < z < 4.11. The same evolution rate is adopted for super Lyman 
Limit Systems. The evolution rate of Damped Lya Absorbers over the range 2 < 2 < 4 is consistent with that of 
Lyman Limit Systems, but poorly constrained by observations. 




20 40 60 80 

b (km s-') 



100 



FIG. 5 Di stribution of Ly a forest Doppler parameters for 
QOOOO-26 (|Lu et aZ.I.[l996bl 'l. A lognormal distribution (solid 
line) provides a statistically acceptable fit to the distribution. 



and the general class of DLAs: the hydrogen in the DLAs 
is essentially neutral, while the sLLSs are sufficiently pen- 
etrated by the UV metagalactic background as to be par- 
tially ionized. It is therefore reasonable to distinguish the 
DLAs as an entirely unique class of absorber. 

A great number of surveys of QSO intervening 
absorption systems have been carried out, building 
up a large inventory of their statistical properties. 
Systematic surv e ys of th e Lya forest were co n ducted by 
IWevmann all ^9791) ISargent all (Il 98dl). iBechtold" 



19941) ■ iHu et all (Il995l). IStorrie-Lom bardi et al. 



were conducted by 'Sargent et alU 19891 ) , iLanzetta et al. 
(fr995b) and O'Mcara et al. (200j). Surveys of Dampe d 




(Il986l) 



Lyg Absorbers were co nducted b y iWqlfe et al. .. 

ILanzetta etldl (Il99lt). iLu et all (Il993l). iLu and Wolfd 
(I1994D. ILanzetta et all (ll995bD . IWolfe et all (I1995D 



Ellison et all (|2002| ). I Colbert and Malkan 



Akermari et all 

(I2OO6D and iJorgenson et 



|2005D 
|Jorg( 



IStorrie-Lombardi and Wolfd (I2OOOI) . Ellison et al. 



Srianand et al 
] 



(I2005D . 

( 20061) . A 



iRao et al. 

radio survey of intervenin g abso rption systems was 
carried out bv iGupta et all ( 20071 ). Surveys for galax- 
ies associated with intervening absorption systems 
were conducted for galaxies associate d with diffuse 
Lyg absorption systems a. t z < 1 by Lanzetta et all 
(ll995al).IChen et aZ.I(l2005bl) . lMorris and Jannuzil (|2006l ). 
an d iRvan-Weberl (|2006|), wit h Damped Ly g Abso rbers 
bvlChen an d Lanzettal (I2003D. ICooke et all (|2005l ). and 
ICooke etaL (2006), and with metal a b sorption systems 
at z < 1 by Tumlirison et all ( 2005 ). Prochaska et all 
(|2006D . IStocke e t al! ('2006a') and'Kacprza k et all (|2008D . 
and at z > 1 by Ade lbcrgcr et al. (2005 ). 

The neutral hydrogen column densities of the ab- 
sorbers are measured to range from roughly 10^^ — 
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22 , 



range 

Lower column density systems may exist, 



but are difficult to detec t. An upper cut-off at 3 — 5 x 
10^^ cm~^ is suggested bvlProch agka et all ( 200 5V It was 
early recognized by Tytler, (T9873) that the distribution 
function of the column densities is a near perfect power 



law, rfjV/diVHT P C N~f, with = 1.5 - 1.7 (iHu et 



I1995I : iKim et all . l2002al: iTvtleii Il987al ). A recent de- 
termination for absorption syste ms in the redshift rang e 
0.5 < z < 1.9 is shown in Fig. |5| (jjanknecht et al I I2OO6I) . 
Although t here may be small dev i ations from a perfect 
powe r law (iGiallongo et al. j Il993l : iMeiksin and Madaul . 



I1993I : iPetitiean et a/.l . 199S ). the nearness to a single 



power law over such an enormous dynamic range strongly 
suggests a single formation mechanism. 
iTvtler et all (l2004d) . Surveys of Lyman Limit Systems The measured Doppler velocities b range over about 



19961). iKim et all (l2002a^ ■ IBechtold et all (|2002t ). and 
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10 < 6 < lOOkms" 
concentrated between 15 



with the vast majority 
60 km (lAtwood et al 



Kim 



et al 



1985[ ICarsweh et all Il99ll: I Hu et aLl 1199a - 
19971: iKi rkman a nd Tvtleil . Il997a': 'L u et aH 
Ranch et al.^ .1992) . A typical distribution is shown in 



1996b : 



Fig. El using the line list for QOOOO-26 from E u et al\ 
(|1996bl ). Temperatures may in principle be inferred from 
Eq. PH)) . but doing so is hampered by two difhculties: 
1. the systems may be broadened by a kinematic com- 
ponent and 2. the absorption features may be blends 
of more than a single system. Evidence for kinematic 
broadening is found when metal features are also de- 
tected (see below). In general, there is no unique fit 
to an absorption feature, particularly in the presence of 
blending: several statisti c ally ac ceptable fits are possi- 
ble (.Kirkman and Tvtleil . [T997a[) . and these will change 
as the signal-to-nois e ratio or spectral resolution changes 
(jRauch et all Il993l ). Indeed, Eq. ((H]) shows that each 
absorption feature itself may be regarded as the blend- 
ing of an infinity of smaller features. It is only because 
of dumpiness of the IGM that the features may be local- 
ized (Eq. [16|). yet internal structure is still visible when 
multiple narrower metal features are detected in a single 
Lyo! system. 

Broadening is also expected from the line finding and 
fitting procedure. The systematic influence of procedures 
used to locate and fit the absorption lines on the result- 
ing distributions of Doppler parameters has received lim- 
ited attention. The potential usefulness of the Doppler 
parameter distribution for extracting physical informa- 
tion about the IGM, such as its temperature distribu- 
tion, merits further study of the effect of fitting algo- 
rithms on the derived distribution. Artificially large 
^-values, for example, have been reproduced through 
Monte Carlo simulated spectra with much narrower gaus- 
sian d i stributions for th e 6-values (Kirk man and Ty tler. 
Il997al: ILu et all Il996b( ). The Monte Carlo models as- 
sume Poisson distributed line centers. Allowing for clus- 
tering of the line centers (see below) may lead to even 
broader distributions. 

The measured distr ibutions may be fit by a lognor - 
mal function (Mciksin^iL^, l200ll: IZhang et all Il997t ). 
As an illustration, the best fitting lognormal function 
/(^) = exp(-0.5^2)/y2^^ .^yjjgj.g ^ ^ _^ 2.13 ln(6), 

provides a statistically ac ceptable descript ion of the h 
distribution of QOOOO-26 (iLu et all Il996bl ) for systems 
lying between the restframe Ly/3 and Lya wavelengths 
of the QSO, excluding a small region possibly influenced 
by the proximity effect (see § IIV.C.l|l . The fit distribu- 
tion is shown in Fig. [51 The KS test yields a probabil- 
ity for obtaining as large a difference as found between 
the predicted and measured cumulative distributions of 
Pks — 0.361. A lognormal distribution will result when 
the error on a quantity is proportional to the magnitude 
of the quantity. For example, consider an absorption sys- 
tem with an intrinsic Doppler parameter b. If the itera- 
tions of the nonlinear fitting procedure of the remaining 
features in the spectrum each perturb the first feature 



by an amount A6 oc &, then A In & will change by a ran- 
domly distributed amount. Fitting several features will 
then result in a sum of random changes to A In b. By 
the central limit theorem, a normal distribution for In b 
will result. A compar ison between the measured 6- values 
and their errors from lLu et aZI (|l996bh shows a positive 
correlation. The Spearman rank-order correlation coeffi- 
cient is Ts = 0.31, with p{> Ts) = 9.3 x 10~^. A similar 
result is found for HS 1946-f 7658 (jKirkman and Tvtleil . 
Il997all . for which a correlation between the Doppler pa- 
rameters and their errors is found with = 0.28 and 
p{> Ts) = 7.4 X 10~^. A lognormal distribution again 
provides an acceptable fit to the Doppler parameter dis- 
tribution. The results for this and several other distribu- 
tions measured in high resolution Keck spectra are shown 
in Table [TTl In all lognormal distribution provides 

a good fit. Although not conclusive, these results suggest 
that the measured Doppler parameters may in part be in- 
duced by a lognormal-generating stochastic process. It is 
noted that the process need not arise entirely from the 
line-fitting, but could also originate from the physical 
processes that gave rise to the structures that produce 
discrete absorption features. 

Uncertainty in the origin of the Doppler parameter dis- 
tribution leaves the physical interpretation of the Doppler 
parameters with some ambiguity. Although their rela- 
tion to the gas temperature is not straightforward, the 
Doppler parameters may usually be used legitimately to 
set upper limits on the gas temperature. (Even an up- 
per limit is not always guaranteed, since noise affecting 
a weak line may leave it too narrow, or even produce 
an artificial narrow line.) For pure Doppler broaden- 
ing, the range 20 — 60 km s^^ corresponds to tempera- 
tures of 2.4 X 10^ — 3.8 X 10"' K, the range expected from 
photoionization a t the m oderate ovcrdensities expected 
for the absorbers (jllui an d Gncdin, 1997; Meiksin, 1994). 
Cooler temperatures are possible, however, particularly if 
the gas has been expanding sufficiently fast for adiabatic 
cooling to be appreciable. 

The element of indeterminacy in the measurements of 
the line parameters has made it difficult to search for evo- 
lution in the distribution of Doppler parameters. Evolu- 
tion is especially interesting in the context of late He ii 
reionization, which may have occurred at z ~ 3 — 3.5. 
Evolution of the Doppler parameter distribution, how- 
ever, could also result if aborption systems with different 
physical characteristics dominate the detected and fit ab- 
sorption lines at different epochs. Numerical simulations 
suggest the systems that give rise to a given H i column 
density range shift to structures of different gas densities, 
and therefore different gas temperatures, and different 
sizes as the Universe expands and as the intensity of the 
photoionizing UV background evolves. The difficulty of 
interpreting any change in the Doppler parameter distri- 
bution in terms of sudden heating is exacerbated by an 
increase in line blending with increasing redshift. Using 
an analysis based on the lower cut-off in th e Doppler pa- 
rameter distribution. ISchave et ai\ ( 2000bl ) report evolu- 
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TABLE II Lognormal fits to the Do ppler parameter distr i bution s measured for tfie Lya forest seen in tiie spectra of QOOOO-26 



Lu et al. 



Hu et al 



HS 1946+7658 l|Kirkman and Tvtleil , 11997^ ). and Q0014+813, Q0302-003, Q0636+680 and Q0956+122 
Tiie fits are represented in ttie form = exp(— 0.5^^)/\/27r, wiiere ^ = 7 + 5\n{h). 



QSO name 


z range 


no. lines 


(6) (kms^^) 


7 


5 


Pks 


QOOOO-26 


3.42 <z< 3.98 


334 


32.8 


-7.15 


2.13 


0.36 


Q0014+813 


2.70 <z< 3.20 


262 


33.9 


-8.86 


2.58 


0.99 


Q0302-003 


2.62 <z< 3.11 


266 


33.9 


-8.75 


2.56 


0.97 


Q0636+680 


2.54 <z< 3.03 


312 


29.3 


-7.24 


2.23 


0.30 


Q0956+122 


2.62 <z< 3.11 


256 


31.6 


-8.05 


2.42 


0.53 


HS 1946+7658 


2.42 < z < 2.98 


399 


32.5 


-7.18 


2.16 


0.50 



tion in the inferred IGM gas temperature over the range 
2.0 < z < 4.5, peaking at z ~ 3, consistent with the on- 
set of He II reionization at this redshift. Si milar results 
are ob t ained frorn a sepa rate analysis by iRicotti et al\ 
(|2000l ). iKim et al\ (|2002bl) also suggest evolution in the 
gas temperature based on an increase in the lower cut- 
off of the Doppler parameter distribution from z — 3.8 to 
z = 3.3. But the value s for the cut-o f f they derive are dis- 
crepant with those of ISchave et al! ('2000b ), which casts 
uncertainty on either interpretation. K im et al\ ( 2002b| ) 



show there is considerable sample variance in the cut- 
off at the same redshift, and attribute the discrepancy 
withfSchavc et al. (2000t) to this variance. On the other 
hand, Le hner et ali (|2007l ) find no significant evolution in 
the 6-distribution over t he redshift r ange 1 .5 < ^ < 3.6 
using the same sample as iKim et ~ai. l (|2002aD . This result 
is confirmed by the nearly identical parameters obtained 
for the best fit lognormal distributions bet ween the high- 
est an d lowest redshift samples in Table HIl iLehner et al\ 
(|2007D do argue for the presence of a second broader pop- 
ulation at z^O.4. fjanknecht et al. (2006) find no evi- 
dence for evolution in the Doppler parameters over the 
redshift range < z < 2. 
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FIG. 6 Evolution of the number density of Lya systems with 
13.10 < logjQ A^'hi < 14.00, including best-fits to the solid 
and open circles over the ranges indicated, as represented by 
the solid and dashed lines, respectively, taking dN/dz cx (1 + 
z)^. The d otted line extrap o lates the lower redshift fit to 
z = 0. From lJanknecht et al\ (|2006t ). (Fi gure reproduced by 
permission of Astronomy & Astrophysics.) 



Damped Lya Absorbers show a more complex physi- 
cal structure. Measurements of the 21cm hyperfine ab- 
sorption signature produce a range of spin temperatures, 
from values as low as Ts — 200 K to lower limits of 
Ts ^ 9 000 K, with the lower v alues occurring only at 
z<3. iKanekar and Chengalud (2003) provide a tabula- 
tion from the literature. The presence of high ioniza- 
tion metal species show there is a second warmer com- 
ponent pr esent as well, with T ~ 10"' K (see below). 
IWolfe efa l. (2005) suggest a two-component interstellar 
medium consisting of a warm neutral medium (WNM) at 
T ~ 8000 K in pressure equilibrium with a cold neutral 
medium (CNM) at T ~ 100 K. 



2. Evolution in the number density of the Lya forest 

The number of absorption systems per unit redshift 
increases with redshift. Some evolution is expected as 
a result of the expansion of the Universe. For a proper 
number density na{z) of absorption systems at redshift z, 
with proper cross section aa{z), the expected number of 
absorbers per unit proper length is dN /dip — na{z)aa{z). 
The proper line element is related to redshift according 
to dlp/dz = c/[7?(z)(l + z)], where H{z) = HqE{z) is the 
Hubble parameter (see ij lH.AI above). For a flat Universe 
= 0) and standard cosmological parameters, E{z) ~ 
0.55(1 + z)3/2[l + 2.3/(1 + z)3]i/2. The evolution in the 
number density is then given by 



A summary of the absorption line properties, as well 
as of the physical characteristics discussed below, is pro- 
vided in Table HI 



^ ~ (2100Mpc)na,c(z)(Ta(2)(l+2)'/' 
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FIG. 7 Evolution of the number density of Lya systems with 
logj^Q A'hi > 13.64, including best-fits to the solid and open 
circles over the ranges indicated, as represented by the solid 
and dashed lines, taking dN/dz oc . The additional low 

redshift dashed line is the best fit to the open triangles. The 
open square, left-pointing solid triangle, and stars correspond 
to sys tems with log^Q A''hi > 14.00. From iJanknecht et al\ 
12006"). (Figure reproduced by permission of Astronomy & 
Astrophysics.) 



Lyman Limit Systems, with dN/dz inc reasing by about 
a fac tor of 2 from z = 2.5 to z = 4 ( Prochaska et all 
I2005D . 



As will be discussed in § I VII below, numerical simu- 
lations show that the structure of the IGM evolves, al- 
though in the comoving frame the structure is remarkably 
stable. Because dN/dz is fit for a fixed H i column den- 
sity range, the diminishing gas density towards decreas- 
ing redshifts will result in fewer systems satisfying the 
column density threshold, so that dN/dz will decrease 
towards decreasing redshift. Evolution in the ionizing 
background will also affect the number of absorption sys- 
tems lying above the column density threshold, and this 
is a major factor in the evolution of dN/dz. The slow- 
down at z < 1.5 in the evolution is in fact attributed 
predominantly to a reduction in the intensity of the ion- 
izing background: as the ionizing rate decreases, fewer 
systems will slip below the column density threshold than 
under pure density evolution. As a result, the decrease in 
dN/dz towards lower z slows down. The difference in the 
rate of evolution between low c olumn density and high 
column density systems found bv IJanknecht et all ( 20061 ) 
(compare Figs. [5] and [7]), shows that structural evolution 
in the IGM must also play a role. 



where na^c{z) — na{z){l + z)^^ is the comoving number 
density of systems. For a constant comoving number den- 
sity and fixed proper cross-section, only moderate evolu- 
tion is expected, dN/dz oc (1 -I- z)^/^. At low redshifts, 
this gives a reasonable description of th e evolution. Us- 
ing Hubble Space Telescope observations. IWevmann et all 
(fl998) find for all system with equivalent widths above 
0.24A, dN/dz = (32.7±4.2)(1 -hz)0-26±0-22 foj. ^ < ^ 5_ 

Subsequent higher spectral resolution surveys using HST 
have extended the statistics to sy stems with equiva- 
lent w idths below O.lA at z < 0.1 (iPenton et al. 2000, 
|2004| ). The large survey of lOanforth and ShuU (2008. ) 
obtains dN/dz = 1291:5 down to O.OBOA at z < 0.4 with 
(z) ~ 0.14. 

For z > 1.5, however, the number density evol ves con- 
siderab ly faster than theconstant comoving rate (|Tvtlerj . 
Il987b[ ). iKim et al\ (^2002a^ obtain dN/dz = 6.1(1 + 
^■j2.47±o.i8 systems with column densities in the range 
13.64 < logio -^Hi < 17. This corresponds to consider- 
able evolution in the product na,c{z)(Ja{z). A smooth 
transition is found over the redshift range 0.5 < z < 1.9 
by [janknecht et all ( 2006[ ). with a dependence of the 
evolution rate on column density. High column density 
systems corresponding to Lyman Limit Systems evolve 
some what more slowly than lowe r column density sys- 
tems. IStengler-Larrea et all ( 19951 ) find for systems with 
TL > 1, dN/dz = 0.25;°;J°(1 + z)i-50±0-39 over the 
redshift range 0.32 < z < 4.11. B ased on a larger, 
higher redshift sample, iPeroux et all ([2003.) find some- 
what more rapid evolution, with dN/dz = 0.07:|:o;i3(l + 

2 4r; + 7B 

z) ■ -0 65 over the redshift range 2.4 < z < 5. Damped 
Lya Absorbers show evolution comparable to that of the 



3. Characteristic sizes and spatial correlations 

Multiple QSOs in proximity to each other on the sky 
provide an invaluable means for estimating the physical 
sizes of the regions in the IGM that give rise to the ab- 
sorption features. Especially useful have been multiply 
imaged lensed quasars. Using high resolution spectra of 
a pair of images of a lensed QSO with a separation of 
3 arcsecs, ISmette et all (1995) detected a large number 
of absorbers coincident between the two lines of sight. All 
the clearly detected coincident systems show comparable 
equivalent widths between the lines of sight, suggesting 
the same physical structures are being probed. These 
authors set 2a lower limits of about 70 kpc for weak ab- 
sorption systems (log^o -^hi < 14), and of about 15 kpc 
for strong systems (log^g A^hi > 15), assuming the clouds 
are spherical. For disk-like systems, the lower limits are 
about twice as large. 

Other effort s resu lt ed in s imilar lower limits 
(iBechtold et all . Il994l : iMcGilll . Il990l : ISmette et all 
Il992 ). Multiple lines of sight with wider separations 
of about 1 arcminute sug g est st i ll larger sizes o f 0.3- 
1 Mpc (ICrotts and FaM . 'l998': 'Pinsha w et~ai\ . \l99A 



iD'Odorico et all 11998k Pctitjean et all Il998l). and 



possi bl y as large as 4 Mp c (comoving) (jWilliger et all 
I2OOOI ). ISmette et all ( 19951 ) report a likely size of about 
10-30 kpc for a single Damped Lya Absorber, compa- 
rable to the lower limit obtained by direct H i imaging 
measurements of a separate system by iBriggs et all 
11989). 

A difficulty with interpreting the larger Megaparsec- 
scale sizes for the Lya forest systems is that the equiva- 
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lent widths of systems with common redshifts along the 
neighboring lines of sight do not always match well over 
these scales. An alternative interpretation of the co- 
inciding redshifts is that they indicate underlying spa- 
tial correlations between the absorption systems rather 
than their spatial coherence. Measurements of the auto- 
correlation function between absorption systems along 
single lines of sight sugg est correlations on sca l es of a few 
hund r ed km s~^ scales dCristiani et al. u et ali . 

119951 : iMeiksin and Bouchetl . Il995[) . An anti-correlation 
at separa t ions o f 500— 1000 kms~ (ICristiani et all Il995l : 
IHu et all . Il995 : Meiksin and Bouchetl . Il995 ) has been 
disputed by Kim et aJT l 19971 ) on the basis of the limited 
lines of sight analysed, but the issue appears not fully 
resolved. In the presence of correlations, the error on the 
correlation function depends on the 3 and 4-pt functions, 
which are not well-measured. A more recent analysis 
at z < 2 detects only weak positive correlations on the 
few hundred km s~^ sc ale, and no clear anti-correlation 
( Janknecht et a/.l . |2006[ ) . Transverse positive spatial cor- 
relations over comparable scales at (z) = 2.14 were ob- 
tained bvlCrotts and FanS (119981 ) and at 0.4 < z < 0.9 by 
iPetrv et all ( 2006 ). Coherence in absorption lines along 
lines of sight separat ed by ^ 2h~^ Mpc at 2.6 < z < 3.8 
has been detected bv lCasev et al\ ( 20081 ). 

A characteristic size of 70 kpc for systems with iVni — 
10"'^^ cm~^ suggests a large fractio n of the baryons are 



conta ined within the Lya forest ( Meiksin and Madaul . 
I1993D . Taking the size to correspond to a characteris- 
tic line-of-sight thickness, the corresponding H i num- 
ber density is tt-hi — 5 x lO^^m"'^. Ionization equi- 
librium in a UV background producing an ionization 
rate of Fhi — rHi,-i2lO^^^ H i atoms per second 
(§ I IV. A. 21 below), gives a neutral fraction of xhi 

(1.2nHiaA/rHi)^/' ^ 1.5 x IQ-^N^^l^T^ 



where aA ~ 4 x 10" 



'19 , 



-0.75 



0.37p-l/2 

1,-12' 

is the radiative re- 



combination rate for gas at a temperature T4 = T/10 K 
and iVHi,i4 A'hi/10^^ cm~^. The ratio of the baryonic 
mass density of the Lya forest to the critical Einstein- 
deSitter density Pcrit{z) = Pcrit{0)H{z)'^ / H§ is 



1.4mH 

Pcrit 



aw Nmfdlp^ ^ 



dNmdz xm \ dz 



\Amy 



8nG <9W Nm 



-(l + z) / dN^ 



'3ci/(z)^ ' 'J ™'dNiiidz xm 

~ 3.0 X io-^Noh-^n;^^/^T°-^'^rK^^_^2 

X (l + z)^-i/2ijY^Wax\ 

V A^HI.min / 

^ OMT°''r'£_,,, (28) 

for dN/dz = A^o(l + 2)'^ with A^o = 6.1 and 7 = 2.47, 
and dN/dNm oc N^l'^ with iVm.min = lO^^ '^'* cm-^ and 
A^Hi.max — lO^^cm"^, evaluated at z = 3 with — 0.3 
and h — 0.7. This is comparable to estimates for the 
baryon density of the Universe (see § III. B. 41 below) . The 
result is not conclusive, however. Since the size estimates 



are in the transverse direction, a much smaller value for 
r^Lya is possible if the structures are flattened, down by 
the aspect ratio of the thickn ess to breadth of the systems 
( Rauch and Haehneltl . [19951) . 



4. Deuterium absorption systems 

High column density Lya absorption systems offer the 
opportunity to measure the primordial deuterium abun- 
dance D/H, a sensitive indicator of the cosmic baryon 
density, or equivalently the nucleon to photon ratio 77, 
according to Big Bang nucleosynthesis. To date, a rel- 
atively limited number of accurate determinations of i] 
from the deuterium abundance have been possible. A 
line of sight must meet several c riteria to be us eful for 
obtaining a clean measurement ( Kirkman et all [2003): 
the hydrogen column density must be sufficiently large 
that deuterium is detectable; the velocity width of the 
system must be sufficiently narrow so as not to encroach 
on the 82 km s~^ offset deuterium feature; the structure 
of the system must not be so complicated that a weak 
interloping H i feature could be mistaken for deuterium; 
and the background QSO must be sufficiently bright to 
obtain a high signal-to-noise ratio measurement. These 
combined criteria whittle down the number of accept- 
able QSO lines-of-sight to about 1% of those available at 
z = 3. 

Despite the difficulty, several accurate deter- 
minations have been made (iBurles and Tvtler . 



| l998al[bl: ICrighton et all \2004 [Kirkman al 
Tevs hakov et aL[. 12002^ [O'Meara et all [2006 



2003; 



2001 



iPettini and Bowenl l200lh . The curren t best estimate 
for D /H is logioD/H = -4.55 ± 0.04 (jO'Meara et all 
[2006[ ). Using standard Big Bang nucleosynthesis mod- 
els, this translates to a baryon density parameter of 
rith"^ = 0.0213 ± 0.0013 ± 0.0004, where the Icr errors 
refer to the error on the D/H measurement and the 
uncertainties in the nuclear reaction rates, respectively. 
This corresponds to a nucleon to photon ratio of 
T] = (5.8 ± 0.3) X 10"^". This compares favorably to the 
determination from the 5-year WMAP measurements of 
the CMB fluctuations, combined with Baryonic Acoustic 
Oscillations and Type la supernovae constraints, of 
rj ~ (6.2 ± 0.2) X lO^i" (iKomatsu et all [looi). 



5. Helium absorption systems 

The UV spectral capability of the Hubble Space Tele- 
scope made it possible for the first time to measure 
intergalactic He 11 Lya A304. Using the spectrum of 
Q0302-003 a t Zem ^ 3.286 observe d by the Faint Ob- 
ject Camera, [Jakobsen et al. I ([1994D found thcII > 1-7 
(90% confidence) at z = 3.2, re vised to a 95% confidence 
interval of 1.5 < THeii < 3.0 by Hogan et a/1 ( 1997t) and 
to thoII > 4.8 bv [Heap et al\ ( 2000t) after more detailed 
modeling of the absorption trough. Because the same 
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sources which ionize the H i should also ionize He i, vir- 
tually all the intergalactic helium will be in the form of 
either He ii or He iii, depending on whether or not there 
were sources hard enough to fully ionize helium. 

Remarkably, using the spectrum of HS 1700+6416 
(zom — 2.7) observed by the Hopk ins Ultraviolet Tele- 
scope {HUT), Davidsen_eL2l] (|l996l ) measured the more 
moderate optical depth of THeii — 1-00 ± 0.07 over the 
somewhat lower redshift range z = 2.2 — 2.6. This led to 
the speculation that He ii was ionized around z ~ 3. Sub- 
sequent measurements of the spectrum of HE 2347—4342 
(zem — 2.885) have shown the He ii Lya optical depth to 
be very patchy at z ^ 3, as if a sec ond reionization epoc h 
were being observed, that of He ii ( Reimers et aZ.I . [l997[ ) . 

Ob servations of PKS 19 35-692 (zcm =i 3.18) using 
HST (lAnderson et aLl . ll999D and of HE 2347-4342 using 
the Far Ultraviolet Syectroscopic Exp l orer have detected 
a He II Lyg forest (iKriss et al\ . l200ll : IShull et all \2004 
IZheng et all . l2004b ). permitting more detailed assess- 
ments of the fluctuations in the He ii to H i column den- 
sities. Measurements of the He ii Lya forest also provide 
an opportunity to probe the underlying I GM velocity 
field. U sing combined FUSE and VLT data. lZheng et all 
(l2004bD were able to identify isolated absorption systems 
in underdense regions with common H i Lya and He ii 
Lya features. A comparison between the fit Doppler pa- 
rameters gives the ratio 6hoIi/&hi = 0.95 ± 0.12. Accord- 
ing to Eq. (fTTj) . this indicates a velocity field in under- 
dense regions with kinematic motions that dominate the 
thermal. 

Subsequent measurements continue to confirm the 
patchy nature of the He ii Lya optical depth, includ - 
ing HST observations of HE 2347-4342 (ISmette et al. 
20021) SDSS J2346-0016 (zcm 3.50) (IZheng et al. 
2004aD . and Q1157-F3143 (zom =i 3.01 ) (| Reimers et al 
20051) ■ and a FUSE observation of HS 1700-^6416 



( Fechner et "aLl. l2006D . As some of these references point 
out, and as will be discussed in § IIV.C.2I below, the in- 
terpretation in terms of reionization is not unambiguous: 
the patchiness may also arise from fluctuations in the 
radiation field. 



6. Metal absorption systems 

The earliest identifications of intervening absorp- 
tion systems were for elements he a vier than he- 
lium, so called metals ( Bahcall et all . Il968t) . Com- 
mon are carbon, nitrogen, silicon and iron, but to 
date the list extends much furthe r, including oxy- 



gen, magnes i um, ne on, and sulfur (ICowie et al. j_ 1995 ; 



Ellison et al!, 'l999a''b'; 'Koehler et al!. 'l99 6t iLu et al. 



1998a: Mcvcr and York, 1987- Richter eioL, 2004 ; 



Sargent et all. Il980l: ISavage all . 20051: Schave et al. 



2000at ISimcoe et all l2002l 12004 ISongaila and Cowie , 
19961) . The metals provide invaluable insight into the 
structure of the IGM in several different ways. The 
widths of metal absorption features allow direct esti- 



mates of the temperature of the IGM and its small-scale 
velocity struc ture. The n arrow widths of C iv features 
were used bv lYork et al. I (Em to demonstrate the ab- 
sorbers had the characteristic temperatures of photoion- 
ized gas rather than coUisionally ionized. The metal sys- 
tems probe the impact star formation has had on the 
IGM, as the metals were most likely transported into the 
IGM through galactic winds, or introduced in situ by lo- 
cal small-scale regions forming first generation. Popula- 
tion III, stars. In the nearby Universe, they indicate the 
presence of shocked gas, as may accompany the forma- 
tion of galaxy clusters. As such, the metals in principle 
document the history of cosmic structure and star for- 
mation. Ratios of the metal column densities may be 
used to constrain the spectral shape of the metagalactic 
ionizing UV background, which in turn puts constraints 
on the possible sources of the background and their rel- 
ative contributions. Lastly, metal absorption lines offer 
a unique opportunity to search for variability in the con- 
stants of nature. 

Tables of the resonance lines most easily detectable 
in intervening absorpt ion systems are provided b y 
iMorton and YorkI (|l989l ) and lVogel and Reimerd (|l995l ). 
Searches for intervening metal absorbers have been 
facilitated by the exploitation of atomic line doublets for 
some of the stronger species. The most common doublets 
used are Mg ii AA2796, 2803 A and C iv AA1548, 1551 A. 
These are among the strongest lines detected in the 
Lya forest. They also have the important advantage of 
wavelengths redward of Lya, placing them outside the 
Lya forest and avoiding this potential source of confu- 
sion. Searches for Mg ii absorption systems are by far 
the most common because of their established connection 
to galaxies. Mg ii ab s orber survey s were carried out by 
iLanzetta erall (119871) .iTvtler et a/I (ll987l) . [Sareent et al! 



1988bD . ISteidel and SargentI (119921). \Boisse et al. 



1992D. lAldcroft aZ.I (I1993D. iLe Brun eit a/.l (Il993t) 



ChurchiU e 7aZ.I (Il999l) . [ Nestor ei aZ.I j2005l). lYork et all 



( 20061). iNestor et g^. 



(I2006D . 

Sur- 



.2006). .Lvnch et all 

iFox et aZ.I (|2006( ) and iNaravanan et all ((2003)- 
veys of C iv absorpt ion s ystems were c onduc ted by 
ISargent et dl (|l988a[) and iMisawa et all (|2002l ). Ab- 
sorption features by O vi AA1032, 1038A are more 
difficult to identify since they lie within the Lya forest, 
and so are not easily distinguished from the forest. 
Searches have been successful, however, and have more- 
over discovered a second IGM environment. While at 
high redshifts the O vi sy stems probe the component of 
the IGM co ntaining C iv (iBurles and Tvtleiill996ll2002l: 
iKirkmaiTand T vtleil ll997bf ). at low redshifts they are as- 
sociated with a warm-hot int ergalactic medium (WHIM) 
jDanforth and ShuU', '2005', '2008'; 'Richter et al!, '20o4 
Savage et al., 2002; Tripp and Savage, 2000; Tripp et qlj, 
20081 ). that is most likely coUisionally ionized rather 
than photoionized. A search for O vi absorption sys- 
tems as an i ndicator of a WHIM at high redshift was 
conducted bv lSimcoe et Ml (|2002D . Efforts are currently 
underway to detect an intergalactic x-ray absorption 
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line forest using the XMM-New t on and Chandr a x-ray 
sateUites (iKaastra et oZI . 120061: iNicastro efal\ . 120051: 



iRasmussen et Detections of low redshift 

O V III have been r e porte d by iFang et all (|2007l |2002| ) 
and lWilhams et all pOOTi ). 



Attempts to infer the origin of the metals, whether 
distributed by galactic outflows or introduced locally 
by Population III stars, have led to various efforts 
to estimate the range in ambient metallicities and 



to determine their spatial distribution (lAracil et al. 


2004; Barlow and Tvtled, 


1998 


; iBergeron et al 


,. ,1994 


Bergeron and Stasinska, 


1986 


: Cowie and Songaila 


1998 


; Cowie et all. 1995; 


Ellison et all 2000; Fieri et al. 


2006 


; IScannapieco et al. 




2006bl: Simcoe et al. 


, 2ooe 


Songaila, 200 iL l2005l 2006 


; Songaila and Cowie 


. Il996h 



Updated solar abundances useful for modeling the 
abundan ces of the metal abs o rption systems are pro - 



vided by Grevesse and Sauvall (ll998^ ■ iHolwege: 
lAllende Prieto et all (|200ll . l2002h and lAsplun ' 



2001), 



20031). 



Early attempts to tease out weak metal features were 
based either on a spectral shift-and-stack approach to 
generate a high s ignal-to- noise ratio composite C iv line 
(iLu et all . IT998d) or throug h a pixel-by-pixel sta t istical 
analysis of the spectra ( Cowie and Songaflal . Il998l : 
ISongaila and Cowig Il996). A compar i son of these two 
techniques is provided by lEllison et all ( 2000t ). 

An alternative is to perform long integrations of bright 
QSOs to obtain high resolution, high signal-to-noise ra- 
tio spectra. With signal-t o-noise ratios of 200 -300 per 
spectral resolution element. lSimcoe "et "^ (|200l detected 
metal lines in 70% of the Lya forest systems. Ab- 
sorption systems with H i column densities as low as 
iVni > lO^^'cm"^ show C iv features, while systems with 
A^Hi > 4 X 10^'^cm^^ show O vi features. Although the 
inferred metallicities are model dependent (fixed in part 
by the assumed spectrum and intensity of the ionizing 
background), the measurements suggest metallicities as 
low as 3 X 10""* that of the Sun. No evidence for a metal- 
licity floor could be discovered. The metallicity distribu- 
tion inferred is consistent with a lognormal distribution 
with a mean of 0.006 solar. 

The mechanism and degree of metal mixing in the 
diffuse IGM are unclear. While the mixing of metals 
in stellar interiors by convection and diffusion may pro- 
duce a homogeneous composition, preserved in supervo- 
vae ejecta and winds, it is far from clear that the mixing 
of the ejecta will result in uniform metallicity. The mix- 
ing process is still poorly understood in the Interstellar 
Medium, and much more so in the IGM. Just as stirring 
a drop of milk in a cup of coffee will result in the mixing 
and intermingling of the two fluids, so may the stirring 
produced by dynamical instabilities such as the Kelvin- 
Helmholtz and Rayleigh- Taylor mix metal enriched stel- 
lar ejecta with the surrounding primordial gas. Insuffi- 
cient mixing, however, rather than a uniform mixture, 
will instead result in striations or patches of high metal- 
licity gas entrained withi n the primordia. 1 gas of essen- 
tially zero metallicity. iSchave et all ( 2007t) find evidence 
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FIG. 8 Portion of the Lya forest in the spectrum of 
Q1422 -I- 231 at 2 = 3.6 (lower panel) with the correspond- 
ing Ly/3 systems (middle panel) and A1548A systems of the 
C IV AA1548, 1551 A doublet (upper panel) shifted to coin- 
cide. Absorption systems saturated in both Lya and Ly/3 are 
indicated by the long vertical dotted lines. The short dotted 
lines show where the correspond ing C iv A155lA line shoul d 
occur for these systems. From ISongaila and Cowid l|l996t ). 
(Figure reproduced by permission of the A AS.) 



for just such inadequate mixing for diffuse absorption sys- 
tems (with a median upper limit of 13.3 on logj^g-^Hi), 
at z ~ 2.3. By combining C iv measurements with up- 
per limits on H I, C iii. Si ii. Si iii, N v and O vi, they 
place robust constraints on the metallicities and physical 
properties of the C iv absorption systems, with a median 
lower limit on the metallicities of [C/H] > —0.42 and a 
median upper limit on the absorber sizes of i?^ 1.5 kpc. 
They obtain typical median lower and upper limits on the 
hydrogen densities of nn > 100 m~'^ and tih < 1000 m~'^, 
a range that suggests a cloud size of about 100 pc. The 
clouds, however, would likely be transient, either quickly 
dispersing if freely expanding or shorn apart by dynami- 
cal instabilities if moving relativ e to a confining medi um, 
on a timescale of about 10^ yrs. ISchave et all ( 2007t ) ar- 
rive at a picture in which new clouds are continuously 
created through dynamical or thermal instabilities. In 
this picture, nearly all the metals in the IGM could be 
processed through just such a cloud phase. The dispersed 
clouds would retain their coherence as small patches, but 
of too low column density to be detectable. The resulting 
IGM metallicity would persist as patchy, not mixed. 

Matching the metal absorption lines to the cor- 
responding H I Lya features frequently reveals mul- 
ticomponent substructure obscured by the width of 
the Lya feature. An example from QSO 1422 -|- 
231 is shown in Fig. [51 The C iv systems show ev- 
idence for spatial clustering on scales of up to ^ 
300kms~^ (iD'O dorico et aZ .I. Il998al: iPichon et aLl . 120031: 
iRauch et all . 11996: Songa ila and Cowid . Il996t) . This 
is consistent with a 2a lower limit of ~ 30 kpc for 
the clustering of C iv and Mg ii systems obtained by 
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ISmette et al. I (Il995h based on a statistical analysis of co- 
incident metal systems in neighboring lines of sight to 
a lensed QSO. An analysis of a triply imaged quasar at 
z = 3.9 bv lEllison etai\ (|2004[ ) shows that high ionization 
systems like C iv and Si iv absorbers retain coherence 
over scales of ^ 100 kpc. By contrast, low ionization sys- 
tems like Mg II and Fe ii absorbers show coherence only 
over much smaller scales. Assuming a spherical geome- 
try, they are estimated to have sizes of 1.5-4.4 kp c (95% 
c.i.). Using similar measurements, iRauch et al\ ( 200ll ) 
find C IV systems at 1.5^ 3 are coherent over scales 
of about 300 pc, with column densities along neighbor- 
ing lines of sight differing by 50% for separations of about 
1 kpc. Detectable levels of C iv extend over scales of sev- 
eral kiloparsecs. 

Damped Lya Absorbers show a wide range of met- 
als with both low and high ionization states, in- 
cluding boron, carbon, nitrogen, oxygen, magnesium, 
aluminum, silicon, phosphorus, sulfur, chlorine, ar- 
gon, titanium, chromium, manganese, iron, cobalt, 
nicke l, copper, zinc, germanium, a r senic an d kryp- 
ton (iDessauges-Zavadskv et al\ . 120061 : iLu et al\ . Il996at 
IWolfe et aLl . 12005 ). The metallicities evolve from about 
0.03 solar a,t z A to 0.15 so lar at z ~ 1, with few above 
0.3 solar (jWolfe et al\ . l2005h . No DLA is found with a 
metallicity smaller than 0.0025 solar, distinguishing them 
as a population from the Lya forest. By contrast, su- 
per Lyman Limit Systems are found to have metallicities 
ranging from 0.1 solar to super-solar abundances with as 
much as 5 times solar in zinc ( Kulkarni e t al, 2007). As 
both these classes of systems are likely related to galax- 
ies, further discussion of them is deferred to § I Villi 

The presence of the metal lines in absorbers of all H i 
column densities in principle offers an opportunity to ob- 
tain a nearly model-free measurement of the gas temper- 
ature. Allowing for kinetic broadening of the Doppler 
parameters following Eq. (jlip , the measurements of the 
Doppler parameters from features corresponding to two 
or more elements permits the thermal and kinematic 
broadening components to be separated. (In principle, 
instrumental broadening could be incorporated into the 
kinematic term.) Given total Doppler parameter mea- 
surements bi and bj for two species i and j of masses rrii 
and rjij , the temperature may be solved for as 



T 



2kB 



(29) 



Applying this method to C iv and Si iv features de- 
tected in three high redshift QSOs, and assuming that the 
C IV and Si iv features trace the same syst ems with the 
same temperature and velocity structure, iRauch et all 
(Il996h obtain a typical temperature of T « 38000 K 
and typical kinetic motions of fckin ~ 6 km s~^ for a se- 
lect subsample of especially well-fit features. The large 
C IV column densities of the systems (typically exceed- 
ing 10^^ cm~^, M. Ranch, personal communication), sug- 
gest large H i colu mn densities of iVni > 10^^ cm~^ 
(jSimcoe a/] . l2004) . 



A few other Lya systems selected for quiescent veloc- 
ity fields as part of a program of deuterium abundance 
measurements yield somewhat cooler temperatures. The 
Doppler param eters of C iv and Si iv measured by 
iKirkman et all (,2000 ) in a system with logj^Q A^hi — 16.7 
yields a temperature of T w 30000 - 32000 K. For a 
Ly man Limit System with logj^Q A^hi = 17.1 at z = 0.7 
(jTvtler et all [l9991 the Doppler parameter measured 
for Mg II compared with that of H i gives T w 31000 K 
and fokin ~ 13kms~^, assuming the hydrogen and metals 
sample the same velocity field, which is not an obvious as- 
sumption to make. For a super Lyman Lim it System with 



logio A^Hi = 19.42, lO'Meara et 'all (|200l[ ) find from a si- 



multaneous fit to O I, N I and H I T = 1.15±0.02 x 10'' K 
and &kin = 2.56 ± 0.12 km s~^. In all cases, the tem- 
peratures obtained are consistent with photoionization 
heated gas. Although direct gas density measurements 
are not possible, models generally suggest that higher 
column density systems correspond to higher densities, 
especially at these high values. The temperatures sup- 
port the expected trend of decreasing temperature with 
increa sing density fo r sufficiently dense absorption sys- 
tems (|Meiksinl . fl993) . 

The small kinetic contributions to the line broaden- 
ing are in contrast to what is found for He ii absorp- 
tion systems , which appear to be dominated by kinetic 
broadening ( Zheng et all . l2004bf ). The He ii features, 
however, were chosen for their clear isolation in the 
spectrum, and so may probe physically underdense re- 
gions, in contrast to the higher column density systems 
above. An attempt to estimate temperatures from a 
larger sample of He ii and H i systems has not been able 
to produce secure values, but nonetheless shows that only 
~ 4 5% of the lines are dominate d by kinematic broaden- 
ing ( Fechner and Reimersll2007[) . 

Intergalactic metal absorption systems may be used 
to probe the time variability of fundamental physical 
quantities over cosmological timescales. The excitation 
of ionic fine-structure levels by the CMB allows metal 
lines in the Lya forest to be exploited for measuring the 
evolution of the CMB temperature ( Bahcall and Wolj . 
Il968l ). Because of the possible influence of other excita- 
tion mechanisms, the measurements are strictly speaking 
upper limits. Existing excitation tempe rature measure- 
ments are T = 7.4 ± 0.8 K at z = 1.776 (ISongaila et all . 
I1994D . T = 7.9 ± 1.0 K at 
T = 15.7 ± 3.5 K at 



1. 9731 (iGe et al. 



T 



2.34 ( Ge et al . 2001). and 



K (95% c.i.) at z = 3.025 ()Molaro et all . 



I2Q02D . These values are all consistent with a linearly 
evolving CMB temperature of T{z) — To(l -I- z), ex- 
pected for a homogeneou s and isotropic expa nding uni- 
verse, with To = 2.725 K (|Mather et a/.l . Il999f ). 

Detected resonance line multiplets have been used to 
place constrai nts on the time var iation of fundamen- 
tal constants. iBahcall et al. 1 (I1967D used the fine struc- 
ture wavelength splittings of the Si ii lines near A1260A 
and A1527A and the Si iv lines near A1394A to show 
that the fine structure constant (a = [l/47reo]e^/?ic) 
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varies by less than 5% between z = and z = 1.95. 
Tentative evidence for an evolving fine structure con- 
stant between 1 < z < 1.6 of Aa/a (— 1.9 ± 



0.5) X 10-5 was reported by IWebb et all (|l999|) us- 



ing measurements of Mg ii and Fe ii transitions. The 
value was later revised to Aa/a ~ (—5. 7 zt 1.1) x lO"* ' 
over the redshift range 0.2 < z < 4.2( Murphv et all . 

suggestive of a positive detection of a varia- 
tion. There is by no m eans concensus on t he de tec- 
tion. IChand et all (|2004[ ) and ISrianand et"aH (|2004l ) re- 
port Aa/a = (—0.6 ±0.6) x 10"^ over the redshift range 
0.4 < z < 2.3. Their analysis method has been con - 
tested bv lMurphv et all (l2007l) andlMurphv et all (|2008l ). 
and defended by Srianand et all ( 20071 ). In an inde- 
pendent effort using a mo dified technique designe d to 
eliminate some systematics, iLevshakov et al. 1 (|2007D ob- 
tain Aa/a = (5.4 ± 2.5) x 10"^ between z = and 
z = 1.7. Including the hyperfine splitting of the hy- 
drogen 21cm line permits constraints on the variation of 
X = a^gpfi, where fi is the electron to proton mass ra- 
tio and gp is the proton gyromagnetic ratio. The cur- 
rent best limits are, for a linear fit over the redshift 
range 0.2 < z < 2 . 4, i;/ a: = (-1.43 ± 1.27) x lO'^^yr 



( Tzanavaris et aL . 2005f). Combining H i 21cm data and 
OH 18cm data, iKanekar et al ] llool) obtain a limit on 
the effective combination F = gp{a'^ / fj.y-^'^ of AF/F = 
(0.44±0.36'''^*±1.0"y''') xlO-5 between z = and z « 0.7. 



C. Direct flux statistics 

While the absorption line parameters relate most di- 
rectly to the physical characteristics of the absorption 
systems, measuring similar parameters in numerical sim- 
ulations can be very computationally demanding, partic- 
ularly when building up sufficient statistics to make such 
a comparison meaningful. More straightforward compar- 
isons may be made using the flux statistics directly. The 
statistical measures in most use are the mean transmitted 
flux, the transmitted flux per pixel probability distribu- 
tion, the flux power spectrum and its Fourier transform, 
the flux autocorrelation function. An additional analysis 
based on wavelets permits a quick means of gauging the 
distribution of widths of the absorption features. 



1. Mean transmitted flux 

By far the simplest s tatistic is the mean t r ansmi tted 
flux through the IGM. lOke and Korvcanskvl ( 19821 ) in- 
troduced two measurements. Da, the mean flux of a 
background source removed by scattering between the 
Lya and Ly/? transitions in the restframe of the QSO (but 
avoiding emission line wings), and Db, the mean flux re- 
moved between Ly/? and the Lyman edge. These defini- 
tions are straightforwardly generalized to the fraction of 
flux removed between each successive pair of Lyman tran- 
sitions. Removing the effect on higher orders by lower 



orders, however, is not straightforw ard but may be ac- 
complished in a statistical sense (.Press et "^■ 119931) . The 
measurement Da may be evaluated as Da = (1 — fi/ci), 
where ft is the measured flux in pixel i, Ci is the cor- 
responding predicted continuum of the background ob- 
ject, normally a QSO, and the average is carried out 
over a given wavelength (or corresponding redshift) in- 
terval between Lya and Ly/3 in the restframe of the 
QSO. Considerable care must be taken in measuring 
the contribution at a given redshift from Lya alone. 
Systematics that would bias the result if not corrected 
for include: 1. a low H i fraction near the QSO due 
to its own ionizing intensity (the "proximity effect"), 
2. contamination by metal absorption systems, 3. the 
large stochastic influence of Lyman Limit Systems and 
Damped Lya A bsorbers, the latter of wh ich produce a 
divergent mean (jMeiksin and WhEi,[200l), and 4. evolu- 
tion of the absorption in the region measured, including 
over individual pixels of the spectrum. A detailed dis- 
cussion IS provided bv lTvtler et~all (|2004al) . 

For a spectrum that may be decomposed into indi- 
vidual absorption features with equivalent widths w(x) 
depending on some set of parameters x, such as column 
density or Doppler parameter, the mean of the transmit- 
ted flux %, = exp(— Tjy) is given by 



{%) = (exp(-T^)) = exp(-r/), 



(30) 



with the effective optical depth r; due to the lines given 

by 



n 
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dAA(x) 
dX 
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(31) 



where dN'{x)/dX is the number of systems per unit 
rest waveleii g th pe r unit vo l ume d x in parameter space 
( Press et al I 119931 : ISpitzed . Ulii), liere represented in 
terms of the number per (proper) equivalent width per 
unit redshift d'^N/dwdz. The variance (T^) - {Tf of 
the transmission depends on the details of the line pro- 
files. Express ions f or idealized examp les are provided by 
iMiinchl llMi) and lPress et~ai\ (Il993h . 

It is instructive to consider a comparison between 
Eq. (PT|) and the Gunn-Peterson optical depth (|15p . Us- 
ing Eq. (I27p for the number density evolution given by 
dN/dz = J dwd'^N /dwdz and Eq. (f24l) for systems on 
the linear part of the curve-of-growth, r/ may be recast 
in the form of the Gunn-Peterson optical depth 



n 



3 Af„r„,(n(z)) 



H{z) 



(32) 



for a mean distributed scattering gas density {n{z)) — 
'^T^ fiuna{z)(Ja{z)N [z] = ?>TT fluQ a{z)n{z) , whcrc N ^ nL 
is the mean column density of the absorption systems 
with mean atomic density n and line of sight thickness 
L, and Qa = na<JaL is the porosity of the absorbers. 
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equivalent to the spatial filling factor of the absorption 
systems for Qa < 1- For a column density distribution 
varying like dN/dNm cx A^h/ '', the effective absorber 
optical depth (excluding DLAs) is dominated by satu- 
rated systems, with F{to) « 1 — 2. In this case, Eq. ([2T|) 
gives w/Xo — 2{b/ c)F{tq) « 36/c. If the absorber veloc- 
ity widths scale like the Hubble expansion across them, 
b « fHL, Eq. dSIl) gives instead of Eq. 



Tl ~ 2Qaiz)f{z)F{To) ~ 3Qa{z)f{z). 



(33) 



In practice the effective optical depth will scale some- 
what between the linear and saturated line limits. It 
will not in general scale like the Gunn-Peterson optical 
depth Eq. as is often asserted in the literature. For 
dN/dz cx (1 -I- z)^, Eq. ((3T|l instead gives the dependence 
r; cx (1 + z)^^'' . This is corroborated by estimates of 
the mean H i transmission due to Lya scattering from 
all components, whether or not they may be modeled as 
individual Voigt profile absorption systems. More gener- 
ally, the mean transmission may be expressed as 



{%) = exp(-reff). 



(34) 



From an analysis of Lya scattering by the IGM in the 
redshift range 2.5 < z < 4 .3 using moderate resolution 
spectra, iPress et al\ (|l993h obtained Toff — A{1 + z)^^"' 
with yl = 0.0175 - 0.00567 ± 0.0002 and 7 = 2.46 ± 0.37. 
T he evolution expon ent 7 agrees closely with the finding 
of lKim et all ( 2002a) for dN/dz a bove. A similar result 
was obtained bv lKim et d\ (|200lh using high resolution 
spectra: A = 0.0144-0.00471 - y and 7 = 2.43±0.17(as in- 
ferred bv lMeiksin and White! (120041 ) ). iKim et al\ (l2002aD 
give the slight variation based on an expanded data set of 
A = 0.0032±0.0009 and 7 = 2.37±0.20. Because the im- 
plied probability distribution for (T) is skewed, the sta- 
tistical expectation value for the mean transmission must 
be computed from these expressions for Tp.ff using Monte 
Carlo integration (iMeiksin and Whitel . l2004l) . A good fit 
to the evol ution of TeS ba sed on the data compiled by 
iMeiksin and Whiti (12004), exte n ded t o lower redshifts 
using the data of lKirkman all ( 20071 ). is given by 



Teff = 0.0164(1 -I- z)i-i; 
= 0.00211(1 + z)^-^; 
= 0.00058(1 + z)'^-5: 



0<z<1.2 (35) 
1.2 < z < 4 
4 < z < 5.5 



(lMeiksinl . [2006al) . The evolution is found to vary from 
relatively slowly at low redshifts, as expected for discrete 
saturated features, to the more rapid evolution of the 
Gunn-Peterson effect at high redshifts as the porosity of 
the absorption systems approaches unity and increasingly 
diffuse gas dominates the effective optical depth. 

A dip i n the evolution of r^ g at z ~ 3.2 was re- 
ported by iBernardi et al. I ((2OOI) based on a sample of 
1061 moderate-resolution QSO spectra measured by the 
SDSS. The feature was interpreted as evidence for sud- 
den heating due to the onset of He 11 reionization. This 
interpretation is not supported by recent simulations 



( McQuinn et al. . 20081). T he feature was not confirmed 
bv [McDonald et all (|20(36l ) using an alternative analy- 
sis method applied to 3035 SDSS QSO spectra. Using 
a sample of 86 hig h-resolution high signa l -to-noi se ratio 
spectra, however. .Faucher-Giguere et al. (|2008b ) find a 
feature similar to that reported by Bernardi et al. ( 20031 ). 
The existence of the feature appears still not to be en- 
tirely secure, but tentatively suggests that a physical 
change in the state of the IGM occurred at z ~ 3.2, possi- 
bly reflecting a change in the H i ionization rate, perhaps 
owing to the growing influence of QSOs IIV.A.2"]) . 



2. Statistics based on pixel fluxes 

More detailed transmission information may be ex- 
tracted from the use of the full spectrum, including the 
probability distribution of the transmission per pixel, 
the flux power spectrum, the autocorrelation function of 
the transmission, and a wavelet decomposition. Because 
these statistical measures are best interpreted in conjunc- 
tion with numerical simulations, they are discussed in 
more detail in § IVIl Here they are only defined. 

The pixel transmission distribution (PTD) and its 
variants, like the trans mission fiux probability distribu- 
tion function (TFPDF) (jMcDonald et a/.l . l2000f) . simply 
quantify the frequency of occurrence of a given fiux value 
in each pixel of the spectrum, normalized by the contin- 
uum, oft en expressed in frequency b ins. It was intro- 
duced bv iJenkins and Ostrikerj ( 199l[ ) to search for un- 
derlying intergalactic absorption in excess of the amount 
accounted for by discrete absorption syste ms. Several 
tabulations are available in the literature: 



iRauch et all 

(Il997bl). iMcDonald et all (l2000l). iBecker et all {200% . 



iDesiacaues et all ~{ 2007 ) and Kirkman et al. (2007). 

Second order statistics include the flux autocorrelation 
function and the flux power spectrum, in analogy to sim- 
ilar statistical descriptions of the clustering of galaxies. 
The flux transmission au t ocorre lation function was intro- 
duced bv lZuo and Bondl (|l994f ) to investigate properties 
of the Lya forest. For a velocity separation v, it is defined 
by 



^^iv) = {f{v')fiv'-v))/{fr^i, 



(36) 



where f(v) denotes the pixel fiux at velocity v, 
and (/) is the measured mean flux at the rel- 
evant redshift. ^'^ Several estimates of are 



This differs somewhat from the definition of IZuo and Bondl 
(1994), who in place of (f)^ in the denominator, used the more 
general form {f{v')){f{v)). In most of the literature, it is as- 
sumed the mean flux changes negligibly over the velocity sepa- 
rations for which it is computed. This will break down for suf- 
ficiently large velocity separations that cosmological expansion 
becomes a factor, but in this case the physical interpretation of 
the autocorrelation function becomes less transparent as it would 
correlate transmissions at different cosmological epochs. 



20 



provided in the literature: Zuo and Bondl (Il994) 
iLiske e t ali (l 2000D. iMcDonald et al\ (|2000D. ICroft et al\ 



(12002b '). Bec ker et al 



(I2004D. iD'Odorico et al\ tOO(t) . 

iMcDonald et Id] (|2006l ) and lKirkman et all (|2007l) . The 
results generally show significant correlations over scales 
of a few Megaparsecs (comoving) . Similar to the flux au- 
tocorrelation function is the flux cross-correlation func- 
tion along neighbo ring lines of sight t o QSOs near each 
other on the sky 
l2006t iLiske et all 



Becker et all 120041 : ID'Odorico et all 



2000) ■ Cross-correlations of similar 



strength to the autocorrelations are found, extending 
over similar spatial scales. 

The flux power spectrum has received considerable at- 
tenti on as a mea i is of constraining cosm ological parame- 
ters (jCroft e^"a/ ], |1998'; 'Hu i et Qi.U 200l'^ . It is a measure 
of the variance in the amplitude of the Fourier transform 
coefficients of the transmitted flux. Defining 6f = f — {f) 
and its Fourier transform 



1 



Av/2 



Av/2 



dvSf{v)e'^\ 



(37) 



over the velocity interval Av, the flux power spectrum is 



PAk)^iAv){\sy\')/{fy^ 



(38) 



The dimensionless form kPp{k)/7r is often more conve- 
nient. It forms a Fourier transform pair with the flux 
autocorrelation function through 



kPpik) k 



Av/2 



J-Av/2 



ikv 



(39) 



The reader should be aware that several different con- 
ventions for PF{k) exist in the literature, some without a 
mean flux normalization, some differing by a factor of 2tt, 
and some referring to an effective 3D flux power spectrum 
extracted from the ID flux power spectrum in a manner 
analogous to the relation between the ID and 3D mass 
power spectra. Many tabulations of Ppik) arc a vailable: 
ICroft al\ (Il998i) ICroft et al\ (119991). IMcDonald et al\ 
(|2000l ). ICroft et al\ (|2002bD and lKim et al\ (|2004allbl ). A 
substantial effort was made to measure Pf(^) from the 
spectra of high reds hift quasars discovered by the Sloan 
Digital Sky Survey (|McDonald et all\200^ . 

A shortcoming of the Fourier transform is that it loses 
spatial information. In the case of a spectrum show- 
ing the Lya forest, the positions of individual absorp- 
tion features are lost. A wavelet transform is analogous 
to a Fourier transform, but retains local information, 
giving in effect a description of the power on different 
scales as a function of position. A wavelet is a square- 
integrable function, 'ip{x), defined in real space such that 
-ipjk = 2^/'^'ip{2^x — k) (where j and k are integers) forms 
an orthonormal basis for the set of square-integrable 
functions. It satisfies J^^dx'ip{x) = 0, and is normally 
concentrated near x — k2~^ . The wavelet coefficients are 
defined by wjk — J dxf{x)ipjk{x); the original function 
may be recovered through f{x) — J2j k'^3k''Pjkix)- The 



wavelet transform effects a localized multiscale rendition 
of the data. Although not commonly used, wavelets are 
in principle a powerful tool for analysing the Lya forest, 
producing a statistical description similar to the distri- 
bution of 6-values. The discrete wavelet transform, like 
the fast Fourier transform, may be computed extremely 
rapidly, permitting a large number of spectra to be an- 
alyzed much more quickly than by performing a full ab- 
sorption line analysis. This is an advantage over full 
Voigt profile line fitting for assessing the mat ch of sim- 
ulations to observations ([M eiksin, 2000; Mc iksin et al. 



20011: iMeiksin and White. 2001: ITheuns et all l2002a : 



Theuns and Zaroubii . .2000: ■Zaldarriagai |2002| ). 



III. IONIZATION EQUILIBRIUM 
A. Ionization 

1. Hydrogen and helium 

The absorption signatures of the IGM are produced 
by atoms and ions with bound electrons, either through 
resonance scattering or photoelectric absorption. The 
amount of scattering or absorption depends on the ioniza- 
tion structure of the IGM. In this section, the photoion- 
ization of hydrogen and helium is discussed. In hot en- 
vironments like galaxy clusters, coUisional ionization can 
dominate. CoUisional ionization equilibrium may play a 
particularly impor t ant role in t he physics of the WHIM 
( Cen and Ostrikeil I1999L 120061 : iDave et all l200ll . I1999D 
and is treated more generally in the section below on 
metal ionization. 

Photons with energies exceeding the ionization poten- 
tial of a bound electron in a hydrogen atom will ionize 
the neutral hydrogen atoms in the IGM at the rate per 
neutral atom 



HI 



J "'f HI 

diy- — , 
n.pi' 



(40) 



where a™ is the hydrogen photoelectric cross section, i/t 
is the threshold frequency required to ionize hydrogen 
(the Lyman limit), and u^, is the specific energy den- 
sity of the ambient radiation field. The specific energy 
density is related to the specific intensity of the radia- 
tion field Ii,{r,t,h) by inJ^/c, where J^(r,t) — 
(l/47r)^ dni^{r,t,ri) is the angle-averaged specific in- 
tensity. 

Free electrons will be radiatively captured by the pro- 
tons at the rate per proton rieaHii(T), where anii, often 
referred to as the Case A radiative recombination coef- 
ficent a^(r), is the total rate coefficient for radiative 
capture summed over recombinations to all energy lev- 
els. Expressed in terms of the ionization fractions xhi = 
n-Hi/"^H and a;Hii = '^hii/"h, where nn = nm + nmi is 
the total hydrogen number density, the ionization rate 
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equations become 

dxm 
dt 

dxmi 

dt ~ dt ' 

The time-derivatives are lagrangian, so that Eqs. 
vahd in the presence of velocity flows. 

In equilibrium, the neutral fraction will be 

neaA{T) 



cq 



HI 



(41) 
T|) are 

(42) 



For a pure hydrogen gas, — nmi, and the rate equation 
for xm becomes 



dt 



'"rec-'-HI 



(r 



HI 



2t^^l)xm 



C = 0, (43) 



where t^cc = I/I^-hc^a) is the total radiative recombina- 
tion time for the gas. The equilibrium neutral fraction 
may be expressed as 

1/2 



"HI 



1 



-0+1 - 



1 



- 1 



(44) 



where (jj — FHi^rcc denotes the number of photoioniza- 
tions per neutral atom over a recombination time. It is 
instructive to consider the case of constant cf). Eq. (HH)) 
then has the general solution 



cq 

a;Hi = -Thi 



V 



(1 - ^0 ^r]) exp(77r) - 1 



(45) 



where 77 



-2-x' 



HI' 



t/trcc, and So = a;Hi(0)-XH'*j(0) 



is the initial deviation of the neutral hydrogen frac- 
tion from its equilibrium value. For a fast photoion- 
ization rate {(p ^ 1), x'^^ ~ (/)~^ ^ 1 and xm — 
+ <5o exp(— Fnit). An initial deviation Sq from equi- 
librium decays exponentially. The neutral fraction ap- 
proaches its equilibrium value on the photoionization rate 
timescale F^j, regardless of whether equilibrium is es- 
tablished through net photoionizations or net recombina- 
tions. Ionization equilibrium will be maintained provided 
any change to the gas (such as source variability or the 
expansion of the gas) occurs on a time scale long com- 
pared with the photoionization time scale. By contrast, 
for a slow photoionization rate 0^1, ~ 1 and the 
neutral fraction approaches equilibrium on the recombi- 
nation timescale, xm — 1 + (5o exp(— t/trcc) for \So\ ^ 1 
(<5o < 0). 

The photoionization equations for helium are analo- 
gous to Eqs. (|1T|) . Denoting the respective He i. He 11, 
and He 11 fractions by ^Hei) ^^Heii and ccHeiii; they are 
given by 



dx 



He I 



dt 

dXHc II 

dt 
(^aJHoiii 



dt 



— XHoiEhcI + a:HoII'^eQ!HoII, 
rfl'Hc I dX}ic III 

Jt Jt ' 

XHe IiEhs II — a;He III^^ettHe III , (46) 



where auc ii and auc iii are the total recombination rates 
to all levels of He i and He ii, respectively, and FhcI and 
FhoII are the respective photoionization rates. The total 
electron density is Ue — nmi + ?^HcII + 2nHciii- 

The Case A radiative reco mbination coeff icient for hy- 
drogenic ions is provided bv ISeatonI (jl959l ). from which 
the coefficients for q;hii(T) and aHciii(^r) may be de- 
rived. (In general, the radiative recombination coefficient 
az{T) for a hydrogenic ion of atomic number Z is related 
to that of hydrogen by az{T) = ZaHn{T/Z^).) For el- 
ements He and beyond, the recombination coefficient is 
the sum of two rates. One is the usual radiative capture 
term, in which e~ + X""*"^ X" -|- 7 for a species X 
with n electrons. The other arises from dielectronic re- 
combination, in which a fast incoming electron excites an 
ion, which becomes trapped in a weakly bound autoion- 
izing state upon capturing the electron, then decays by 
the emission of a resonance line photon: e~ + Ar"+^ 
AT"* A" -I- hpiy^i. The probability for emission of the 
photon is suppressed by the factor exp(— /ipj/ui/fceT); 
usually excitation is followed by the reverse reaction. 
A form for the radiative recombination rate coefficient 
coi itribution to aaeuiT) for He 11 ^ He i is provided 
by IVerner and Ferlandl ( 19961 ) , while a form for the di- 



electronic term is provided bv lAldrovandi and Pequigno^ 
(1973'). The recombination coefficients are listed in Ta- 
ble mil Recombination to He i may also occur through 
charge transfer, He+ + H° ^ He" + H+ + AE, where 
AE is the energy defect. Where hydrogen is highly ion- 
ized, this is a negligible effect compared with radiative 
recombinations. For partially ionized hydrogen, how- 
ever, as for example within a Damped Lya Absorber, 
or the IGM prior to its complete reionization, the ef- 
fect may pla y an important role. The ra te coefficient is 
provided bv iKingdon and Ferlandl ( 19961 ). The reverse 
process, charge transfer ionization, is suppressed by a 
Boltzmann factor and so is normally negligible except at 
high temperatures. 

The exact photoionization cross section for hydrogenic 
atoms with atomic number Z in the ground state is 
(,Heitler. . , 1954.1 



, , Ao fVT\^ g4-(4tan 'e)/e 



(47) 



for hpi/ > hpi^T — 13.60Z^eV, the threshold energy for 
ionization, and where Aq = (2^7r/3e^)a7rao = 6.30 x 
10^^^ m^, where oo = [47reo]?i^/TOee^ is the Bohr radius, 
a — [l/47reo](e^/?ic) is the fine structure constant, and 
e = {v/vt — ^Y^"^- For most practical applications, how- 
ever, an approximate form is preferred. The ionization 
cr oss sections are we ll-approximated by the form given 
by[Osterbroc3([l98^, 



a{v) = ot 



(1-/3) - 



-s-1 



(48) 



The parameters for hydrogen and helium are listed in 
Table IIVI For results requiring very high accuracy, more 



22 



TABLE III Case A recombination (a^ [m'^s ^]) and recombination cooling {(3i [J m^s ^]) coefficients. 



anil 2.065 x lO^^^T^s (6.414 - i In T + 8.68 x lO^^T^ 



aHcii 3.294 X 10" 



• T 'I 
, 15.54; 



1 + (t5^) 



1 + 



„ / -9.4x10" 

1.9 X 10"^ 1 + 0.3e T — 



-4.7x10-" 3 
; T T"'? 



QHelll 8.260 X 10-"r-3 


V.107-ilnr + 5.47x 10" 






(3uu 2.851 X 10~'"'ri 


5.914 - i In T + 0.01184r^ ' 





/3hcii 1.55 X io-39r"-^«''^ 



1.24 X 10" 



1 + 0.3e- 



-4.7X 10" 3 

e T r"2 



Puoui 1.140 X 10-^^r2 ( 6.607 - i In T 7.459 x lO^^Ts 



TABLE IV Ionization cross-section parameters used in Eq. 


(|48)). From'Osterbrockl (|l989ri. 






aT[m ^] 


[Hz] 


/3 


s 


ani 6.30 x IQ-^^ 
OHei 7.83 X 10-22 
OHell 1.58 X 10'22 


3.282 X lO^'^ 
5.933 X lO^'^ 
1.313 X 10^** 


1.34 
1.66 
1.34 


2.99 
2.05 
2.99 





precise fi ts to the phq t oioniz ation cross sections are pro- 
vided bv lVerner 6^0!! (|l996l ). including results for an ex- 
haustive list of metals. 



2. Metals 

The large number of ionization states of metals results 
in more complex processes involved in their photoioniza- 
tion. These include autoionization and Auger electron 
production in addition to dielectronic recombination. In 
an autoionization, an excited ion with one or two elec- 
trons lying in the photoionization continuum is produced 
following a photoionization event in which a deep shell 
electron is ejected. The excited ion spontaneously ion- 
izes through a radiationless transition and the ejection of 
multiple Auger electrons. Although there has been con- 
siderable progress in accurate estimates of all the relevant 
rate coefficients, many remain poorly determined either 
experimentally or theoretically, requiring interpolation 
(and sometimes ex trapolation) from c oeffic ients known 
for related species ()Ferland and Savinl . l200lf ) . These un- 
certainties can affect the predicted metal abundances in 
the absorption systems. The effect of uncertainty in the 
dielectronic recombination rat es on c ommo nly measured 
metal features is discussed bv lSavinI (|2000D . Some of the 
methods of approximation are described in the literature 
for the rate coefficients (see below). 

Sophisticated computer codes have been developed to 
solve for the equilibrium ionization states. A widely used 
publi cly-distributed code is CLOUDY ( Ferland et all 
[1^998) , which solves the radiative transfer problem for 
plane-parallel and spherical geometries. The photoion- 



ization of i ntervening absorption syst ems has been de- 
scribed by IVerner and lakovlevl (1993), whose formula- 
tion is followed here. 

Denote by Xij the fraction of a particular element i in 
the j"^ ionization state, so that j runs from 1 (neutral) 
to Zi + 1 (fully ionized) for an element of atomic number 
Zi. Then X]7"=T^^y = 1- The general rate equation for 



HeiCij + Oy + Ty-i + dy_l) + Ptj 



Zi + 1 
k=]+2 

+ nHI?'5-i + "HclT-fj-i 



ne(cy_i -t-ay_i) +Pij-i 



'^HIIC^-I + '^HoIlCy'li 



i-2 

1^1 



(49) 



where 



j is the collisional ionization rate 
(IVoronovl Il997l). Oil the autoionization rate 
( Arnaud and Rothenfluel . Il985[ ). Pij the pho- 
toionization rate without post -ioni zation 

Auger-electron ejection ( Verner et al\ . Il996t 
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Verner and Yakovlevl. 1995). 


r,-, the radiative re- 


combination rate (lAldrovandi and Peauignot. 1973 


Arnaud and Rothenflue. 19851: 


ShuU and van Steenbere 



the dielec- 



troni c reco i nbinat ion rate ( Nussbaumer and Storevl . 
1198.1 Il986l hmih . c] 
ionization with H ii 



the rate of charge t ransfe r 
(Kingdon and Ferland], Il996f) . 
the rate of charge transfer ionization with He ii 



( Arnaud and Rothenfiuel . Il985[ ). r.^ the rate of charge 



transfer recombination with H i (Kingdon and Ferlandl, 



1996i), and r^" the rate of charge transfer re combination 



with He I ( Arnaud and Rothenfiuel . Il985t ). Photoion- 
ization from deep shells followed by Auger-electron 
ejection, resulting in a change from an initial ionization 
state 7 to a final state k > j + 2, occurs at the rate 
Qijk I Kaastra and Me"w3 . Il993l ). The convention that 



corresponds to recombination to ionization state j 
and Cij corresponds to coUisional ionization of state j 
is used. Because of the low metallicity of the IGM, the 
electron density may normally be well-approximated by 
the electrons lost only by hydrogen and helium. 

Using the above notation for the various species and 
taking i — 1 ior hydrogen and z = 2 for helium, Eqs. (HT1) 
and (|46p are modified, allowing for coUisional ionization 
and the charge transfer reactions, to 
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dxi2 
dt 
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dxii 
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P21 + ne(c21 + 021) + 9213 



EE 

i#2 3=2 



„Hc 



dx 



22 



dt 

dX23 

dt 



-X22 ne(r2i -1-^21) + Ey^ 

1^2 j=l 
dX21 dx23 

~dt df' 



"'^3^1] 



-X23 



+X22 



«e(?'22 + ^22) + niir^2 +a;2l9213 
"e(c22 + 0.22) +P22 + ni2C^2 



(50) 



where the first sum of the double sums is over all elements 
except the excluded one indicated. 

Noting that the photoionization rates are proportional 
to the energy density of the radiation field, it follows 
from the structure of Eqs. (|49)) and ((50|) that the equilib- 
rium fractions are determined by the ratio U = nj/nn 
for a given spectral shape of the radiation field with pho- 
ton number density n^, and for a given set of elemental 



abundances. The ratio U is know as the ionization pa- 
rameter. Larger values for U generally produce higher 
ionization stages in the metals, although these also de- 
pend on the spectral shape of the radiation field. Com- 
parisons with measured ion column densities may then 
be used to constrain the natures of the radiation field 
and the absorption systems. 

Fitting formulas to coUisional ionization, radiative re- 
combination and diele ctronic recombination ra t e coeffi - 
cients are provided by IShull and van Steenberel (|l982at ). 
Improved values for many of these rates have since 
then become available, and have been referenced above. 
Power law fits to coUisional and radiative recombina- 
tion rates involvin g hydrogen and helium are provided by 
lAbel et al\ (|199<D . Fitting formulas to more recent rates 
for the elements hydrog en through nickel are provided 
bv lMazzotta et oZI (1998), along with equilibrium ioniza- 
tion fractions for metals in a hot environment, for which 
coUisional ionization will normally dominate photoion- 
ization. A comprehensive compilation of the literature 
is available at D. Verner's Atomic Data for Astr ophysics 
website, as described in iFerland et al\ ( 1998t ). Since 
then, improved estimates have been obtained for some 
rates, including laboratory measured values of dielec- 
tronic recombination rates ( Savin et al. 1 (|2006D and ref- 
erences therein) . A comprehensive review of the atomic 
physics relevant to x-ray lines, as will arise from a war m- 
hot IGM, is provided bv lKallman and Palmeril ( 20071 ). 



B. Thermal equilibrium 

In the presence of heat transfer at the rate (G — L)/n 
per particle, where G and L are the thermal gain and 
loss functions per volume and n is the particle density, 
the second law of thermodynamics requires the entropy 
per particle s to change at the rate 



(51) 



where T is the temperature of the system. For an ideal 
gas, the entropy per particle is s = (7— l)~^fcB ^i^{p/p'^) + 
So, where p is the gas pressure, p is the mass density, 7 is 
the ratio of specific heats at constant pressure to constant 
volume (7 = 5/3 for a monatomic gas), fcs is Boltzmann's 
constant, and so is an arbitrary additive constant. For 
numerical computations, it is often useful to express the 
energy equation in its entropic form Eq. (|5ip . but through 
the entropy parameter 

Se^^. (52) 
It follows from Eq. ([?T|l that 



dSE 
"dt 



= (7-l)p-^(G-L). 



(53) 



http: / /www. pa. uky.edu/~vcrncr/atom. html 
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The gas temperature is related to Se by 



T 



-Sep 



,7-1 



(54) 



where rh is the mean mass per particle. 

Photoionization will provide heat through the excess 
energy absorbed by an electron above the threshold en- 
ergy hpVT required to ionize a given atom. For a single 
species i of density n^, the photoionization heating rate 



(55) 



where hpv^ is the ionization potential of species i. The 
total heating rate from all species is 



G — Ghi + GhcI + GhcII- 



(56) 



The contributions due to the photoionization of the met- 
als is normally negligible because of their low abun- 
dances. 

Cooling is provided by recombinations, coUisional ex- 
citation of the excited levels in neutral hydrogen, and 
inverse Compton scattering off CMB photons. At high 
temperatures, free-free losses and collisionally excited 
line radiation from atoms and ions other than neutral 
hydrogen may become important contributions as well. 
As for heating, all species contribute to cooling, giving 
the total cooling rate 

L ~ Lhu + Liicn + Lhciu + LeH + Lc + Liinc + Lf f . (57) 

For a single species, recombinations radiate the elec- 
tron energy ~ k^T as photons at the rate Li = 
nenij3i{T). Recomb i natio n cooling coefficients j3i{T) are 
provided bv lSeatonI ( 1959( ) for recombinations to H i and 
He II. An expression for the radiative contribut ion to 
/ ^Heii for recombinations to He i is provided by iBlackl 



(Il98lh . The dielectronic con tribution may be approxi - 
mated as = 3Ry a'^l\i (iGould and Thakuii (l970l ) 



using the second part of the expression for ancii in Ta- 
ble [HTI which expresses the dielectronic component. The 
total recombination cooling coefficients are provided in 
Table Uni 

The cooling rate due to the collisional excitation of H i 
by electrons is given to an accuracy of 3% over the tem- 
perature range 4000 < T < 12000 K by the expression 
(ISpitzeil .lT978): 



L^H = 7.3 X 10"^2 J 



ngJiHie 



-118400/T 



(58) 



For temperatures relevant to the IGM, cooling losses due 
to collisional ionization of hydrogen and collisional ex- 
citation and ionization losses from helium are negligi- 
ble. Similarly, for the low metallicities encountered in 
the IGM, radiative losses from metals are normally neg- 
ligible, although they may possibly be important in the 
WHIM. 



Cooling is also produced by the Compton scattering 
of CMB photons off the electrons, which is es pecially im- 
porta nt at high redshifts. The cooHng rate is ( Wevmainil . 
I1965D 



^a^CMB(^)^B [T - Tcmb(^;)] , 



(59) 



where ctt = {&TT/3){e^ /[ATTeo]meC^y = 6.65 x 10' 
is the Thomson cross section, a is the radiation density 
constant, nie is the electron mass, c is the speed of light, 
and TcMB is the temperature of the cosmic microwave 
background. 

In the WHIM, cooling by thermal bremsstrahlung ra- 
diation, or free-free emission, due to the acceleration of 
a charge in the Coulomb field of another charge may be 
significant. For an ion of number density rij and ch arge 
Zi, the cooling rate is ( Rvbicki and Lightmanl . [i979^ 



Shapiro and Kani 



3 (| 19871) 



gs = 0.79464 

10^ K and gs = 



provide the fitting for- 
0.1243 logio(T/Z2) for T/Zf < 
2.13164- 0.1240 logio(r/Z,f) for 



Lff = 1.43 X 10""*° J m^s-^T^/^nen,ZfgB, (60) 

where gs is a frequency average of the velocity- 
aver aged_Gaunt_^£toi\ and ranges between 1.1 and 
1.5. 
mula 
3.2 X 

T/Zf > 3.2 X 10^ K. Cooling rates including collisionally 
excited metallic line losses and their ind i vidual emissivi- 
ties are provided bv lCaetz and Salpeteij ( 19831 ) for solar 
abundances. Non-solar values may be interpolated with 
the zero-metallicity limit for which atomic cooling is de- 
termined only by the hydrogen and helium losses. 

Photoelectric heating by x-rays is automatically in- 
cluded in Eq. (|55p . X-rays, however, may introduce two 
other sources of energy exchange: gains through Comp- 
ton heating and losses through the effect of secondary 
electrons. 

Sources of high energy radiation, such as a hard x-ray 
background, may inject addi tional energy into th e IGM 
through Compton h eating fArons and McCrav, 196g 
Madau and Efstathiou . 1999: Rces, 1969; Recs and Scttj, 
1970l ). For an IGM of electron density and tempera- 
ture T, an x-ray background with energy density will 
heat the gas through Compton scatterings at the rate 



Gr7 = 



rUeC' 



/•OO 

■ / dvu^{hpv - AkBT). 
Jo 



(61) 



For photons with energies exceeding 100 keV, Compton 
scattering takes place in the relativistic regime, requir- 
ing the Klein-Nishina cross section in stead of the Thom- 
son. The heating rate t hen b ecomes ( Blumenthall . 1197^ 
iMadau and Efstathioul . ll999l ) 



Gc = -n-eCCTT 



1 





-lOx" 



dxu„ 



2x- 3 



2x 

blx^ + 93x2 + 512; + 9 



ln(l + 2x) 
(62) 



3(2-|-3a;)3 

where x = hpv/meC^ . The Compton heating rate Gc 
should then be added to Eq. 
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The secondary electrons ejected by x-rays have only a 
small effect on the temperature of highly ionized gas. In 
partially ionized gas, however, their contribution to the 
energy losses may become appreciable. They play a par- 
ticularly important role in the cooling of gas ahead of ion- 
ization fronts during the Epoch of Reionization. Expres- 
sions for the energy of a primary electron deposited as 
heat, further ionizations, and coUisional excitations, for 
both hydrogen and helium, in a partially ionized medium 
are provided bv lShuU and van" Steenberd (I1985D . 



IV. THE METAGALACTIC UV BACKGROUND 
A. Mean energy density of the UV background 

1. Origin of the UV background 

Since the ionization structure of the IGM is determined 
by the specific radiation energy density Ui, of the meta- 
galactic ultraviolet photoionizing background, a complete 
understanding of the structure of the IGM requires an es- 
timate for Wi,. In turn, measured properties of the IGM, 
like its optical depth to hydrogen and helium ionizing 
photons or the ratios of ionization fractions of metals, 
place constraints on the shape of the metagalactic UV 
background. 

In a homogeneous and isotropic expanding universe, 
Eq. (O gives for the integrated specific energy den- 
sity produced by sources with a total proper emissivity 



1 



dl 

(63) 

where Tcs^i', z, z') is an effective optical depth due to 
absorption by the IGM and i^' = + z')/{l + z). 
For discrete absorption systems, Eq. pTjl may be used, 
adapted to photoelectric absorption by choosing x = Nm 
and taking w — J dX"[l — exp(— r^")], where — 
i^{l+z")/{l+z). Then, on replacing dA"aW/(aiVHi9A") 
by dz"d^JV/{dNmdz"), 

Teff(z/, 2, z') = / dz" I dNiii ^^^ ^ ^J l-exp(~T^//)], 







' dNmdz 



where 



(64) 



T^n = iVHiaHi(i^")+^HciaHci(i'")+^HciiaHcii(!^")- (65) 

The emissivity consists of two contributions, iso- 
lated radiation sources, such as galaxies and QSOs, 
and diffuse emission from the IGM itself. Survey es- 
timates for the contributions from isolated sources are 
discussed below. There are four principal contribu- 
tions to the diffuse emission capable of photoioniz- 
ing the IGM ( Haardt and Madaul . Il996l ): 1. recombina- 



Lya recombination radiation, 3. He ii two-photon con- 
tinuum emission, and 4. He ii Balmer continuum emis- 
sion. The specific recombination emissivity associated 
with each of t hese processes may be ex pressed generally 
in the form (|Haardt and Madaul . Il996l ) 



er{i', z) — hpi>f{i>)B{i',z) 



dz 
din 



(66) 



where 



B(v,z) = C rfjVm PrccPcsc(t^) 

Jo dNmoz 

X /°°rf^'^u;ab.(^'), (67) 
Jut hpiy' 

where lyr refers to the photoionization threshold for ion- 
izing the ground state of an ion of atomic number Zi. The 
function B{v,z) expresses the escape rate of recombina- 
tion photons produced by discrete absorption systems in 
terms of the product of the fraction Wahs(v) of photons of 
energy hpv absorbed by a given ion, the probability Proc 
that a recombination will result in the particular radia- 
tive transition considered, and the probability pcac that 
the emitted photon will escape the absorption system 
that produced it. These factors may be computed from 
detailed ra diative transfer computatio ns using a code like 
CLOUDY (jHaardt and Madaul . fl99l . 



The capture of electrons into the n L level of a hy- 
drogenic ion results in the emission of a continuum pho- 
ton distribution sharply peaked at the photoelectric edge. 
The rate of direct recombinations may be computed from 
the Milne relation between the free-bound recombination 
cross section afb{v) for an electron moving at velocity v 
and the bound-free photoion ization cross-section a„(z^) 
(jRvbicki and Lightmanl . Il979l) 



2gn 



9i 



(68) 



tions to the ground state of H i and He ii, 2. He ii 



where g„ = 2n^ is the statistical weight for the atom with 
electron in the bound state n and gi = 2 is the statistical 
weight for the ion. The bound-free cross-section from 
state n in a hydrogenic atom of atomic number Zi may 
be expressed as 

647r n 2 f^T\^ ^ ,„r^\ 
"^"(^^ = ^^""0 [—) Gn, (69) 

where G„ is a Gaunt factor and a = [l/47reo]e^/?ic is the 
fine-structure constant. The Gaunt factor for the ground 
state is 

--(4tan-ie)/e 

G,=8.3'/\l + er' \_^_,^,, , (70) 

where e = (j//j/t - 1)^/^ (cf Eq. gT]). The pho- 
ton released following recombination will have en- 
ergy hpiy — hpVT + (l/2)meW^, where nie is the 
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mass of an electron and hp is Planck's constant. 
For a Maxwellian velocity distribution — 
exp(— meW^/2fcBT) for the elec- 
trons, the free-bound emissivity due to recombinations 
to state n in a hydrogenic ion of atomic number Zi and 
number density n; is then 



nen,,hpV(jfb{v)vfM{v)-^ 



An 

— neUi 



X exp 



hi 



dv 
dv 

3/2 



2'KmekBT 
hpv — hpVT 



hpv an{v)2n 



(71) 



Recombinations to He ii produce hydrogen-ionizing 
photons [hpv > 13.6 eV). The recombinations gener- 
ally occur to high energy levels, followed by a cascade of 
bound-bound transitions. Recombinations that populate 
the 2^P level of He ii result in a He ii Lya photon of en- 
ergy 40.8 eV, which photoionizes neutral hydrogen in the 
Lya forest systems. Recombinations that populate the 
2^5* level produce He ii 2^5* l^S two-photon emission 
with a total energy of 40.8 eV and a spectrum peaking at 
20.4eV. Direct recombinations to 2^S and 2^P produce 
He II Balmer continuum emission, for which the thresh- 
old photoionization energy is the same as the hydrogen 
Lyman limit. 



2. Sources 

QSOs were early recognized as a promi- 
nent source of metaga lac tic ionizing p hoton s 
(|Arons and McCravl Il970al lbl: iRees and Settil . Il970l ). 
Since then, there have been numerous estimates of 
their contribution to the UV metagalactic back 
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(|Bovle et all [1988, ') 



lQ0.4(l-/3i)(M-Af*) IQ0A{1-I32)(M-M') ' 

(72) 

where 0(M, z)dM describes the spatial density of QSOs 
at redshift z in the absolute magnitude range M to 
M -\- dM, (p* is a spatial normalization factor (normally 
expressed in comoving units), and M* is a characteris- 
tic break magnitude. In a Pure Luminosity Evolution 
(PLE) model, (/)* is held fixed while M* is allowed to 
evolve. In a Pure Density Evolution (PDE) model, M* 
is held fixed while (p* is allowed to evolve. Although 
PLE and PDE models well fit t he data over limited 
redshift ranges. iFan et all ( 2001b[ ) show that the steep- 
ness of the bright end of the luminosity function has 



evolved between z < 2 and z > 3.6. Estimates of the 
QSO lumino s ity fu nction at z < 2.1 are provided by 
ICroom et d\ (I2004D ba sed o n the 2dF redshift survey, 
and by Richards et al. (200 ^, com bining data sets from 
the 2dF and SDSS. iRichards et al\ ^06) use the SDSS 
to measure the bright end of the QSO luminosity function 
and its evolution over the redshift range < z < 5. 

A long-standing uncertainty affecting estimates of the 
QSO contribution to the UV background at high redshifts 
(z > 2) is the poorly determined number of low lumi- 
nosity QSOs. iHunt et al . (2004) have measured the low 
luminosity end of the luminosity function at z « 3. An 
estimate of the QSO luminosity function over the redshift 
range 3 < z < 6 is provided by Meiksin (2005i" ) based on 
fitting PLE and P DE mo dels to the resul t s of Fan et aJ\ 
pOOlbL.Fan et~d\ (|2004l ) and lHunt et~al\ (l200j) An es- 
timate of the bolometric QSO luminosit y function over 
the re dshift range < z < 6 is provided bv lHopkins et al\ 
(|2007D . 

The QSO emissivity is given by = 
J dMip{AI)L,^{AI), where L^{M) is the mean spe- 
cific luminosity corresponding to a QSO of magnitude 
A'/. The spectral shape of QSOs, however, is uncertain 
at frequencies above the H i Lyman edge. Results based 
on HST data using the Faint Object Spectrograph 
(FOS) give an intrinsic spectral distributio n fi, ex i/~°Q 
with QSO spectral index aq = 1.76 ± 0.12 (jTelfer et al\ . 
120021 ) over the restframe energies 0.8 — 1.8 Ry. The 
sample contains QSOs predominantly in the redshift 
range 0.3 < z < 2.3. Distinguishing between radio-quiet 
and radio-loud QSOs gives aq = 1.57 ± 0.17 and 
aq = 1.96 ± 0.12, respectively. Results based on FUSE 
data give instead aq — 0.56lg'3g (without di stinguishing 
betw een radio-quiet and radio-loud QSOs) (jScott et al\ . 
I2004D . The origin of the discrepancy is unclear. The 
FUSE sample is at lower redshift (z < 0.67) and the 
sources have lower luminosities, so the difference may 
reflect a hardening with decreasing reds hift and /or 
lumin osity, for which there is evidence ( Scott et al\ . 
l2004[ ). The spectral range used for the power-law fits 
to the FUSE spectra is restricted to restframe energies 
0.8 — 1.4 Ry, similar to the range for the HST sample. 
The spectral corrections due to absorption by the 
intervening IGM, however, are smaller for the FUSE 
sample QSOs than the higher redshift HST QSOs, so 
that the FUSE results may be more reliable. 



The angle-averaged intens ity Ju 



,/47r of t he UV 



background determined by iHaardt and Madad (|l996f ) 
based on an earlier assessment of the QSO contribution is 
shown in Fig. [5] at redshifts z = 2 and 3. The steps in the 
background are produced by photoelectric absorption of 
the source radiation by intervening absorbers at the H i 
and He ii Lyman edges. The comparison with the purely 
absorbing IGM estimate shows the formation of emission 
"horns" due to the reprocessing of ionizing radiation by 
the IGM into redshifted H i and He ii Lya photons. The 
redshifted He ii Lya photons enhance the flux of H i- 
ionizing photons as does the continuum emission result- 
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TABLE V Estimates of the H i mean Lya transmitted flux (co l 2) and H i-io nization rate Fjjf^ (in 1 '^'^ s ^) required to 
match the mean Lya flux according to lMeiksin and White! (|2004l ) (col 3) and to lBolton et al\ (|2005l ) and lBolton and Haehneltl 
(|2007l ) (col 4). Columns 5 and 6 are the H i-ionization rates rgf°'^^^ and rgf°'™| ^ (in 1 0"^^ s'^) from QSOs using the Pure 
Luminosity Evolution (PLE) and Pure Density Evolution (PDE) models of iMeiksinI (j2005l ). Column 7 gives the H i-ionization 
rate F^^^' using the emissivity from iHopkins et cU. ] poo?) and the same attenuation lengths as for the PLE and PDE 
models. For the PLE and PDE models, the indicated redshift entries z = 3.89, 4.00, 5.00 and 6.00 correspond to the ranges 
3.6 < z < 3.9, 3.9 < z < 4.4, iA < z < 5.0 and z > 5.7, respectively. 
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FIG. 9 The UV metagalactic ionizing background inten- 
sity Ju , normalize d by the intensity Ji at 1 Ry, from 
iHaardt and Madaul ^9M) at z = 3 (thick line) and z = 2 
(thin line) for a power-law QSO source spectrum short- 
ward of the H I Lyman edge. The threshold energies required 
to ionize several commonly detected ions are indicated. Also 
shown are estimates of the UV background as constrained 
by metal absorption systems aX z = 2.73 (dashed lin e) and 
z = 1.79 (dotted line). From lAgafonova et ai\ (j2007l l. (Fig- 
ure reproduced by permission of Astronomy & Astrophysics.) 



ing from direct recombinations to the ground levels of 
H I and He ii. At higher redshifts, the UV background 
more strongly resembles the spectral shape of the sources 
because the opacity of the IGM becomes so large that 
the UV background is increasingly dominated by local 
sources. 



Combining Eq. (j40| with Eq. ([63|) , the photoionization 
rate may be expressed in terms of the emissivity as 



HI 



e-L{z)aT 
/ip(3 + aMG 



(73) 



where £l{z) is the emissivity at the Lyman edge, ajviG is 
the spectral index of the UV background just shortward 
of the Lyman edge, VeS is an effective lengthscale over 
which the distribution of sources contributes to the radi- 
ation field at any given point, the photoionization cross- 
section is approximated as Qi, ~ aT{v /vt)~^ , and Zip is 
Planck's constant. It is well known that in a uniform in- 
finite Euclidean universe with point sources of radiation 
distributed throughout, the effective lengthscale is infi- 
nite and the local radiation field diverges (Olber's para- 
dox). Cosmological effects, in particular the finite age 
of the sources and the effects of redshifting, will guaran- 
tee convergence, and r^s describes an effective horizon 
to the sources. Attenuation by a photoelectrically ab- 
sorbing IGM will shorten the effective length even fur- 
ther, and will dominate at high redshifts. These two 
extremes, a cosmologically limited lengthscale and an at- 
tenuation limited lengthscale, may be combined into the 
approximate expression r^s — 1/(1/Zh + '^/fo), where Zh 
is the lengthscale limitation due to cosmological expan- 
sion and the finite age of the sources, and tq is the length- 
scale due to attenuation by the IGM. For sources turning 
on at z = Zon with (comoving) emissivity parametrized 
as e.u{z) — 61,(1 + z)'^ {v /vl)~°'^ , the effective (proper) 
horizon to the sources is 



OLMG 



(1 + ^) 



-3/2 



l-[(l±z)/(l±Zo„)]^/'+"^-^ 



(74) 



9/2 + as - 7 

for large z (|Meiksin and Whitd.l20 03l). Values for the at- 
tenuation length tq are provided by iMeiksin and Whitg 
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(|2004i) based on numerical simulations (see § IVI.C.ip . 
The values obtained lie between the medium and high at- 
tenuation models of iMeiksin and Madaul ( 1993h derived 
from direct absorption line counts. 

Resulting values for the QSO contribution to the H i 
ionization rate are given in Table IVl at selected redshifts. 
The estimates are derived from the QSO emissivities of 
iMeiksinI (|2005[ ) and lHopkins et al\ (120071) using Eg. (1751) 
and th e attenuation lengths rp from Meiksin and Whita 
(|2004D . allowing for an extrapolation of the measured Ly- 
man Limit System distribution to z > 4 (see § IVLC.ll 
below). Results for both the Pure Luminosity Evo- 
lu tion model and the Pure Density Evolution models 
of iMeiksTnl ( 2005f ) for the QSO luminosity function are 
given. The estimate b a sed o n the QSO luminosity func- 
tion of iHopk ins et al. I (I2007D uses the emissivity shown 
in their Figure 11. All the estimates agree within the 



Star-forming galaxies are an important potential 
source of ionizing radiation. While the spectra are likely 
too soft to produce appreciable He ii- ionizing photons, 
the hot young stars in star-forming galaxies could in 
principle provide even more hydrogen-ionizing photons 
than QSOs. An estimate of the contribution of star- 
forming galaxies to the UV background requires mak- 
ing a variety of assumptions regarding the Initial Mass 
Function, populations and ages of the stars as well as 
the numbers of faint galaxies below the survey lim- 
its. Some of these uncertainties are mitigated by using 
the results of galaxy surveys that measure the intrin- 
sic flux from the galaxies close to the Lyman edge (the 
range 1300— I8OOA is typical), since the stars that pro- 
duce the ionizing radiation are the same as those dom- 
inatin g the detected light of the galaxies ( Madau et al\ . 
I1999D . Even then, however, uncertainties in the spec- 
tra of the massive stars that dominate the radiation 
at early times introduce a comparable level of uncer- 
tainty into the estimate of the co mbined ionizing emis- 



sivity from the galaxy population ( Bruzual and Chariot) . 
120031 ). Estimates for the conversion between the mea- 
sured flux at 1500 A (restframe) and the predicted flux 
shortward of the Lyman discontinuity are provided in 
Table IVTl assuming a Salpeter Initial Mass Function with 
-Slower = l-^o and Muppcr = 100 M0 and either solar 
[Z = 0.02) or 1/20 solar [Z = 0.001) me tal abundances 
based on the STARBURST99 models of lLeitherer et al\ 
(|1999D .^^ A similar val u e of f„(1500)//j912-) ~ 6 is 
quoted bv iMadau et al\ ( 1999( ). The values for the ra- 
tio /i.(1350)//i.(912-) are about 5-10% larger, and for 
/^(1700)//,.(912-) about 5-10% smaller, than the val- 
ues in Table |Vll The lower mass cut-off Miowcr = 1 Af© 
appears appropriate for some starbursts; in any case, 
choosing a lower value only alters the normalization to 



the total star formation rate as the lower mass stars con- 
tribute negligibly to the spectra except at very late times 
( Leitherer and Heckmanl . ll995h . 

Another major uncertainty is the escape fraction 
fcsc{i^) of ionizing radiation, as this must multiply the 
emissivity based on galaxy counts and finally governs the 
contribution of the star-forming galaxies to the UV back- 
ground. There are few direct measurement s of the escape 
fraction, and virtually all are upper limits. iPutman et al\ 
( 20031 ) infer a global escape fraction from the Milky Way 
of 1-2% from Ha spectroscopic measurements of nearby 
high velocity clouds. Based on a detection of Lyman con- 
tinuum flux from FUSE obse rvations of a local starburst 
galaxy, iBergvall et al\ (|2006t) suggest an escape fraction 
of 0.04 — 0.1, subject to spectral modeling uncertainties. 
iFernandez-Soto eF 'd\ (|2003., ) place an upper limit of 4% 
for galaxies in the redshift range 1.9 < z < 3.5. Similar 
up per limits based on FU SE observations were o btained 
bv iHeckman et all (l200ll) and [Peharveng et all (|2001). 
IShaplev et all ( 2006 ) report the detection of UV radia- 



tion shortward of the Lyman edge in a sample of Lyman 
Break Galaxies at z « 3 with a corresponding escape 
fraction of 14%, although the estimate is sensitive to the 
assumed calibrating population synthesis models and cor- 
rections for the IGM optical depth, and it may also apply 
only to the more luminous galaxies. 

In addition to direct observational constraints, sev- 
eral attempts have been made to estimate the es- 
cape fraction from theoretical models allowing for in- 
homogeneities in the interstellar medium of galax- 
ies, the influence of supernovae, and the flow of ion- 
izing radiatio n thro ugh chimneys produced by large 
scale outflows (ICiard i et ai, 2002; Clarke and Ocy, 2003 



Dove and Sluilll. [l994l; iDove et all \200(k iFuiita et all . 
20031 : IWood and Loebl . l2000D . Escape fraction values 



ranging from less than 1% for quiescent disks to more 
than half in an inhomogeneous IGM are obtained. The 
models suggest that high star formation rates, such as in 
Lyman Break Galaxies, may sufSciently fragment the in- 
terstellar medium of the galaxies and/or create sufficient 
channels through supernova-driven superbubbles that es- 
cape fractions in excess of 10% may be achieved, partic- 
ularly if star forming regions are distributed throughout 
the disk of a galaxy. 

A composite luminosity function of star- forming galax- 
ies detected by the Lyman-brea k method^^ at z » 3 
was derived bv iDickinso^ (|l998l ). iMadau et all (|1999D 
combined this with determinat ions of the Lyman-b reak 
population at 2.75 < z < 4 by iMadau etM (|l996l) and 
iMadau et all ( 1998t ) to estimate the emissivity of star- 
forming galaxies and their contribution to the metagalac- 
tic UV background. 



^° The Lyman edge moves into the optical range at ^ 3, per- 
mitting the identification of galaxies with large Lyman breaks, 
The code and model results a re available at as is expected for g alaxies with high star formation rates 

|http://www. stsci.edu/science/starburst997] l lLeitherer et ami999l ). 
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TABLE VI Ratio of /^(1500 A)//^(912- A) for starburst galaxies under going continuous star formation with a Salpeter stellar 
Initial Mass Function, ove r the indicated t imes for solar (Z = 0.02) and 1/20 solar (Z = 0.001) metallicity. Based on the 
STARBURST99 models of iLeitherer ef al\ ||1999D . 



z 


1 Myr 


20 Myr 


100 Myr 


300 Myr 


900 Myr 


0.02 


2.7 


4.9 


6.2 


6.7 


6.8 


0.001 


2.6 


4.6 


6.2 


6.9 


7.5 





Establishing the star- formation history of the Universe 
is currently an area of intense activity. A revised UV lu- 
mi nosity function of Ly ma n Break Galaxies was pr ovided 
by ISteidel &t a/] (|1999D . iGiavalisco et al\ (|2004t) mea- 
sured the luminosity density of Lyman Br eak Galaxies 
in the redshift interval 3.5 < z < 6.5, while [Steidel et al\ 
( 20041 ) estimated the star-formatio n history of galaxies in 
the redshift interval 1.4 < z < 2.5. lYoshida et al. (200S) 
have measured the luminosity function of Lyman Break 
Galaxies in the SUB ARU Deep Field at 2; « 4 and 5. 
iBouwens et al\ ( 20061 ) provide a comprehensive analysis 
of the UV luminosity function of galaxies at 2; « 6 and 
the evolution of the luminosity function to z w 3. 

Despite the heroic efforts to establish the UV lumi- 
nosity density and cosmic star formation rates at high 
redshifts, the contribution of the galaxies to the meta- 
galactic UV background is still uncertain due to several 
factors. In addition to the ionizing radiation escape frac- 
tion and UV-to-Lyman limit conversion factor, the effects 
of internal galactic extinction, uncertain ty in the amount 
of intergalactic extinction by the IGM ( Meiksinl . l2006af ) 
and the minimum source luminosity on the UV luminos- 
ity density and its evolution create further uncertainty in 
the net contribution of the galaxies. 

The error bud get may be illu s trated using the results 
for 3 < z < 6 ( Bouwens et al\ . l2006l : IGiavalisco et al\ . 
|2004[ ). The uncertainty in the measured luminosity 
density of the sources is 20% or larger. Another fac- 
tor of 3 uncertainty results from the u ncertain extinc- 
tion correction to the number counts ( Bouwens et al\ . 
|2006[ ). The uncertain lower luminosity limit of the 
sources introduces a factor of 2 uncertainty to the inte- 
grated contribution of the sources to the UV background 
(|Bouwens eFoZI . l2006l ) . The conversion of the flux from 



the measured bands to shortward of the Lyman edge (see 
Table IVI[) introduces at least another 20% uncertainty. 
The attenuation length r p is uncertain by about 30% 
( Meiksin and Whitd . 1200 j ) . and the uncertain evolution 
rate and birth time of the sources introduces another 10% 
uncertainty in the effective length scale to use in Eq. ([75)1 . 
Since these uncertainties are primarily systematic, the 
factors must be multiplied. Altogether, taking middle 
values of the expected ranges still leaves a combined un- 
certainty of about a factor of 4 for the contribution of the 
galaxies to Thi, excluding the uncertainty of the escape 
fraction. 

IGiavalisco et al\ ( 20041 ) obtain for the comoving lumi- 



nosity of sources at 1500A, logj^Q pisoo — 

19.25 ± 0.05 

(in units of WHz^^Mpc^^) at 3.4 < z < 4.1 and 
logioPisoo — 19.0 ± 0.1 at 4.6 < z < 5.2. Applying 
the conversion factors in Table IVII gives for the cor- 
responding emissivities shortward of the Lyman edge 
logioP9i2- - 18.4 ± 0.06 at 3.4 < z < 4.1 and 
logiQP9i2- ~ 18.2 ±0.1 at 4.6 < z < 5.2. Using these 
values in Eq. (|73D , combined wit h the a ttenuation length 
estimates of iMeiksin and Whit j (|2004l ) (see S lVI.C.ll be- 
low), gives rHi,-i2 — 2.7/csc at 3.4 < z < 4.1 and 
Thi. -12 — 2.1/esc at 4.6 < z < 5.2. Comparing with the 
required values for rHi.-i2 in Table IVl f see ? IVI.C.l() . al- 
lo wing for the con tribution from QSOs in the PLE model 
of iMeiksinI (|2005l ) (the PDE model estimates are even 
smaller), and allowing for a factor 4 uncertainty in the 
galaxy contribution imposes the restrictions on the es- 
cape fractions of 0.004 < fcsc < 0.9 at 3.4 < z < 4.1 
and 0.02 < /esc < 0.6 at 4.6 < z < 5 . 2. Sim ilar escape 
fractions are inferred by llnoue et Ml (|2006D . Galaxies 
thus plausibly make up the missing ionizing radiation, 
but the uncertainties are still too great to preclude ad- 
ditional sources, which may actually dominate. Since 
re-radiation by the IGM itsel f may contribute some 20- 
40% of the total emissivi ty ( Haardt and Madaul Il996l : 
iMeiksin and MadaullT993l ) . these restrictions on /esc may 
be further eased. The emissivities are based on source 
counts uncorrected for the effects of internal extinction 
on color. Corrections may increase the em i ssivit ies by 
as much as a factor of 3-4 ( Bouwens et al I I2OO6D . cor- 



responding to a similar reduction factor in the required 
escape fractions. 



At 5.5 < z < 6.5, iBouwens et ^ (Hooi) obtain a 
comoving integrated UV emissivity at A = 1350A of 
logio P1350 — 19.0 ± 0.25 at 5.5 < z < 6.5, with the error 
range covering the uncertainty in the lowest luminosity 
of the contributing sources. Using the same procedure 
as above, this translates into the rather large ionization 
rate of Phi, -12 — 20/csc- Correcting the colors for inter- 
nal extinction increases this value by about another 50%. 
Consistency with the residual ionization rate at z ~ 5.5, 
after subtracting the QSO contribution (which is small), 
gives 0.004 < /esc < 0.09 (or somewhat less, allowing for 
extinction corrections), a very plausible range, suggesting 
that these sources may in fact dominate the UV ionizing 
background at high redshifts. 
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FIG. 10 The probability density j(u) for the fluctuations in 
the energy density u of the UV background for a mean num- 
ber of sources within the attenuation volume of iVo = 1, 10 
and 100 for a QSO luminosity function (solid lines). The cor- 
responding asymptotic limits for ^ 1 are also shown 
(dashed lines). The peak in the distrib utions shift to lower 
ul{u) as A^o decreases. Adapted from iMeiksin and White! 
dioOS). (Fi gure reproduced by permission of the RAS and 
Blackwell Publishing.) 



B. Fluctuations in tlie UV background 

The discreteness of the sources of radiation will pro- 
duce local fluctuations in the UV background. Several 
physical effects contribute to the fluctuations. Any sin- 
gle source will dominate the radiation field sufficiently 
nearby. The result is a region of enhanced photoion- 
ization around the source and a reduction in the local 
opacity of the IGM. Known as the "proximity effect," 
the diminished strength of absorption by the Lya forest 
near an ionizing source like a QSO may be used to con- 
strain the intensity of the UV background, as is discussed 
in the next section. If the opacity of the IGM to ioniz- 
ing photons is sufficiently small that a large number of 
sources distributed over large spatial scales contribute to 
the radiation field, the fiuctuations will be small. If the 
attenuation length, which is the characteristic distance 
over which the radiation from sources is diminished by 
IGM absorption, is smaller than the characteristic clus- 
tering scale of the sources, then the radiation field will be 
modulated by the large scale spatial distribution of the 
sources. In the limit of such a short attenuation length 
that the number of sources contributing to the radiation 
field in any given region is subject to large Poisson fiuctu- 
ations, the radiation field will exhibit substan tial fiuctu- 
ations (|Meiksin and Whitd . [2005 [Zuol . Il992at) . This will 
become a particularly appreciable effect near the reion- 



ization epochs of H i and He ii, but may be an important 
effect well after the reionization epochs as well if the IGM 
opacity is sufficiently large at high redshifts. 

The magnitude of the fiuctuations may be estimated 
by considering the effect of Poisson fiuctuations in the 
source numbers on the integrated intensity. The atten- 
uation length ro is equivalent to the mean free path of 
an ionizing photon, the distance over which the mean 
number of absorbers encountered is unity. The probabil- 
ity that a photon survives a journey of a distance r is 
then exp(— r/ro). In general the distance will depend on 
frequency, but most of the absorption takes place near 
the photoelectric edge at frequency vt- Provided the 
distances are sufficiently small that cosmological effects 
may be neglected, the generation of the radiation field 
may be considered in the Euclidean limit. Any individual 
source of luminosity a distance r from a given loca- 
tion will then contribute a flux exp(— r/ro)/47rr^ for 
V > Vt- Allowing for a luminosity distribution function 
$(Ly) per unit luminosity per unit volume with mean 
spatial density of sources between minimum and maxi- 
mum luminosities Lmin and Lmax oi n — J dL^ ^(L^), 
and defining Nq = (47r/3)rQn to be the mean number of 
sources within an attenuation volume, the distribution 
of the dimensionless energy density u = u^/u*, where 
u* = ( £*/c)/47rr'n for a cha racteristic luminosity L*, be- 
comes (jMeiksin and Whitd . [200al 
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where x = L^/L* and (/>(a;) is the normalized luminosity 
function of the sources 
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The function 



(77) 



with t{p) defined implicitly through p(r) = r^^e^^. 
For large u, the distribution takes the asymptotic form 
f{u) ^ (3/4)iVo(x'^/^)u~^/^, where x^^"^ is averaged over 
4>{x). It follows that the variance (u^) — (u)^ of the radi- 

1/2 

ation energy density is undefined, diverging like Umax for 

Umax 1 . 

Examples of the resulting distribution function are 
shown in Fig. [TOl A QSO luminosity function of the 
form Eq. ()72|) was used with parameters given by the 
maximum likeli hood Pure Luminosity Evolution fit of 
( Meiksinl . |2005| ) at z = 3 (renormalized to Hq = 
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70kms-iMpc-i): pi = 1.24, /?2 = 2.70, (/.;,„p„. = 
2.6 X IQ-^Mpc-^mag-i, and M* = -27.2 at restframe 
A = 1450A. In anticipation of the discussion of hehum 
reionization below (§ IIV.C.2p . the attenuation lengths 
were fixed at 7, 15 and 33 (proper) Mpc, correspond- 
ing to averages of 1, 10 and 100 sources per attenuation 
volume. When the average is small, the distribution be- 
comes quite skewed towards low values, corresponding to 
few or no sources within an attenuation volume. 

The Poisson fluctuations of sources within an atten- 
uation volume will also give rise to spatial correla tions 
in the radiation field. The correlation function is 
ll992bD 



become 



0(r) 



3No (x)^ r 



du — In 

r/ro ^ 



u — r /tq 

. (78) 

where the averages (. . .) are performed over the dimen- 
sionless QSO luminosity function 0(a;) (eq. [75]). For 
the luminosity function above at z = 3, the coefficient 
(a;^)/3(x)^ ~ 2.6 x 10"*, so that substantial correlations 
are expected on scales r < tq. 



C. Observational constraints on the UV background 

1. The proximity effect 

The number of Lya fo rest systems decreases systemat- 
ically near brig ht QSO s (|Lu et aLl . [T99lt [Murdoch et all 



ly near on 
1 iTvtleil . 



1987bl ). The decrease, known as the 



1198 

"proximity effect," is expected sufficiently near a QSO 
where its ioniz i ng flu x dominates th e UV background 
(|Carswell et all Il987t) . iBaitlik et all (|l988l ) quantified 
the effect in terms of the statistical properties of the 
Lya forest to obtain an estimate of the UV background. 

The factor by which the H i column density in an ab- 
sorption system will be reduced is 1 -I- w, where lo is the 
ratio of the ionization rate F^j due to the QSO at a lu- 
minosity distance between the QSO and the absorber 
to the metagalactic ionization rate Fhi. The ionization 
rate due to the QSO is 



HI 



dv 

hvv Airr'i 



1 



1 + zq rl J hpv 



dl 
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where f j, is the observed flux from the QSO at frequency 
1^' = + Zq) and dL = {1 + zq)c J^'^ dz/H{z) is the 
luminosity distance to the QSO a t redshift zq for a flat 
universe [riK = 0) (lPeeblesl . [l993l) . For a power-law H i 
column density distribution dN/dNni oc N^f , the in- 
crease in intensity will decrease the number of absorp- 
tion systems above a given threshold N^^'^^^^ by the fac- 
tor (A^^hresh)l-/3/[^throsh/(j ^^)]l-;3 ^ (1 + ^)-(/3-l) . 

A mean line density dN/dz = iVo(l + z)'^ would then 



dN 

dz 



No(l + zy[l + Lo{z)]-^^-^\ 



(80) 



Measurements of a; as a function of distance from 
the QSO as inferred from the depression in line density 
dN/dz permits a determination of the metagalactic Fhi. 

Since iBaitlik et all i 1988f) . several estimates of the 
metagalactic ion ization rate have been made using the 
proximity effect. IScott et al\ ( 2000l ) performed a uniform 
analysis of the spectra of 99 QSOs at moderate resolution 
(~ lA) over the redshift range < z < 5 to obtain an 



estimate of Fhi = 1-9 



+ 1.2 



X 10 



-12 , 



"1 over 1.7 < z < 3.8. 



(They also provide references to much of the previous 
literature.) The rates have generally exceeded the esti- 
mates of the contribution from quasars. Several system- 
atic effects could account for some of the discrepancy. 
Redshifts for the quasars measured from low ionization 
permitted lines like the Balmer lines or Mg ii or from 
forbidden hues like [O iii] AA4959, 5007A are redshifted 
compared w ith Lya and C iv by a few to several hun- 
dred kms~^ (iBoroson and Green . ll992HLaor e7aH . ll995l: 
IVanden Berk et all 2001 ). If the higher redshifts reflect 
the true distances to the QSOs, then the QSO flux pro- 
ducing the proximity effect will be overestimated using 
QSO redshifts based on Lya or C iv, which in turn would 
result in an overestimate of the met agalactic ionizatio n 
rate by as much as a factor of two ( Scott et~cil\ . |2000| ). 
Numerical simulations suggest that enhanced large-scale 
structure in the vicinity of a QSO will partly compen- 
sate for the photoionizing influence of the QSO, again 
resulting in an overestim ate of the metagalactic rate by a s 
much as a factor of 2-3 ( Faucher-Giguere et all [20083) • 



Quasar variability will introduce scatter in the esti- 
mates for Fhi and degrade the fit of to ag ainst for 
a mo del assuming a constant quasar flux (iBaitlik et all 
Il988l ). Episodic behavior of QSO luminosities could in- 
troduce a further systematic in the estimates of the meta- 
galactic ionization rate. If a QSO satisfies a survey mag- 
nitude threshold as a result of entering a bright phase, the 
neutral fraction will lag behind, resulting in an overesti- 
mate of the metagalactic rate. From Eq. (|45| for 0^1, 
it follows that the actual metagalactic rate is a factor 
1 — {1 + uj)/ exp(— F^jt) smaller than the estimate as- 
suming equilibrium, where to is measured from the deficit 
in Lya forest absorption lines near the QSO and F^j* — 
Fhi-I-F^j. Time-averaging over a period tq during which 
the QSO is sufficiently bright to enter the survey pro- 
duces the average factor 1— (FHiiQ)"^[l — exp(— F^jtg)]. 
For tg « 10^^ — 10^"^ s, this would result in too low a 
value for u, and so too little suppression of the Lya forest 
line density ove r a distance ctp fr om the QSO, as is fre- 
quently found ( Scott et~ai\ . |200(][ ). An overestimate of 
the metagalactic photoionization rate by several percent 
could result. 
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FIG. 11 Distribution of column density ratio 77 — Nhcu/Nui, 
for A^'hcII measured from FUSE observations and correspond- 
ing A^Hi from Keck observations for absorption features over 
the redshift range 2.30 < z < 2.75 from Fechner et al. (2006). 
Systems detected in spectral bins affected by metal line ab- 
sorption have been excluded. The systems are further re- 
stricted to lie in the range 0.01 < ro,Hi < 0.1 to assure com- 
pleteness and reliable fits. Without this latter restriction, the 
full distribution ranges over —0.2 < logj^Q < 3. The solid 
curve shows the expected 77 distribution allowing for fluctu- 
ations in the He ii-ionizing background for an average of 10 
QSOs per attenuation volume. The dashed line allows for a 
small increase in the mean 77 values within the measurement 
errors. 



2. Constraints from helium ionization 

The ratio of He 11 to H i column densities places a con- 
straint on the ratio of the metagalactic UV background 
at the He 11 and H i photoelectric edges. From Eq. ([1^ 
and the corresponding equation for He 11, assuming the 
helium (like the hydrogen) is highly ionized and so almost 
entirely in the form of He iii, 



_ ^HcII _ Thi ^Hc aYlclll{T) 

Nm FhoII n-a am\{T) 



(81) 



where Fhi and rfjeii are the photoionization rates per 
atom of H I and He 11, respectively. For an approxi- 
mate temperature dependence of anii c>c y~o 'i'2 Yieai 
T = 10^ K (and similarly for He 11) and a hydrogen- 
to-helium number density ratio of 12.9, it follows from 
Eq. ([STjl . Eq.((40|) and the analogous equation for FhcII 
that 
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3 + aL.Hcii 
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Sl ^ 0.42 



F 



HI 



HI 



HcII 



(82) 



where the "softness parameter" Sl — wl.hi/wl^hoH is the 
ratio of the metagalacti c UV energy density at th e H i 
and He 11 Lyman edges ( Madau and Meiksinl . Il994[ ) , and 
ctL.m and aL,Heii are the spectral indices just shortward 
of the respective edges. (Allowing for the full tempera- 



ture dependences of the H 11 and He iii radiative recombi- 
nation rates i ntroduces a weak temperature dependence 
of 77 cx T" 055; IShull et al\ (|2004[ ).) Because of the steep- 
ness of Ui, at energies just above 4 Ry (Fig. [S]) , instead 
of the softness parameter, it is usually more convenient 
to introduce the ionization rate ratio ^' = Fhi/FhcII, as 
it is this ratio which governs 77 directly. 

It follows from Eq. that the ratio of the He 11 to 
H I optical depths at line center is given by 



To.HcII 



1 b 



HI 



To,m 4 6hcII 



77 : 



0.10-^*. 

OHoII 



(83) 



The value of is predicted to be large for a UV back- 
ground dominated by QSOs (Fig. For an intrin- 
sic spectral index of aq = 1.5, ^ « 100 at z » 3 , 
corresponding to 77 4 ()Haardt and Madaul . Il996l : 
iMadau and Meiksinl . Il994l ). while for gp = 2, ^ ^ 340 , 
corresponding to 77 « 140 ( Madau and Meiksinl . Il994[ ). 
Interpolating to aq = 1.8 gives ^ « 220 and 77 « 90. Ad- 
ditional photoionizing sources, like galaxies, will increase 
^ and 77 further. The ratio of optical depths is thus ex- 
pected to be large, so that a strong He 11 Lya signature 
is expected from the IGM. The ratio also depends on the 
ratio of the Dopplcr parameters. For pure thermal broad- 
ening, it follows from Eq. pO|) that feHi/^Hoii = 2, and 
To,Hcii/7"o,Hi ~ 0.2^* will be quite large. But if velocity- 
broadening dominates (see Eq. [TT]), then 6hi/&HcII ~ 1 
and to,hoIi/''o.hi ~ Cl^*. Systems optically thin in H i 
Lya with tq^hi as low as 0.1 will form saturated He 11 
Lya lines. 

Because of line-blanketing, the measured optical 
depths will scale according to Eq. (|3ip . A large 
effective optical de pth d ue to blanketing i s ex- 
pected (iGiroux et all \l99A iMadau and Meiksinl . 1 1994) : 
11993^ 



iMiralda-Escude . 
by observations 



(see 



This 

Hlb^I) . 



has been confirmed 

^ The estimate of 

IMadau and Meiksinl (|1994 requires y > 200 for t h e 95% 
lower limit thoII > 1.5 at z = 3.2 of lHogan et al. 1 (I1997D 
if the absorbers are thermally broadened, but 4" > 1000 
to reach the 95% upper hmit of thcH < 3.0. For velocity- 
broadened absorbers, however, for which IZheng et al\ 
(^04b) find evidence, the upper limit may be satisfied 
by ^ < 500. If the lower hmit thcII > 4.8 found by 
iHeap et all ( 2000l ) apphes, however, then 4" > 1000 is 
required even for velocity broadened absorbers. 

There is additional reason to believe values of 4* > 
1000 occur at z > 3. While the detection of a high He 11 
Lya optical depth was expected, the wide range in fluc- 
tuations in 77 came as a surprise. Some of the systems 
indicate values as large as 77 > 1000 a,t based on 

follow-up optical observations of the H i L ya systems cor- 
responding to the He 11 systems detected (iFechner et all 
[2OO6; Rci mers et al\ . [2005 -. Zhen g et al\ . l2004bD Com- 
bining A^Heii measurements from FL^S'-E' observations over 
the redshift range 2.30 < z < 2.75 with measurements 
of A^Hi for t he correspond i ng H i systems using Keck 
observation, iFechner et al\ ( 20061 ) obtain a very broad 
distribution for 77, as shown in Fig. 1111 Systems de- 
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tected in spectral bins affected by metal line absorption 
have been excluded. The systems are further restricted 
to those with line center optical depths in the range 
0.01 < To, HI < 0.1 to assure completeness and reliable 
fits. Without this latter restriction, the full distribution 
ranges over —0.2 < log^Q 77 < 3. A tren d of decreasini 
n with decreasing rcdshift is also found (jFechner et 
[2OO6; Zheng et ai, 2004b). Only fluctuations producing 
an extremely soft spectrum with ^E* > 2000 could account 
for the largest 77 excursions. Such high values could arise 
from strong local H i-ionizing sources, or from a broad 
distribution in the He ii-ionizing background with occa- 
sional low excursions. Higher valu es of 7? are found in 
H I voids, regions with thi < 0.05 (jFechner et hood 
IShuU et ali^OM), again consistent with a local overden- 
sity of hydrogen-ionizing sources compared with He 11- 
ionizing sources producing the H i voids. The trend is 
also consistent with an increase in the relative attenua- 
tion length of hydrogen-ionizing photons compared with 
He ii-ionizing photons, so that more distant hydrogen- 
ionizing (but not He ii-ionizing) sources are able to add 
to the radiation field in the H i voids. 

Monte Carlo realizations of the expect e d fiuc tuations 
in THeii were performed bv iFardal et all (1991), allow- 
ing for the discreteness of the absorption systems and 
of the quasars. At z = 3, they find a range of values 
1-3^ THeii ^ 2.3, peaking at thoii ~ 1.9. Since the He 11 
column density scales inversely with the local intensity of 
the UV background (provided most of the helium is ion- 
ized to He III), the contribution of source discreteness to 
the distribution in the He 11 to H i column density ratio f] 
may be predicted from Eq. ([75)1 as f{ri) = f{u)\du/dri\ = 
f{u)u'^/[T]{{u))(u)], where 77 = ■r]{{u)){u) /u.^'^ The flat- 
tening of f{u) for 7i <C 1 shows that the average of 
77 oc 1/u is undefined, diverging logarithmically at the 
upper end due to regions of low ionization. 

Using the distribution shown in Fig. [TU] for an average 
of A^o = 10 QSO sources per attenuation volume, the re- 
sulting distribution for 77 is shown in Fig. 1111 The distri- 
bution function shown by the solid line is computed as the 
average of two distributions matched to the measured dis- 
tributions with mean and spread of logj^g V — ^-70 ± 0.32 
at -g ~ 2.33 and lo gm?? = 1.98 ± 0.32 at z ~ 2.68 
( Fechner et Ml, I20O6D . The dashed line uses means in- 
creased by 0.25(7 to logiQ 77 = 1.78 and 2.06. The mea- 



Because the H l-ionizing radiation field will be much more nar- 
rowly distributed, its fluctuations will contribute much less to 
the fluctuations in r] and so may be neglected, except for large 
values due to the proximity effect. The joint H I and He II prox- 
imity effect would need to be modeled separately, including the 
effect on changing the attenuation length, which may be espe- 
cially important for He ll-ionizing photons when the attenuation 
length and range of the proximity effect are comparable. 
The presence of diffuse radiation will limit the maximum values 
of rj, as will the approach to near full recombination to He II. The 
maximum value, however, will still be so large that no average 
will be well-determined. 
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FIG. 12 Comparison of the observed values of Si iv/C iv vs 
C ii/C IV with model predictions using CLOUDY with an in- 
put power law spectrum. Various breaks are allowed 
across the He 11 Lyman edge: dashed line - no break; dot- 
ted line - factor of 2; dash-dotted line - factor of 10; long 
dashed line - factor of 100; solid line - factor of 1000. The 
left panel shows results for a solar abundance ratio of Si/C; 
the right panel uses a ratio of Si/C 3 times solar. The dif- 
ferently shaped points indicate different lines of sight. Small 
point s indicate 2: < 3.1 absorbers , while large indicate z > 3.1. 
From lSongaila and Cowi3 (|1996| ). (Figure reproduced by per- 
mission of the A AS.) 



sured values are well described by the model, except 
for a small tail at logj^Q 77 ~ 2.6. These values may be 
produced by patches of only partially ionized helium, 
as would remain duri ng the epoch of He 11 reionization 
before its completion ( Reimers et aLl , [2005HZheng et all . 
j2004b). Alternatively, they may result from local H i- 
ionizing sources reducing the H i column density, but 
not the He 11. Effects of large scale structure and ra- 
diative transfer on the fluc tuations were explo r ed us - 
ing numerica l simulation s by Mase lli and Fcrrara (2005^), 
iBolton et all (|2006>) and lTittlev and Meiksin. (,2001) (see 
g IVLC. 31 below). 

The measured values for rj may be compared with 
the expectation for a UV background dominated by 
QSOs. By Eq. ((82)) . the mean values for 77 corre- 
spond to 140 < < 270. From Eq. * ~ 
(e™/eL'")(aT (^SiAfff")- Approximating the 



emissivity as oc v , taking ~ 70 Mpc from 



Eq. (fT05)) below and 



^HcII 



15 Mpc for 10 QSOs in 



an attenuation volume at the He 11 Lyman edge f ^IV.Bp . 
gives 5* ~ 150 — 300 for 1.5 < as < 2.0, producing a good 
match to the required average value of 5" for a plausible 
emissivity spectral shape due to QSO sources. 



3. Constraints from metal ionization 

The presence of metal absorption systems permits fur- 
ther constraints on the shape of the UV metagalactic 
background. The ionization edges for several transitions 
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of measured metal ions, primarily carbon, nitrogen, oxy- 
gen and silicon, lie in the energy range 1-10 Ry, as shown 
in Fig. [H The measured column density ratios of ion- 
ized states of these elements have been used to distin- 
guish a UV background dominated by hard sources like 
QSOs or Active Galactic Nuclei ( AGN) from a back- 
ground dominated by galaxies/stars (|Madaul . ll99lh . The 
inferred ionization fractions may also be used to con- 
vert the measured column densities of the detected ions 
compared with the H i column densit y into physica l 
densities, sizes and total metallicities (iBergeron etal. , 
1994 ; [Bergeron and Stasinskal . Il986t iDonahue and Shul] , 
19911 ). The findings suggest Lya forest systems have 



metallicities less than 0.01 solar, while Lyman Limit 
Systems have metallicities from less than 0.01 solar up 
to 0.1 solar. Some systems show evidence for multiple 
phases with densities different by at least a factor of 10 
(|Bergeron et al\ . \l994 ]. Sizes of typically tens of kilopar- 
secs are found, with total hydrogen densities of roughly 
10^ - lO'^m"^. 

Based on a compa rison of Si iv to C iv vs C ii to C iv, 
ISongaila and Cowid ([l996.) find systematically different 
behaviors for systems at z < 3.1 compared with z > 3.1, 
as shown in Fig. [121 The Si iv to C iv ratio for a given 
C IV to C II ratio cannot be satisfied for any softness pa- 
rameter without enhancing the Si/C ratio to about 3. In 
this case, the high redshift systems require a break of a 
factor of 10 to 100 at the He ii Lyman edge, while little or 
no break is required by the lower redshift systems. The 
large break at z > 3 is similar to the finding based on the 
He II to H I ratios found at z ~ 3 discussed above. Re- 
fined estimates of the UV background, including a plau- 
sible contribution from starburst gal axies, suggest the 
Si/C ratio need not exceed twice solar ( Giroux and Shulll . 
1997). Based on obse rvations of metal systems at z > 3, 
Savaglio et d\ (|1997[ ) and D'Od orico et al. ( 1998a) sim- 



ilarly find a much larger ratio of Si iv to C iv column 
densities compared with the expectation for a QSO or 
AGN dominated UV background. They find a large soft- 
ness parameter between the He ii and H i Lyman edges 
of Sl > 1000 is required to bring the overabundance of 
silicon over carbon to below a factor of 10 over solar. 

Using 19 high quality QS O spectra oy e r the red- 
shift range 1 . 5 < z < 4.5, ISchave et~al\ (|2003l ) and 
lAguirre et al\ ( 20041 ) perform a pixel optical depth anal- 
ysis of C IV, Si IV and H i in conjunction with simu- 
lated spectra. They find evidence for a highly inho- 
mogeneous distribution of carbon with a cosmic mean 
of [C/H] = -2.80 ± 0.13 at z ~ 3 and no evi- 
dence for evolution over the redshift range 1.8 < z < 
4. 1, adopting the ga l axy-Q SO metagalactic spectrum 
of iHaardt and Madaul (|2001i ). rescaled to reproduce the 
measured mean Lya transmitted flux. Introducing a 
large break in the metagalactic background at energies 
above 4 Ry, the He ii photoelectric threshold, for all 
redshifts indicates a decrease in metallicity with time, 
which could only be physical if large quantities of metal- 
poor gas were accreting onto the metal-enriched absorp- 



tion systems or replacing them in the surveys as they 
accreted onto galaxy halos, possibilities that perhaps 
should not be discounted. Introducing the large break 
only at z > 3.2, as would occur if He ii were ionized 
late, however, yields a near constant (at most a weakly 
increasing with time) carbon abundance. They find no 
redshift evolution in the ratio rsiiv/TciVi and the silicon 
and carbon data require no sharp spectral change in the 
UV background at 2: ~ 3. They find a strong correlation 
between rsiiv/To Tv and totv, and argue the redshift ef- 
fect reported bv ISongaila and Cowid ( ~996h may result 
from a different range in rciv probed as a function of 
redshift, possibly due to different quality spectra for sam- 
ples at z < 3 vs 2 > 3. They also infer a Si/C ratio of 
[Si/C] = 0.77 ± 0.05, with no indication for evolution in 
the ratio. 

In a study of metal absorption systems in the spectra of 
z ~ 2.5 QSOs with i ntervening C iv and O vi detections, 
ISimcoe et all ^oM) find roughly 30% of the Lya forest 
is enriched to abundances below [C, 0/H] < — 3.5. The 
systems are adequately modeled using a UV background 
dominated by QSOs or AGN having a power-law source 
spectrum h'~^^^, assuming a solar abundance ratio of car- 
bon to oxygen. If instead the absorption systems are 
overabundant in oxygen by [C/0] = —0.5, as found for 
metal-poor Galactic halo stars, then it is necessary to 
include a galactic radiation component. 

Modeling metal absorption systems at z ~ 2 using 
CLOUDY and an assumed UV metagalactic background 
with contributions f rom both QSO sources and galaxies, 
ISimcoe et all ( |2006( ) infer a trend of decreasing sizes and 
increasing densities from high to low ionization states. 
The largest, lowest density systems, with sizes of 0.1- 
10 kpc and densities of 2.5^ logj^g "-h ^; 3.5 (in units of 
m~'^), are seen exclusively in C iv or Si iv . These sizes 
are comparable to the minimum inferred sizes based on 
coincident absorption systems along neighboring lines of 
sight (§ III.B.6p . The lower ionization species C 11 and 
Si II have sizes of 0.1-1 kpc and densities of n ~ 10* m~'^, 
while Mg II and Al 11 systems have inferred sizes of about 
100 pc and densities of nu ~ 10^ m~^. The smallest and 
densest structures are Fe 11 absorbers, with sizes of about 
1 pc and densities of nn — 10^ m~'^. 

Comprehensive analyses of metal absorption systems 
between 1.4 < z < 3 that are optically thin at t he H i 
Lyman edge were car ried out by Agafonova et al. ( 20051 ). 
iReimers et d\ (|2006l) and I Agafonova et all l^20m . They 
employ a Monte Carlo Inversion procedure which allows 
for a random distribution of gas density and velocity 
along the line of sight. The procedure assumes a constant 
metallicity for each system, that the gas is optically thin 
to the ionizing UV radiation, assumed uniform, and that 
the gas is in thermal and ionization equilibrium. The 
ionizing UV background is expressed in parametric form. 
All parameters of the model are then varied to obtain the 
best fit to the metal absorption features, using CLOUDY 
to compute the ionization state s and temperatur e s. Th e 
method is described in detail bv lLevshakov et'al. I (|2000D . 
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lAgafonova et all ( 20051 ) and references therein. Applied 
to individual systems, a range in He ii to H i column den - 
sities is inferred from the data (jAgafonova et 
suggestive of a variable UV radiation field. Since the 
variability is generally based on systems optically thick 
to He ii-ionizing radiation, the method does not distin- 
guish between a UV radiation field that is variable on 
the metagalactic scale external to the absorbers and a 
radiation field that differs between the absorbers as a 
consequence of internal absorption. 

These authors find a range of hydrogen densities from 
the analyses, with values as low as nn < 10^ m^'^ to 
as high as rtH > 10^ m~'^, similar to the findings of 
ISimcoe et al\ (I2006D . The line-of-sight sizes of the metal 
absorbers range from a few parsecs to over 100 kpc. 
Systems at z ^ 3 are best fit by imposing a strong re- 
duction in flux between 3 Ry and 4 Ry, as shown in 
Fig. [HI Such a reduction could be accounted for by 
strong He ii Lya absorption. At z < 1.8, a much shal- 
lo wer break at 4 Ry is r equire d compared with the model 
of iHaardt and Madaul U^M), 

suggesting there may be 
fewer high column density systems filtering the He ii ion- 
izing radiation than assumed by the model. The ra- 
diation field at z < 3 is dominated by sources with 
QSO/AGN spectra. Only a very small soft spectrum 
component from galaxies may contribute, suggesting an 
escape fraction of fcsc < 0.05. Large fiuctuations in the 
UV background are also required. At z^3, these are 
likely produced predominantly by the effects of radiative 
transfer through an IGM optically thick to He ii-ionizing 
photons. The optical depth of He ii-ionizing photons is 
dominated by systems with 15 < logj^Q A'hi < 16, while 
much rarer systems with logj^g TVhi > 17.2 dominate the 
H i-ionizing photon optical depth. As a result, the mean 
free path of He ii-ionizing photons is much shorter than 
for H i-ionizing photons, and relatively fewer sources con- 
tribute to the He ii-ionizing radiation background in any 
given region. By 2, however, log^g A^hi < 16 systems 
are sufficiently sparse that the IGM should become op- 
tically thin to He ii- ionizing radiation to large distances, 
and the fluctuations are instead likely produced predom- 
inantly by the spread in the spectral indices of quasars. 
Fluctuations in = -A^hcIi/-^hi and an anti-correlation 
between r/ and TVhi may be accounted for by the filtering 
through the IGM of radiation from QSO sources with a 
variety of spectral shapes and by differences in the mean 
densities of the absorbing clouds. The progressive hard- 
ening of the UV background towards lower z would also 
account fo r the decrease in the median values for r i as z 
decreases (|Fechner et all 120061: IZheng et all l2004bl ) . 



The modeling has some limitations, especially for sys- 
tems with large values of rj. A system with logj^Q A^hi < 
17.2 is optically thin to H i-ionizing radiation but thick 
to He ii-ionizing radiation, so that a break at 3-4 Ry 
must be introduced to account for the radiative transfer 
internal to the system. The assumption of thermal equi- 
librium assumes adiabatic cooling or heating are negligi- 
ble. The timescales are also long for the more rarefied 



systems. For a temperature of T « 20000 K, the cool- 
ing time due to radiative recombination to H i exceeds 
10^ yr for nni < 200 m^'^, which is becoming long, par- 
ticularly if flows occur on shorter timescale s so that adi- 
abatic gains or losses are not negligible. ISimcoe et all 
( 20041 ) also point out that the O vi abundances are sen- 
sitive to the poorly determined soft x-ray background. 
Still, the studies have lent considerable insight to the 
physical nature of the metal absorption systems and the 
background UV radiation fleld that ionizes them. 



V. EARLY MODELS 

Several models have been suggested for describing the 
physical structure of the Lya forest absorption systems. 
All of the early models have been supplanted by numeri- 
cal cosmological simulations, which provide a largely suc- 
cessful description of the Lya forest and associated ab- 
sorption systems within the broader context of cosmo- 
logical structure formation. Vestiges of the early models, 
however, survive. In this section, elements of those mod- 
els are discussed that continue to provide physical insight 
into the findings of the simulations, or which may still 
offer viable alternative descriptions of some absorption 
systems. 



A. Pressure-confined clouds 

The earliest physical model for the absorbers was 
one in which the absorption features aris e from distinct 
cloud s pressure-confined by a hot IGM ( Sargent et all 
Il980l ). Denoting the cloud temperature and hydrogen 
number density by Tc and Uc, and those of the IGM 
by TiGM and tiigm, the condition of pressure equilib- 
rium is UcTc = jt-igmTigm- Allowing for clouds photoion- 
ized by a QSO/AGN-dominated UV background, small 
enough to be Jeans stable and able to withstand ther- 
mal evaporation by the hot confining medium, while re- 
quiring the IGM to meet the Gunn-Peterson constraint, 
ISargent et all (|l980( ) found for the cloud parameters at 
z « 2.4, ric = 100 m-3^ 3 x 10^ K, njGM = 10 m^^ 

and TiGM = 3 x 10^ K. These values were somewhat 
modified bv lOstriker and Ikeuchil (|l983[ ) and others. 

The pressure-confined models suffered from a serious 
shortcoming: a wide range of sizes was required to repro- 
duce the range in observed H i column densities. The 
higher column density systems would need to be so large 
that they would become unstable under their own grav- 
ity. Jeans stability limits the size of a self-gravitating 
cloud of density ric and temperature Tc to be sufficiently 
small that a sound wave can cross it in less than the time 
for it to gravitationally collapse. This requirement is nec- 
essary to ensure the internal gas pressure gradient of the 
cloud is able to balance its weight. For an isothermal gas. 
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this corresponds to a maximum size of 



A , = 



1/2 



f '^igmTigm 



5.2 X 10- ( ^^i?^^ ( i^^l m(84) 



In ionization equilibrium with an ionization rate Fhi ~ 
10~^^s~^, Eq. (j44|) gives a neutral hydrogen density for 
a cloud at Tc = 3 X 10"^ K of nm ^ 0.002 m^^. This gives 
a maximum column density for a Jeans stable cloud of 
-^Hi.max — 10^'^ cm^^. The observed distribution extends 
to values two orders of magnitude larger before the ab- 
sorbers become Lyman Limit Systems. A further failure 
of the model is that a system with a column density of 
Nui = 10^''cm~^ would have a typical size of 17 kpc. 
This conflicts with the sizes inferred from coincident ab- 
sorbers along neighboring lines of sight to Icnscd QSO 
images, which give cloud sizes in excess of 70 kpc for 
systems with Nm < lO^cm'^ fq ilLB.Sp . 

Elements of the pressure-confined model, however, 
may still describe some absorption systems. The clouds 
could, for example, be confined by a wide range of 
pressures in different environ ments, which are them- 
selves gravitationally confined (iBaron et a/.l . Il989l ). The 
association of metal absorption systems with galax- 
ies suggests they may be pressure-confined systems 
within galactic halos, as may be Lyman Limit Systems 
(|Mo and Miralda-Escud"i.[l996.) . 



B. Dark matter minihalos 

A more viable physical model for the absorbers is grav- 
itational confinement by a dark matter halo. iReed ( 19861 ) 
recognized that the Lya forest is a necessary consequence 
of the Cold Dark Matter (CDM) scenario for structure 
formation. The same normalization of the CDM primor- 
dial matter power spectrum that recovers the clustering 
of galaxies predicts an abundance of low mass minihalos 
with virial velocities comparable to the sound speed of 
photoionized gas. The gas trapped in these halos will 
produce absorption features with H i column densities 
and line w i dths c omparable to those observed ( Ikeuchil . 
119861 : IReesl [l986l) . 

To fix parameters, consider the formation of a mini- 
halo through gravitational collapse. The formation of a 
dark matter halo may be approximated by the spheri- 
cal collapse model. The evolution of the radius r{t) of 
a spherical shell of mass M in an Einstein-deSitter uni- 
verse, which provides a good descript ion at high red shifts, 
is given by the parametric solution (|Peeblesl ll980t) 



r = r\/(l — COS0), t 



rv 



(6* -sin 61), 



(85) 



where ry is the virial radius of the shell and Uciic = 
(GAI/rvY^^ is the circular velocity of the shell at the 



time of virialization. For a top hat initial spherical per- 
turbation, the central dark matter density at the time 
of virialization is poM = (9/8)(37r -I- 2)^pcrit(2;y), where 
Pcrit(-Zy) is the critical Einstein-deSitter density at the 
virialization epoch at redshift zy- Since the gas is dis- 
sipative, its core radius will be smaller than the virial 
radius of the dark matter halo. Approximating the gas 
as being in hydrostatic equilibrium within a uniform den- 
sity dark matter halo, the gas will have a gaussian profile 
p{r) = p(0) exp(— r^/r^) with core radius 



/ k^T/m 
V27rGpDM/3 



1/2 



13(1 + Zv)-3/2T4/2/i-1 kpc. 



(86) 

where T4 is the gas temperature in units of 10^ K. A 
minihalo virialized at zy = 3 will be only a few kilo- 
parsecs in size. A minihalo of mass 10^ Mq would have 
an adequate circular velocity to confine the gas. The 
H I column densities produced will range roughly over 
13 < logiQ iVni < 16. Higher column den sity systems be - 
come self- gravitating and Jeans unstable ( Meiksinl . ll994l ). 
and will likely produce stars. 



iBond et al\ ( 19881 ) point out that the gas will in gen- 
eral be in a dynamical state, either collapsing with an 
increasing core density or expanding after photoioniza- 
tion. The observational signatures and the r mal p roper- 
ties of minihalos were explored bv iMeiksinI ( 1994 ). Gas 
associated with the minihalos gives rise to absorption 
features with iVni $ 10^'^ cm~^ on scales of 100 kpc at 
2; = 3. Even though the gas is in motion, the line profiles 
are found to remain nearly Doppler in shape, although 
broadened by the gas motion. In the idealized case of an 
isothermal gas in free expansion, so that the gas profile 
is gaussian with a linear velocity profile v = Vc{r/rc), 
the Voigt profile is exactly preserved for lines of sight 
passing through the cloud, but with the Doppler param- 
eter broadened to 6 = {2kBT / mn + v^Y^^ . Typical tem- 
peratures in the inner region of minihalos with circular 
velocities of Vc — 40 km s~^, which give rise to systems 
with Nm > 10^5 cm-2, ar e 30 - 40 x 10"^ K consistent 
with the measurements of lRauch et al. I (|l996h for metal 
absorption systems expected to be associated with com- 
parable H I column densities (see § III. B. 61 above) . In the 
inner 10-30 kpc, temperatures of T ~ 30 — 60 x 10^ K are 
predicted, with substantial deviations from thermal equi- 
librium due to adiabatic heating of the infalling gas. In 
the outer 10-100 kpc, the gas is significantly adiabatically 
cooled by Hubble expansion. Temperatures half these 
values are expected in smaller halos of Vc — 20 km s~^. 

Further consequences of the minihalo model, in- 
cluding number counts, effects of photoionization, and 
ablation due to ram pressure on minihalos moving 
relat ive to their environments, have been explored 
bv llkeuchi aLl ( 1989a. bh. iMurakami and Ikeuchil 



(I1990D. iMiralda-Escude and Reed 
iMurakami and Ikeuchil (|19931I1994[k 



(|1993D . 



and 



37 



C. Caustics and sheets 

A gravitational perturbation will generically form 
sheets since the enhanced density along a given direc- 
tion will accelerate the collapse along that direction. 
Once a growing density fluctuation becomes non-linear, 
it will produce a caustic. The caustic formation is well- 
described by the Zel'dovich approximation ((Zel'dovich . 
Il97nh 

x(q, t) = q - i?(t)V*(q), v(q, t) = -aI)(t)V*(q), 

(87) 

where x is the comoving position of a particle with ini- 
tial comoving position q, v = ax is the peculiar velocity, 
^'(q) describes the initial deformation of the density field, 
the gradient is with respect to q, and D{t) is the growth 
factor for line ar perturbations D{t) = (a/a) / da/ a? 
( Peeblesl . Il993l ). The physical radius is related to the co- 
moving radius through r = a(<)x, where a{t) = 1/(1 + z) 
is the expansion factor of the Universe at the epoch cor- 
responding to redshift z. Conservation of mass gives 
for the density, up until the time of caustic formation, 
p{x,t) — p{q)/\cfix/d'^q\, where the denominator is the 
determinant of the Jacobian of the coordinate transfor- 
mation. For ID collapse, the density prior to caustic for- 
mation grows as p{x,t) = pQa{t)~^/[l - D{t)d^'^ / dq^], 
which is exact for the ID collapse of a slab. Grav- 
itational instability ensures that a uniform density el- 
lipsoidal density perturbation will collapse most rapidly 
along its shortest axis, forming a "Zel'dovich pancake." 

Baryons infalling along with the dark matter will pro- 
duce Lya forest absorption features. The line profiles will 
become distorted at the time the caustic form s before 
the baryons co me in to dynamical equilibrium (McGiljl, 
[1990; Meiksinl , Il994 Ikiralda-Escude an d Roes, 199^. 
The H I column densities produced will range over about 
13 < logiQ A^Hi < 15, wi th postshock te mperatures peak- 
ing at 50 - 60 X 10" K (jMeiksinl . Il994l ). 



D. Galactic models 

Ever since the earliest papers on the Lya forest, it was 
recognized that galaxies should give rise to intervening 
absorption systems. Their characteristics and frequency, 
however, are unclear. The interstellar medium of disk 
galaxies is expected to produce damped Lya forest ab- 
sorbers, although relatively few DLAs have been shown 
to arise from disks (see § IVIILA.l] below) . Other sugges- 
tions include co ndens ates in cooling outflows from dwarf 
galaxie s (,WanHl995l). the gas-rich dwarf galaxies them- 
selves ( York et al\ . 19861) . debris from satellite galaxies 
lost through tidal strippi ng by a ce ntral galaxy or ex- 
pelled by su pernovae (Morris and van den BergH Il994 
IWang| . ll993D . or pressure-confii ied clouds in a two-phase 
gaseo us galactic halo medium ( Mo and Miralda-Escudel 
I1996D . The outflow model in particular may give rise to 
distorted absorption line features if a line of sight inter- 
cepts the outflowing diffuse gas. 



Another suggestion is that the absorption lines arise 
from gas in extended disks with too low surface density 
to produ ce many stars, the "Vanishing Cheshire Cat" 
model of ' Salpeted In this model, an mner disk 

gravitationally-confined region of radius about 15 kpc is 
surrounded by an extended pressure-conflned disk of ra- 
dius about 250 kpc. The inner region has a low gas 
surface density compared with most disk galaxies, but 
sufficiently large to produce a DLA feature in any line 
of sight passing through it. The outer region has such 
a low surface density that the hydrogen is photoionized 
by the metagalactic UV background. The outer region 
will give rise to Lya forest features. A transition re- 
gion will produce Lyman Limit Systems. The relative 
abundances of the various H i column densities is ac- 
counted for by the progressively diminishing geometric 
cross section for increasing column densities. The transi- 
tion between gravitational confinement and pressure con- 
finement may be adjusted to account for deviations from 
a perfect power law in the column density distribution. 
The "Cheshire Cat" nature of the galaxy arises once star 
formation begins, producing a star burst galaxy. Metal- 
enriched gas ejected by the stars will accelerate the cool- 
ing and strengthen the burst, resulting in the expulsion 
of the gas in the central region and quenching further 
star formation. The result will be an outer "smiling" 
gaseous disk, possibly enriched in gas and metals by a 
galactic fountain if material from the starburst falls onto 
the disk, and a very low surface brightness "invisible" 
galaxy in the center. 

Aspects of the model, and the balance of gravita- 
tional and pressure forces in slab models generally , have 
been e xplo red bv [Ch arlton et al. (1993). Charlton et all 
(|l994f ) and lSalpeteTand HoffmanI (jl995[ ). 



VI. NUMERICAL SIMULATIONS 
A. Cosmological structure formation 

1. Dynamical evolution of the dark matter and baryons 

The successful description of the large scale distribu- 
tion of galaxies by the Cold Dark Matter do minated sce- 
nario of cosmological structure formation ( Davis et all 
Il985l ) led to attempts to account for the Lya forest within 
the same scenario . Using a comb ined hydrodynamics 
and gravity code, ICen et all ( 1994f ) demonstrated that 
the baryons will trace similar small scale structure as 
does the dark matter due to gravitational instability in 
the Cold Dark Matter scenario. Several other groups 
soon followed w ith sim i lar results ('Hernauist et a/.','l99fi ; 



Mir alda-Escu de et al... 
Zhang et aZ.I . ll995lll997D" 



9961 : [mdsley and Bond, 1997 



In the gravitational instability scenario of structure 
formation, cosmological structures grow from smaller 
structures through gravity. The growth of structures is 
governed by the coupled gravity-fiuid equations describ- 
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ing both the colhsionless gravitating matter component 
and the coUisional baryonic component. Because of the 
enormous change in spatial scales as the Universe ex- 
pands, it is most convenient to solv e the eq uations in a 
frame comoving with the expansion ljPeeblesl , |l99a) . The 



physical position r and velocity u are then expressed in 
terms of the comoving position x and peculiar velocity v 
according to 

r — a{t)x, u = r — ax + V, (88) 

where v = ax and a{t) = 1/(1 + z) is the cosmological 
expansion factor corresponding to the cosmological epoch 
at the redshift z. 

The dynamical equations for the dark matter compo- 
nent are 



dt 



Vd(x,t), 



dt 



d 1 , / N 

-Vrf = --V0(x, i), 
a a 



(89) 
(90) 



where 0(x, t) is the gravitational potential, and the gra- 
dient is with respect to the comoving coordinate system. 
The fluid equations expressing mass and momentum con- 
servation are 

MEili)+3^p+iv.[p(x,i)v(x,t)] = 0,(91) 
Ot a a 

(9v(x, t) a , , 1- , , „, , , 

^ + - v(x, t) + - [v(x, t) • v]v(x, t) = 

at a a 

-— Vp(x,<)-iv(/.(x,i), (92) 
pa a 

where /9(x, t) and p(x, t) are the gas density and pressure 
at x at cosmological time t. Because the evolution equa- 
tion for the entropy, Eq. (|5ip , is already in the comoving 
frame, it remains the same with the understanding that 
d/dt is replaced by d/dt+{l/a)v{x, t) ■ V (Peebles, ,1993) . 
The heating and cooling terms include the physical pro- 
cesses discussed above in § IIII.BI The evolution of the 
ionization state of a mixture of hydrogen and helium is 
given by Eqs. (|46p . The inclusion of ionization assumes 
a given set of Hi, He i and He ii photoionization rates 
and photoheating rates. These must be computed from 
the ambient radiation field. The initial ionization state of 
the gas is set during the passage of ionization fronts pro- 
duced by the radiation sources that photoionized the gas 
during the epochs of hydrogen and helium reionization. 
This is discussed in § IVIII below. 

The dark matter and fluid are coupled through Pois- 
son's equation for the gravitational potential 



V^^ = 47rGa2[ptot(x,t) - (ptot)], 



(93) 



where ptot = + P is the total mass density of the dark 
matter and fluid combined, (ptot) is the spatial average of 
Ptot over the Universe, and p <C (ptot — (Ptot))c^ has been 
assumed. In terms of the present day matter closure pa- 
rameter VLra, (ptot) = 3-ff^ri,„/(87rG'a^). Eqs. 



are valid for any Friedmann-Robertson- Walker universe, 
including those with non-zero vacuum energy (such as an 
effective cosmological constant). 

These are the same equations that describe the mix 
of dark matter, baryons and radiation that give rise to 
the Cosmic Microwave Background, except for the addi- 
tion of the radiation energy density and pressure. The 
principal difference is that while the density fluctuations 
are linear for the CMB up to the time the baryons be- 
gin to recombine and decouple from the radiation, the 
description of the growth of cosmological structures ex- 
tends the computation of the fluctuations into the non- 
linear regime. Numerical solutions to Eqs. ((89|) - ([93]) 
generally begin at a sufficiently high redshift that the pri- 
mordial density fluctuations are still linear on the scales 
of interest, but after the matter and radiation have de- 
coupled. The most convenient set of boundary conditions 
for solving the PDEs for the coupled dark matter and gas 
on cosmological scales is to assume all the variables are 
periodic over the simulation volume. It is therefore im- 
portant that the simulation volume be sufflciently large 
to represent a fair sample of the Universe. 

No structures will form unless the matter is initially 
inhomogeneous. The inhomogeneity is usually specified 
in terms of a primordial power spectrum of matter fluc- 
tuations, which is the ensemble average of the squared 
modulus of the Fourier modes of the density field. The 
mass density expanded in Fourier modes within a peri- 
odic Cartesian volume Vbox may be expressed as 



p(x,t) = (p) 



1 + ^ i5k(i) exp(— ik • x) 

fc#0 



(94) 



where the Fourier amplitudes (5k (i) are assumed to have 
random phases and to be distributed according to a gaus- 
sian random process. The power spectrum, which quan- 
tifies the amplitude of the density fluctuations, may be 
defined as P{k) — Vbox(|<5kP), where the average is over 
independent statistical realizations of the Fourier ampli- 
tudes. 

The real space analogue of the power spectrum is the 
spatial correlation function ^(r) of the density field 



(p(x' + x)p(x')) = (p>^[i + e(x)]. 



(95) 



where the average is carried out over all spatial locations 
x'. The correlation function is related to the power spec- 
trum through a Fourier transform 



e(x) 



1 



(27r)3 



d^k 



P{k)- 



1 



exp(— ik-x), (96) 



where rip is the mean spatial density of the particles, and 
the Poisson noise contribution due to an average num- 
ber of Tip Vbox particles per simulation volume has been 
subtracted. 

It is convenient to define the amplitude of density fluc- 
tuations 



A^(fc) 



_fc3 

27r' 



(97) 
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The amplitude of mass fluctuations {SM/M)'^ in a sphere 
of radius r f is then given by 



kr 



7 



(98) 



where ji{x) — (sinx)/x^ — (cosa;)/x is a spherical Bessel 
function. A common fiducial normalization measure 
is the mass fluctuation filtered on a scale of = 
8/i"i Mpc, where h = i?o/100kms"^ Mpc"^ The spec- 
trum of matter fluctuations may be expressed in the form 



3+n 



D2(0) 



T2(fc), 



(99) 



where T{k) is the transfer function describing the evolu- 
tion of density perturbations as they cross the present- 
day horizon scale with amplitude 5h, n is the spectral 
index, and D{z) is the growth factor for the density per- 
turbations (see § IV. CI) . The transfer function for a Cold 
Dark Matter d ominated universe was first computed by 
iPeeblesI (|l98lh . A fast numerical code for computing 
the transfer funct i on to high accuracy is provided by 
lEisenstein and Hul (|l999l ). 



2. Methods of numerical simulations 

Several numerical methods have been developed for 
solving the combined gravity-fluid equations. The most 
common methods for solving the dark matter evolu- 
tion equations sample the distribution function of the 
dark matter using massive pseudo-particles and solve for 
their evolution either on a mesh or using a hierarchi- 
cal tree. The most common mesh codes use either the 
Particle-Mesh (PM) or Particle-Particle Particle-Mesh 
(P^M) techniques. Tree algorithms instead arrange the 
particles into a hierarchy of groups and approximate the 
forces exerted by distant groups on a given particle by 
their lowest multiple moments. The fluid equations are 
frequently solved using a finite-difference representation 
of the fluid equations, handling shocks either with arti- 
ficial viscosity or through shock capturing. An alterna- 
tive scheme, Smoothed Particle Hydrodynamics (SPH), 
smooths the fluid equations and samples the smoothed 
fluid variables using pseudo-particles. The algorithm 
is intrinsically Lagrangian in nature. A comparison of 
two popular publicly available codes using thes e meth- 
ods^NZO and GADGET, is provided bv O'Shea et al\ 
(|2005l ) . The paper also serves as a good entry point to the 
vast literature on the subject. A more general overvie w 



of A''-body methods is provided bv iBertschinged ( 1998f ). 

All of the techniques share a common method for 
initiating the density fluctuations and associated pecu- 
liar velocity field bas e d on the Zel'dovich approxima- 
tion ("Efstathi ou et aLl , Il985l ). Substituting = akd — 
—aD{t)W'^{qd), using Eq. ([57]) for the dark matter, into 



the equation for peculiar velocity from Eq. (190]) gives 
1 



*(qd) 



(q^, t) + constant, (100) 



after noting that aD + 2aD = 47rG(ptot)a-D, which fol- 
lows from Eq. I|9ip applied to the combined dark mat- 
ter and gas density. The initial deformation is therefore 
proportional to the gravitational potential (j) (up to an 
arbitrary integration constant). This offers a straight- 
forward prescription for specifying the initial conditions. 
First, choose t he (5k ' s according to a gaussian random 
process: Si^ — \J P{k)g exp{i9), where g is a gaussian ran- 
dom deviate between — oo < g < oo with zero mean and 
unit variance, and is a random phase uniformly dis- 
tributed over < 6 < 2tt. In Fourier space, the gravi- 
tational potential is (pu = — 47rGa^(ptot)<5kfc~^ since the 
Fourier transform of V^cj) is —k'^4>k, where 4>k{t) is the 
Fourier transform of 0(q, t). The gravitational potential 
may then be constructed from 



(q,^) = ^0kexp(-ik • q) 

- -AnGa^iptot) J] ^ cxp(-zk • q). (101) 



The deformation function ^'(qd) then follows from 
Eq. poop . The displacements and velocities of the parti- 
cles then follow from Eqs. ([57)l . 



X = q-^^(5kexp(-ik-q), 
V = -a:^^!^4exp(-zk-q). 



(102) 
(103) 



The optimal form of these equations allows for convolu- 
tion with the differencing operator use d to evaluate the 
gradi ent of the gravitational potential ( Efstathiou et all . 

Once initiated, any of the various numerical meth- 
ods described above may be used to evolve forward the 
dark matter and baryon fluctuations into the nonlinear 
regime. 

The resolution requirements of the simulations de- 
pend on the statistical quantity for which convergence 
is sought. While the demands for convergence on the 
H I column density distribution over the range 13 < 
A^Hi < 15 are not very severe, the early simulations 
tended to over-estimate the width of the absorption fea- 
tures. It is now recognized that good convergence on 
the Doppler parameter distribution and of the pixel flux 
distribution requires at the minirnum good resolution 
of th e Jeans length ( Brvan et all Il999l : [Theuns et all . 
Il998| ). A resolution of 9h~^ kpc (comoving) is about 
the right level at z = 2, but still inadequate by z > 3 
( Brvan et all . \l99dl} . although the rate of convergence is 
sensitive to the gas temperature and assumed cosmol- 
ogy. Slightly coarser resolution is adequate for converg- 
ing on the pixel flux distribution. A large simulation 
volume is required to ensure that the full contribu t ion of 
long wavelength modes is included (jBrvan et all . Il999t 
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iTheuns et al. I, 119981: IWadslev and Bondl . [imJ), requir- 
ing a box size of at least 5h~^ Mpc (comoving). Sim- 
ilar demands are found using pure dark matter pseudo- 
hydrodynamics simulations, for which sequences of large 
high resoluti on simulations have been p erformed to test 
convergence (jMeiksin and Whitel . l200lf l. Matching both 
the resolution and box size requirements demands simula- 
tions of minimally 512'^ zones, and preferably 1024^ zones 
for grid based codes to achieve good convergence on the 
properties of low H i column density systems. Similar 
numbers of particles are required for SPH computations. 
Even with these parameters, converging on the extremes 
in the temperature distribution proves difficult. Such 
large hydrodynamics simulations are only now becoming 
feasible. 

The simulations show that the baryons closely trace 
the distribution of dark matter down to the Jeans length 
scale, when pressure forces prevent further collapse of 
the baryons. This has prompted an alternative approach 
to IGM simulati ons based solely on the evolution of 
the dark matter (ICroft et aLl . ll998HMeiksin and Whit3 . 
I2OOII : iPetitiean et all 1995 ). The gas overdensity is ei- 
ther set equal to the dark matter overdensity or by al- 
lowing for a pressure by treating it as a modification to 
the gravitational potential in the Hydrodynamics Parti- 
cle Mesh (HPM) approach (|Gnedin and Huil . ll998D . The 
gas temperature is either solved for all owing for photo- 
electr ic heating and radiative cooling ( Petit lean et all 
I1995D or through a polytropic equation of state T = 
Tj^ ipl (p)V ■ This al t ernati ve is based on the recognition 
bv lHui and GnedinI ( 1997t) that the temperatures found 
in numerical hydrodynamics simulations are determined 
largely through a balance between photoelectric heating 
and adiabatic cooling losses for structures with overden- 
sities < 5. Comparisons with full hydrodynamical 

computations show that using a PM code or HPM code 
with a polytropic equation of state gives a reasonably 
good description of the prop erties of the IGM, i nclud- 
ing its absorption sta tistics ( Meiksin and Whi"t3 . I2OOII : 
IViel et aLl . [2006ll2002D . although full hydrodynamics sim- 
ulations should still be used for detailed comparisons to 
the data. 

Yet another approach altogether is to generate the den- 
sity fluct uations stochastically from a lognormal d i stribu - 
tion (jH, Il993t iBi and Davidsenl . Il997t iBi et all Il992t) . 
This approach captures many of the qualitative features 
of the structure of the IGM, and provided an early argu- 
ment that the Lya forest itself may be accounted for by 
density fluctuations in a diffuse IGM. 



B. Physical properties of the intergalactic medium 

The results of structure formation simulations show 
that the IGM fragments into an intercon nected net 
work of sheets, filam e nts and halos (ICen et al. 



1994; Katz et all 



Miralda-Escude et all 11996 : 



ious morphologies correspond to different physical den- 
sities and different H i column density absorption sys- 
tems. These dependences are illust rated in Figures [T^ 
and [131 taken from the simulation of'Zha ng et all (|l998l ) 
of a Cold Dark Matter dominated E instein-deSitter uni- 
verse with the radiation field from iHaardt and Madaul 
(I1996D based on QSO sources. The simulation volume 
was 9.6 Mpc (comoving) on a side, and the power spec- 
trum was normalized to erg = 0.7 with h = 0.5 used. 
At z = 3, systems with iVni ^ 10"'^^ cm~^ arise from the 
amorphous underdense regions enclosed by the sheets. 
The underdense regions have characteristic diameters of 
a few comoving Mpc, and occupy most of the volume 
of the Universe. The sheets have densities compara- 
ble to the mean baryon density, and give rise to ab- 
sorption systems with 10^^ cm~^ ^ A^hi ^ 10^** cm~^. Sys- 



tems with 10 cm 



HI , 



10 cm originate in the 



filaments, which arc moderately overdense structures 
(l^;p/(p)^;5). Higher column density systems, with 
Nm ^ 10^^ cm^^, arise in spheroidal nodes with overden- 
sities of several to over a hundred embedded within the 
filaments. 

The mean temperatu res of overdense systems range 
over lO'* K to 3 X 10^ K (ICen et all [1991 'Theuns et 



1998t IZhang et all Il998( ), while [Miralda-Es cude al] 



( 19961 ) find mean temperatures up to 6 x 10^ K. These 
larger values may be due to the hard UV background 
assumed. Much higher temperatures are also found, 
extending up to 10^ — 10'' K, for ga s associated with 
collapsed halos ( Theuns et all Il998l ). There is less 
agreement on the temperatures in the underdense re- 
gions. The time to achieve thermal e quilibrium ex- 
ceeds a Hubble time for underd ense gas ( Meiksinl . 1 1994 
iMiralda-Escude and Reed . Il994[ ). The gas temperature 
will therefore reflect the thermal history of the gas fol- 
lowing reionization, and so is sensitive to the particular 
reionization model. The Lya forest simulations also did 
not t ake into account the heatin g effects of reionization 
itself ([Tittlev and Meiksinl . [200I . 

The thicknesses of the filaments are approximately 
100 kpc (proper), with coherence lengths of about a 
Mpc. The scales are consistent with the observational 
constraints based on neighboring line-of-sight statistics 
discussed in § III. B. 31 above. This is comparable to the 
isothermal Jeans length. For a characteristic hydrogen 
number density of 55 m~'^, corresponding to 5 times the 
mean baryon density at z = 3, and a characteristic tem- 
perature of 2 X 10^ K, the pressure is 10^ Km^'^ and the 
Jeans length from Eq. ([531) is 320 kpc. 

The baryonic overdensity at z = 3 is found to 

scale roughly as p/{p) ~ 3{Nm/lO^'^ ciR-^y^^ , 

with the neutral fraction scaling as x m — 1.3 x 
10-5(ArHi/10"cm-2)i/27^-o.75p - 



HI. -12 



Wadslev and Bondl Il997[ : IZhang et all\l99^ . The var 



([Zhang et all 

Il998f ). The neutral fraction dependency is comparable 
to the estimate in § III. B. 31 above, and results in virtually 
all the baryons' residing in the Lya forest at 2: > 1. 

The structure of the IGM at z <C 1 dif f ers so mewhat 
from its structure at z > 2. [Dave et al. I ([1999D obtain 
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FIG. 13 Isodensity contour surfaces of log baryon overdensity at 2 = 3 for a flat Cold Dark Matter dominated universe. The 
contour levels are logjQ(p/(p)) = —1.0, —0.5, —0.3,0.0,0.5 and 1.0. Low density regions are amorphous structures filling sheets 
at the mean density which intersect at overdense filaments. The filaments intersect at highly overdense collapsed halos. From 
IZhang et al\ l|l998f ). (Figure reproduced by permission of the AAS.) 




FIG. 14 Evolution of the H i column density. The filamentary structure of the IGM is nearly stationary in the comoving frame. 
By contrast, the decline in the physical g as density results in rapid evolution of the H i column density, and the IGM becomes 
increasingly transparent with time. From lZhang et al\ (|l998l ). (Figure reproduced by permission of the AAS.) 
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the relation between overdensity and H i column den- 
sity of p/{p) - 20(iVHT/10 ^ ^cm- ^)' " 



^ -0.42^ similar to 

the finding of IZhang et al\ ( 1998| l at z = 3. A given 
overdensity results in H i column densities smaller by a 
factor of 10-50 at z ^ 1 than at z > 2. Most absorption 
systems at low redshift arise in systems with hydrogen 
densities of tih — 0.1 — 10 m^-^ for H i column densities 
of -/Vh — 10^^ — 10^^ cm^'^. Structures retain coherence 
over scales of several hundred kiloparsecs. The gas tem- 
peratures for structures with p/{p) < 10 are about a 
factor of 2 smaller at z ^ 1 than at z > 2. A hot, 
intermediate-density gas phase also arises from shocks in 
filaments and in galaxy, galaxy group and galaxy cluster 
halos, reaching temperatures of over 10^ K. 



C. Statistical properties of the absorption systems 

1. Mean metagalactic ionization rates 

Almost all that is directly known about the structure of 
the IGM is derived from spectra. The simulation results 
must therefore be translated into spectra for comparison 
against observational data to test the validity of the mod- 
els. The spectra, normalized to a unit continuum level, 
may be expressed as exp(— Tj^), where t^, the IGM opti- 
cal depth as a function of frequency, is given by Eq. (fT4|) . 
Allowing for a peculiar velocity v varying along the line 
of sight. 



(104) 



Setting n = {n){l + S), where 6 is the local density con- 
trast, Eq. p04p is identical to Eq. (fT5)) for the Gunn- 
Peterson effect in the absence of perturbations {v — 
6 = 0). For this reason, the varying optical depth 
of the IGM is sometimes referred to as giving r ise to 
a "fluctuating Gunn-Peterson e ffect" ( Gen et al\ . 1 19941 : 
iMiralda-Escude and Re"e3 . ll993l) . 

An exampl e of a spec t rum g enerated from the AGDM 
simulation o f'Gen ei all (|l994l ) with = 0.6, = 0.4, 
h = 0.65, and normalization erg = 0-79 is shown in Fig- 
ure [151 The temperature variations closely trace the 
baryon density fluctuations. The absorption features are 
slightly offset from the density features due to the pecu- 
liar motion of the gas. Mild overdensities produce weak 
features while large overdensities produce stronger fea- 
tures. The apparently isolated absorption features arise 
not from distinct clouds but from modulations in the 
density field produced by the gravitational instability of 
primordial density fluctuations. 

Any comparison between the simulated and measured 
spectra requires normalization by the uncertain ioniza- 
tion rates for H i and He ii. The uncertainty in the 
normalization rate for He ii is compounded by the un- 
known epoch of He ii reionization, which may have been 
as late as z ~ 3. Late He ii reionization will also af- 
fect the gas temperatures, particularly for underdense 




FIG. 15 A line-of-sight cut through a ACDM simulation of 
the LyQ forest ai z = 3. The top panel shows the gas tem- 
perature (solid line), gas density (thick dotted line) and gas 
pressure, in arbitrary units (thin dotted line). All are shown 
as a function of the spatial coordinate r expressed in velocity 
units through ii = Hr, for H = 512/ikms~^ Mpc~^. The gas 
temperature and density are closely correlated. The middle 
panel shows the peculiar velocity (dotted line), the gravita- 
tional acceleration divided by the Hubble constant (dashed 
line) , and the total hydrostatic and gravitational acceleration 
divided by the Hubble constant (solid line). An infalling pe- 
culiar velocity pattern forms around a collapsing structure 
near v ~ 360 km s~^. The bottom panel shows the transmis- 
sion factor exp(— Ti^) before convolving with a Doppler profile 
(solid line), after convolution with a Doppler profile (dashed 
line), and after adding the eflf ects of inst r umen tal resolution 
and noise (dotted line). From lCen et al\ (ligO^ I. (Fi gure re- 
produced by permission of the AAS.) 



gas. This in turn affects the optical depth through the 
gas. Much of the uncertainty on global flux properties 
may be removed by normalizing to the measured mean 
Lya transmission (exp(— Tj^)) of the IGM. Although an 
error in the temperature of the gas will result in a small 
offset in the inferred ionization rate (see Eq. [42]), once 
normalized, the probability distribution of the emissivity 
is only slightly affected. Altering the temperature of the 
gas as required by the line-widths, either directly or by 
boosting the heating rate, can further improve the flts. 

Various groups have provided estimates for the H i- 
ionization rate Fhi based on matching their simulation 
results to measurements of the mean Lya transmissivity. 
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Result s fr om iMeiksin and Whit e (2004), iBolton et dj] 
(|2005f ) and lBolton and Haehnelti (,2007. ) are shown in Ta- 
ble |Vl All the predicted rates a r e in reasonably good 
agreement. ' Bolton and Hachneltl (|2007l ) provide an ad- 
ditional estimate at z = 2 of Phi 



1.291^" which 



agrees well with the estimate from iTvtler et ali ( 2004a[ ) 
of rHi,-i2 — 1-44 ± 0.11. An additional uncertainty of a 
few to several percent may result from fluctuations in the 
ioniz i ng background jC roft, 20041: iGnedin and Hamiltonl . 
120021: IMeiksin and Whita . ,20041) . A comparison of the 
required ionization rates with the predictions from QSO 
sources shows that unless there is a large population of 
low luminosity AGN so far undetected at high redshifts, 
QSOs alone are unable to provide the required ionization 
rates at z ^ 3. 

Specifying Fhi fixes the entire equilibrium ionization 
structure of the gas (except in isolated optically thick 
regions). This permits the contribution of the IGM to 
the attenuation length of ionizing photons to be deter- 
mined directly from a simulation. In principle, a self- 
consistent solution for the equilibrium structure of the 
IGM may then be obtained, including the magnitude of 
the requir ed emissivity for any given s pectral shape using 
Eq. IMeiksin and Whitel (|2004 ) find for a ACDM 

cosmology the (proper) attenuation length is approxi- 
mated to 10% accuracy by 



1.7 X 10*(l + z)-''-2/i-iMpc 



(105) 



for 2.75 < z < 5.5. Because the simulation underpredicts 
the number of Lyman Limit Syst ems by a large factor, 
these are added using the results of lStengler-Larrea et all 
(I1995D . - t is found that the diffuse gas comprising the 
Lya forest and the Lyman Limit Systems contribute 
about equally to the attenuation le ngth, consistent with 
mode ls based on direct line counts (jMeiksin and Madaul . 
I1993I ). The implied (comoving) metagalactic emissivity 
at the L yman edge requ ired to recover the H i ionization 
rates is (jMeiksinl . I2005D 



(106) 



where Amg = 8.4l4;5;g'5 x lO^^ and 7 = -1.6 - 
0.61n(^MG/8.4 X 10^^) and umg is the spectral index 
of the metagalactic emissivity, e^, — eL{v / vl)~°'"^ . The 
errors on Amg are the la and 2a limits. 

At z > 4, Eq. (|106p is well above the predictions from 
QSO counts. As discussed in g IIV.A. 21 above. Eq. (fTTOl) 
may be matched by the emissivity of galaxies at the Ly- 
man edge for plausible escape fractions of the ionizing 
radiation. 



2. Spectral properties of the IGM 

The predicted spectral properties of the Lya forest 
are sensitive to the cosmological parameters used 
(|Jena et ali l2005t iMachacek ^tdl . I2OOOI : iRauch et all . 




FIG. 16 Comparison of the cumulative flux per pixel distri- 
bution from a ACDM simulation with data based on high res- 
olution Keck spectra of 7 QSOs. The data have been grouped 
into three redshift intervals and the observations corrected to 
take into account evolution in the mean transmission within 
each interval. Noise has been added to the simulated data 
to match the noise p roperties of the measured spectra. From 
IRauch et all l|l997bl ). (Fi gure reproduced by permission of 
the A AS.) 



Il997bl : iTheuns et g/l Il999l: IWadslev and Bondl . [l993). 

Comparisons with the spectral properties of the 
Lya forest have shown impressive agreement f or plausible 
cosmological models. iMachacek et al. 1 (I2OOOI) show that 
the distribution of optical depth per pixel in the models 
is well-fit by a lognormal distribution for values between 
0.02 < T < 4, roughly the measurable ran ge. This pre- 
dictio n is confirmed by the observations of lBecker et all 
(I2OO7D . who find the flux distributions in 55 high reso- 
lution QSOs are well-reproduced by lognormal distribu- 
tions for the underlying optical depths. Detailed compar- 
isons show that the predicted cumulative distributions of 
pixel fluxes agree to a precision of a few p ercent with 
meas urements in individual QSO spectra ((Meiksin et all . 
I2OOII) . 

These results suggest the Lya forest may be used as 
an effective tool for selecting between rival cosmologi- 
cal models of structure formation. An example of a 
comparison between the ACDM simulation above and 
high signal-to- noise ratio high res olution Keck spectr a 
of 7 QSOs is sh o wn in Figur e [TCI (iRauch et all. Il997bl) . 
IMachacek eTall (|2000l) and IMeiksin et all (|200ll ) find 
that the shape of the pixel flux distribution and the 
magnitude of the Doppler parameters vary systematically 
with the normalization on small scales, such as aj filtered 
on the Jeans scale, as defined by Eq. (|98)) with filter scale 
rf — Xj from Eq (|84p . A detailed comparison between 
the predicted flux distributions for several cosmological 
models and measurements from several high resolution 
Keck spectra constrains aj, defined for an isothermal gas 
at T = 2 X 10 ' ^ K, t o the range 1.3 <crj< 1.7 at z = 3 
( Meiksin et al ]. l200lh . although this is somewhat subject 
to the uncertain thermal properties of the gas. The clos- 
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est agreement is found for the larger values of aj. A range 
of ACDM models, however, shows that agreement with 
as inferred from WMAP n ieasurements of CMB fluc tua- 
tions prefers lower values ("Mei ksin and Whitd . l2004 f) , so 
that there is some tension between the higher normaliza- 
tion preferred by the Lya forest simulations and the lower 
normalization preferred by the CMB measurements. 

Voigt profile fitting to the spectra permits compar- 
ison between the predicted and measured H i column 
density and Doppler parameter distributions. Here the 
models fair somewhat less well compared with the pixel 
flux distribution. The models successfully reproduce the 
H I column density distributions to good accuracy over 
a wide range of values of 14 < log, n Nm 5 16 at high red 
shifts (2 < z < 4) (fCen et all Il99l iHernauist et al 



19961: 



.Meiksin et all joOlV IWadslev and Bondl . Il997 : 

Zhang et all Il995l 11991" and at moderate to low red- 
shifts (0 < z < 2) (|Dave efd\ . Il999t iPaschos et d\ . 



1^008.) . A discrepancy is found at higher column den- 
sities, with simulations generally finding fewer systems 
with logj^Q A'hi > 16 than measured. This is particu- 



larly true for Lyman Lim it Systems ( Katz et all Il996t 
iMeiksin and Whitd . \2004) . although extrapolating the 
simulation results to smaller mass halos suggests agree- 
ment is possible, depending on the appropriate minimum 
circular velocity of halos able to produce high H i column 
density systems in the i netag alactic photoionizing back- 
ground ( Gardner et aZ.l . l200lh . A small er discrepancy is 
found for DLAs. The simulations of iNagamine et all 
(|2007l ) suggest the incidence of DLAs may be met if 
they are associated with halos having masses as low as 
Mq. 

It is less clear that all Lyman Limit Systems arise from 
collapsed halos. Allowing for radiative transfer effects. 



iKohler and GnedirJ ( 2007 ) find their simulations produce 
Lyman Limit Systems from material bound to and falling 
onto galaxies rather than from minihalos. It is possible 
the systems have more than a single origin. 

Because the detailed distribution of the absorption line 
parameters is sensitive to the method of line finding and 
fitting (see § III.B.lj) . a definitive comparison between 
the simulations and the measured absorption line prop- 
erties requires identical line analysis methods. Automat- 
ing the line analysis is necessary for processing the large 
number of lines of sight from simulations to build up 
reliable statistics. A few codes are publicly available. 
An industry standard is VPFIT,^^ a sophisticated code 
that handles in addition to H i features a wide range of 
metal lines. A simpler alternative well-ad apted to au- 
toma ted fitting of H I features is Au toVP ( Dave et all 
Il997t ). Based on fits with AutoVP, iDave et all (|l997n 
find that the simulations predict a greater fraction of 
lower Doppler parameter syste ms than measured. Simi- 
lar discrepancies were found bv lBrvan et all i 19991 ) and 
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FIG. 17 Comparison of the cumulative H i column density 
(upper panels) and Doppler parameter (lower panels) distri- 
butions for the Lya forest measured along the line of sight 
to the QSO Q1937-1009 with predictions for six cosmological 
models: a standard CDM model (SCDM), two ACDM models 
with high (ACDMh) and low (ACDMl) baryon densities, an 
open CDM model (OCDM), a tilted CDM model (tCDM), 
and a mixed hot darlc matter and cold dark matter model 
(CHDM). Much better agreement is found for absorption line 
systems with line center optical depth s rp > 1 (right panels ) 
than for ro < 1 (left panels). From iMeiksin et all l|200ll 'l. 
(Figure reproduced by permission of the RAS and Blackwell 
Publishing.) 



http: / /www. ast.cam.ac.uk/ ~rfc/vpfit.html 



iTheuns et aZI (|l999l ). 

A detailed comparison between several high resolu- 
tion QSO spectra and numerical simulations for a variety 
of models, analyse d using AutoVP and a similar code, 
SPECFIT, bv Meiksin et all (j200ll ) shows that it is not 
possible to match simultaneously the pixel flux distribu- 
tions, column density distributions, and Doppler param- 
eter distributions using simulations with optically thin 
reionization. While adjusting ctj may help to match the 
Doppler parameters, it destroys agreement with the pixel 
flux distribution. In Figure I17[ a comparison is shown 
between the predicted and measured line parameter dis- 
tributions for the Lya forest detected in the spectrum of 
the z = 3.8 QSO Q1937-1009. The systems have been 
isolated by their line center optical depths. The compar- 
isons show that the discrepancies are due primarily to 
the low optical depth systems, with line center optical 
depths To < 1. 

The origin of the discrepancies is most likely the tem- 
perature of the gas. Artiflcially broadening the absorp- 
tion features gives excellent agreement between the mea- 
surements and predictions for the ACDMl model for the 
pixel flux and column density distributions, and good 
agreement between the Doppler parameter distributions. 
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Since the optically thin absorbers arise in underdense re- 
gions, the temperature of the gas will retain a memory of 
its heating history. Since the underdense gas occupies 
most of the volume of the Universe, whatever heated 
the gas must have been ubiquitous. The best candi- 
date is late reionization of He ii, in the redshift interval 
3 < z < 4. This would boost the gas temperature, pos- 
sibly by the required amount. The effects of heating by 
late He ii reionization may be mi micked to some extent 
by boosting the He ii heati ng rate f|Brvan and Machaceld . 
I2OOOI: iTvtler et all l2004aD . but the only fully satisfac- 
tory solution requires inco rporating radiative transfer 
(jTittlev and Meiksinl . [20071 ) . 

With the goal of developing a physical description 
for th e origi n of the Doppler parameter distribution, 
iHui and Rutledgi (fl999.) introduce a semi-analytic ap- 
proach based on the growth of the density perturbations 
to derive a Doppler parameter distribution of the form 



■ exp 



64 



(107) 



where ha is a fit parameter. Restricting their analysis to 
systems with column densities 12.5 < logj^g TVhi < 14, 
corresponding to the physical condi tions of the IGM 
for wh ich their analysis best applies, iHui and Rutledg"^ 
( I999I ) find the distribution provides statistically ac- 
ceptable matches to the Doppler parameter distribu- 
tions of the Lya forest along the lines o f sight to 
four QSOs with high resolution Keck spectra ( Hu et ai\ . 
Il995l ). Adequate match e s, howe ver, were not found 
for QOOOO-26 (iLu et al\. | l996bD or HS 1946-1-7658 
( Kirkman and Tvtleii 1997a ). The Doppler parameter 



distributions measured in the spectra of all six of the 
QSOs, without any restriction on the column density 
range app lied, are well- f i t by lognormal distributions (see 
Table nil). iRicotti et all (|200Q) generalize Eq. pUT]) to 



exp 



(108) 



where a is a fit constant, and Lp, <& and 6* are allowed 
mild dependences on the column density. The best fits 
obtained for absorbers in several lines of sight over the 



redshift range 2.5 



4 have 1.6 < a < 2.1, 2.6 < (/? < 



5.9 and 2.0 < $ < 7.2, with values almost independent 
of column density. 

Several attempts have been made to infer the temper- 
ature of the IGM directly from the Doppler parameter 
distribution. One approach is to normalize the distri- 
butions using simulations, invoking additional sources 
of heat input to allow for a range in temperatures 
( Brvan and Machacekl [2003)- Another exploits a poly- 
tropic relation be t ween t he gas densit y and temperature 
(|Hui and Gnedinl . Il997l: ISchave et g/l Il999l ). calibrated 
using numerical simulations. Both methods rely on nu- 
merical simulations that neglect the effect of reionization 
on the tempera ture of the underdense g a s, which may 
be substantial ( Abel and HaehnelO . Il999l : iBolton et al\ . 



12004 iTittlev and Meik"siiil . l2007l ) . Modeling the Doppler 
parameter distributions according to Eqs. (11071) and 
(|108 ) to obtain lo \ ver cut off values, IRicotti et al 



Schave et al\ d^OOOb) find evidence for a tempera- 



and 

ture jump and reduction in the polytropic index towards 
isothermality near z = 3.0. Both groups attribute their 
results to late He 11 reionization. He 11 reionization sim- 
ulations predict a more gradual change (McQuinn et ai\, 
I2OO8D . No significant change in the temperature or poly- 
tro pic index over the redsh i ft ran ge 2.4 < z < 3.9 is found 



tro pic maex over tne reasnnt ran ge z.4 < z < 6.\) is tounc 
bv iMcDonald et al\ (|200ll . I2003D . Although any conclu 



sion using these methods is subject to the accuracy of the 
modeling assumptions made, the possibility of measuring 
an unambiguous temperature jump in this redshift range 
merits further investigation. 



3. Intergalactic helium absorption 



As for the H i absorption features, the numerical simu- 
lations predict the properties of He 11 absorption as well. 
The systems that dominate the measured He 11 optical 
depth s arise from t he diffuse structures in u nderdense re- 
gions (jCroft et aZ.l . ll997HZhang" et al. I. fl995i) . This makes 
predictions for the He 11 optical depths difficult for two 
reasons: 1. high spatial resolution is required to accu- 
rately compute the structures in the underdense regions, 
and 2. the gas temperature is sensitive to the history of 
reionization since the time for the temperature to reach 
equilibrium is comparable to or longer than a Hubble 
time. The latter affects both the ionization fraction of 
He II because of the temperature dependence of the ra- 
diative recombination rate and the Doppler widths of the 
absorption features. Another uncertainty is the epoch 
of full heli um reionizat i on. T he results of IZhang et al\ 
(|1995D and IZhang. et al\ ([199H) suggest matching to both 
the H I and He 11 Lya mean optical depth measure- 
ments at 2.5 3.3 requires an ionization rate ratio of 

~ 250 — 400, comparable to the inferred values from 
line statistics f§ IIV.C.2|) and corresponding to an intrin- 
sic QSO spectral i ndex of 1.8-2.0, all owing for filtering 
through the IGM ( Haardt and Madaul . 119960 . A similar 
estimate is obtained from high r esolution numerical sim- 
ulations bv iBolton et al\ (|2006( ). who find values for ^ 



-1-576 

-151 



at z 



as- 



rising from 139167 at z = 2.1 to 301 
suming a spatially uniform UV background. Most of the 
uncertainty stems from the uncertainties in the measure- 
ments of the He 11 Lya mean optical depths. Allowing 
for spatial fiuctuations in the QSO contribution lowers 
the inferred values for ^ by about a factor of two. Al- 
lowing for fluctuations also reproduces much of the ob- 
served scatter in the He 11 to H i column density ratio 77 
{'jj IIV.C.2)) . except for the extreme excursions at < 10 
and ?7 > 200. 
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4. Metal absorption systems 

By providing the physical scales of density structures 
and their associated temperatures, the simulations per- 
mit estimates of the ionization states of metals for any 
assumed UV radiation field. Several observed features 
are reproduced by the models. The metal absorption 
features arise within the substructure of systems that 
give rise to the H i Lya features, and may result in 
multiple features corresponding to a sing l e H i systera 
(iHaehnelt et a/.l.ll996l:lHellsten aH . llQQTtlRaiich al\ . 
ll997al : IZhang et a^.l . ll997^ . The frequencies of C ii, C iv, 
Si IV and N v systems are comparable to those observed 
at z ~ 3 assuming a QSO/AGN-domin ated UV back- 
groun d and metalhcities of [C/H] ~ —2.5 ( Hellsten et all 
I1997D . The C IV to H I column density ratio decreases 
with increasing H i column density, since higher column 
densities correspond to systems with higher densities and 
so lower ionization parameters ( Hellsten et al\ . ll99T\) . 

iHellsten et all (|l998bl ) rank a variety of metal ions 
by their degree of detectability. The strongest absorp- 
tion lines are the O viAA1032,1038A doublet for low col- 
umn density systems, with A^hi < 10^^ cm^^. These 
lines unfortunately are very difficult to extricate from 
the Lya forest. The measured levels of O vi associated 
with C IV systems at z ~ 3 requires an overabundance 
of [0/C] ~ 0.5, assuming a QSO/AGN-dominated UV 
backgr ound, suggestive o f enrichment by Type II super- 
novae (iDave et al If helium was not fully ionized 
by z ~ 3, the resulting break between the H i and He ii 
Lyman edges requires a much higher overabundance of 
oxygen. At least half the IGM must have helium fully 
reionized to keep [0/C] < 10. 

Longward of the H i Lya line, the strongest lines 
are the C ivAA1548, 155lA doublet for all systems with 
A^Hi < lO^'^cm"^. In high column density systems, the 
regions giving rise to C ii absorption sometimes do not 
coincide with those giving the strongest C iv absorption, 
resulting in velocity offsets between the two. 

At low redshift, the formation of galaxy groups and 
clusters will produce an outer region of high density 
p/{p) ~ 100 and temperatures of T ~ 10^ - 10'' K, which 
will give rise to strong absorption lines of O vi (for tem- 
peratures 5.3 < logior<6.3), O VII (5.3 < logioT<6.7) 
and O VIII (6 ^ logj^g T ^ 7.3) in a warm-hot intergalactic 



let and x-rav (Cen and Ostrikeil 19991: Chen et al. 


20031: iDave et all 


20011: iHcUsten et al. Il998a 


Perna and Loebl. 19981 


Viel et al... ,2003a: Yoshida et al. 



20051 ). Especially of interest is the O viAA1032,1038A 
doublet, which has already been detected in numerous 
syste r ns at low redshifts (jPanforth a nd ShuU, 2005, 
l2008t iTripp et all . l2008l) . For relatively high metalhci- 
ties (at least 0.1 of solar), a few higher ionization species 
features may be detected by Chandra and XMM-Newton 
(|Chen et oi.l . l2003D . These along with lines of carbon, ni- 
trogen and possibly neon and iron may be detectable by 



the proposed x-ray satellites XEU^° iChen et a/.l. 12005 : 
Viel et al\ . l2003a( l and Constellation- Ji?^ (jPetre_eL! 



20071 ). The interpretation of the abundances, however, 
may be complicated by imperfect equilibration between 
the electron an d ion tempera t ures i n hot regions with 
T ~ 10^ K (jYoshida et al\ . Hool), unless collective 
plasma processes, such as scattering off Alven waves, are 
adequate for driving the electrons and ions into thermal 
equilibrium. 

Predictions for low ionization species are much more 
difficult. Observations suggest sizes for C ii, Si ii, 
Mg II and Al ii systems of typicall y below a kiloparsec , 
down to sizes of only about a parsec ( Simcoe et aLl . l206^ . 
Such small, dense clumps are a challenge for any cosmo- 
logical numerical simulation to produce. 



5. The flux power spectrum 

In addition to the direct flux and line parameter mea- 
surements, the models may also be tested by appealing to 
higher order moments of the flux distribution. The one 
which has received the most attention is the flux power 
spectrum, as a vehicle for using the Lya forest to con- 
strain cosmological models. Two broad approaches have 
been pursued. In one, measurements of the flux power 
spectrum are inverted to obtain a meas urement of the un- 
derlying d ark m atter power spectrum ( Croft et oil . Il998l 
119991: lHuil . ll999l : In usser and Haehneltl . [T99 ^ In the sec- 
ond, the flux power spectrum itself is used directly as a 
statistic to be matched by models ( Mandelbaum et all . 
120031: iMcDonald et aLl . |2000L l2005d) . 

The basis of the inversion methods is the rela- 
tion between the ID and 3D matter power spectra 
(jLumsden et al\ . Il989t ) 



poo 

Al^ik) = k dy 



and 



A^(fc) 



,d_ 

dk 



Al^ik) 



(109) 



(110) 



where A^p(fc) is the ID mass fluctuation spectrum, re- 
lated to the matter spatial correlation function ^ (x) by 



kx 



Al^yikx) cosikxx). 



(Ill) 



Eq. (|109p shows that A^j-, depends on all small scales, 
well into the nonlinear range of density fluctuations. Be- 
cause the flux along the line of sight to a QSO probes 
the ID density field of the IGM, it might appear that 



http:/ /www.rssd.esa.mt/mdex.php?project=XEUS 
http:/ /constellationx. nasa.gov 
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any predictions require a computation of the full nonlin- 
ear extent of the power spectrum to all small scales. Even 
the highest resolution simulations a re unable to achieve 
convergence on A^^ even at small k ( Meiksin and Whitd . 
I2OO4) . The statistics of the ID power spectrum are also 
notoriously non-gaussian. 

Such stringent demands for interpreting the flux power 
spectrum are avoided for two reasons: 1. the width of 
the absorption features acts to filter the 3D density field, 
erasing the smallest scale, nonlinear contributions to the 
ID flux power spectrum, and 2. the inversion Eq. (jllOp 
at low k depends only on the behavior of the ID flux 
power spectrum at low k. Numerical convergence on the 
ID flux power spectrum A|,(fc) at small k is achievable 
with sufficiently fine spatial resolution and a sufRciently 
large box size. A comoving box size of 25/i~^ Mpc with 
a force resolution of 30 — 60h^^ kpc (comoving) gives a 
flux power spectrum convergent to 10% on scales k < 
O.Olkm^^s at 2 : > 3, with more severe re quirements at 



lower redshifts ( Meiksin and Whitd . 120041) . The greater 
ease of convergence for the flux power spectrum, however, 
pays the price of insensitivity to the underlying matter 
power spectrum. As a result, a large number of spectra 
are required to measure the matter power spectrum to 
high acc uracy. A precision of 5% requires some 10'' lines 
of sight (|Zhan et aLl . l2005f) . 

Inver sion methods have been furt her developed and ap - 
plied bv lMcDonald et al\ (|2000t ) and lCroft etaJ] (|2002bh . 
who assume the ID flux power spectrum is propor- 
tional to the ID matter spectrum up to a scale de- 
pendent bias factor A% (k ) = & ^(fc)A? n(fc), and by 
iNusser and Haehneltl (|l999l |2000[ ) and IZaroubi et d\ 
( 20061 ). who convert local optical depth measurements 
to local density estimates, assuming a power-law relation 
between the two. 

In principle, the flux autocorrelation function may also 
be used to test models. Simulations, however, are found 
to converge even more slowly on this statistic, with con- 
vergence to better than 10% for separations exceeding 3 % 
of the box size problematic ( Meiksin and Whitd . 1200 . 

Methods that rely on direct comparisons between the 
measured and predicted flux power spectra require a large 
grid of simulations to cover the range in expected cosmo- 
logical model parameters. Because of the considerable 
computational expense involved in running fu ll hydro- 
dynamical simulations. iMcDonald et all (|2005b( ) rely in- 
stead on a large grid of HPM pseudo-hydrodynamics sim- 
ulations (see § IVI.A.2P , calibrated by a few full hydrody- 
namics runs. These are computed using 256^ cells in box 
sizes of 10h~^ Mpc (comoving) or s maller. This is well 
below the box size recommended by iMeiksin and Whit3 
( 20041 ). who also find that the prediction for A|,(fc) is sen- 
sitive to the method of smoothing at the several percent 
level near the peak and beyond (towards smaller values 
of k). In a comparison of full hydrodynamics simulations 
with HPM simul ations in 30h~^ Mpc (comoving) boxes, 
IViel et al\ (|2006l) find HPM simulations with many more 
grid zones (1200^) than particles (2 x 200^ for gas and 



dark matter) are required to converge to within 5% of the 
full hydrodynamics derived flux power spectrum on scales 
k < 0.03 km"^ s at z = 3 and 4. For 400^ grid zones, the 
discrepancy is several percent and scale-dependent. At 
z — 2, convergence is even poorer, with disagreement as 
large as 20-30%, as shocks play an increasingly impor- 
tant role in heating the gas, an effect which an HPM 
simulation cannot account for. 

A large number of physical systematic effects may 
also affect the f lux power s p ectrura . The impact of 
galact ic winds ( Croft e^"aZI . l2002al : iMcDonald al . 
2005al) and ionization radiation fluctuations (ICroftl . 
20041: ICroft et ad Il999l : iMcDonald eloZI . l2005a : 
Meiksin and White . |2004[ ) are generally found not to 



have much effect on the inferred matter power spectrum, 
except possibly at z = 5 — 6 when the attenuation 
length to ionizing photons becomes so short few sources 
contribute to the local UV background. 

More of a concern are known systematics that af- 
fect measurements of the Lya forest power spectrum. 
The effect of the wings of Damped Lya Absorbers 
and absorption by interve ning metal absorption systems 
(|Croft et a/.l . l2002bl . [l999h must be corrected for to obtain 
high precision measuremen ts of the flux po wer spectrum 
( {McDonald et al. 2005^. [McDonald et 2 I. (2006) pro- 
vide an estimate for the influence of the damping wings 
from DLAs. They also correct for metal absorption, 
uncertainties in sky subtraction, and calibration errors. 
They find clear evidence for contamination by Si iii, and 
devise a method for removing its influence. 

Any measurement of the flux power spectrum and com- 
parison with models also require a careful estimate of 
the QSO continuum. An optimal method for estimating 
the QSO continuum in the conte xt of flux power spectra 
measurements is formulated bv iHui et al\ ( 200ll ). The 
flux power spectrum was measured from high resolution 
spectra of 27 QSOs bv lKim et al\ (|2004allbl ). The largest 
and most precise set of measurements by far is that made 
using a sample of 3035 QSO spectra fro m the Sloan Digi- 
tal Sky Survey (jMcDonald et all\200^ . The evolution of 
the measured flux power spectrum is shown in Figure [T8l 

Despite their differences, the various approaches have 
yielded reasonably good agreement on the inferred con- 
straints of cosmological param eters. Using the inversion 
method, IZaroubi et al\ ( 20061 ). assuming n = 1 for the 
matter power spectrum, obtain cg — 0.92(r2„i/0.3)~°'^. 
A joint analysis with the WMAP CMB data and a prior 
on the Hubble constant yields erg = 0.92 ±0.04 and Sl„i = 
0.30 ± 0.05. Combining the Ly a flux power sp e ctrum 
measu rements from the SDSS ([McDonald et all l2006l 
l2005bD with CMB, superno vae and galaxy clustering con- 
straints, ISeljak et all (|2006l ) find the Lya flux power spec- 
trum increases erg = 0.80±0.03 to erg = 0.85±0.02. They 
obtain an index of n = 0.965 ±0.012 with no evidence for 
a running spectral index, dn/d\nk = — 0.020±0.012. In a 
sep arate analysis of the measu red SDSS flux power spec- 
tra, IViel and Haehneltl (jlooi) obtain as = 0.91 ± 0.07, 
= 0.28 ± 0.03 and n = 0.95 ± 0.04. Using pure 
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k [(km/s)-'] 

FIG. 18 Comparison of the best fitting model (lines) and 
the mean power spectrum measured from the Lya forest in 
3035 SDSS QSO spectra (points with error bars). The model 
includes the effect of Si iii absorption. The curves show a 
sequence in redshifts, increasing from z — 2.2 at the bot- 
tom to ^: = 3.8 at th e top, in steps of Az = 0.2. From 
iMcDonald et al\ (|2006l 'l. (Figure reproduced by permission of 
the A AS.) 



dark matter simulations, iDesiacaues and Nussed ( 20051 ) 
obtain as = 0.7 - 0.9. 



The flux power spectra measurements tend to pull 
the normalization of the matter power spectrum to 
larger values compared with estimates based on other 
methods. This is consistent with the larger values 
of (Tg — 0.9 — 1 preferred by cosmological models 
that match measurements of the p ixel flux distribution 
function and WMAP con straints (|Meiksin et all l200ll: 
iMeiksin and Whitel . [200l . On the other hand, the model 
comparisons all neglect radiative transfer effects. Al- 
lowing for an inverted equation of state over the red- 
shift range 2 < z < 3, in which underdense regions are 
war mer than overdense ones, lowers the allowed range for 
as ( Bolton et aLl . 120081 ). There is growing evidence that 



the equation of state is inverted, or at least not as mono- 
tonic as the power sp ectra and pixel flux distribution 
analy ses have assumed ( Becker et a/.l . l2007tlBolton et all . 
I2008D . Simulations including radiative transfer show a 
broad spectrum of temperatures at a given density, in- 
cluding hot under dense regions, in scenario s with late 
He II reionization ijTittlev and Meiksinl . |2007| ). 



VII. REIONIZATION 

A. The epoch of reionization 

The reionization of the IGM is one of the principal 
unsolved problems of cosmological structure formation. 
Three epochs of reionization of the IGM may be identi- 
fied, one for hydrogen and two for helium. Because the 
same radiation sources capable of photoionizing hydro- 
gen would likely also photoionize He i to He ii, only the 
epochs of hydrogen reionization and full helium reion- 
ization to He III are likely to be distinct. As discussed 
in §§ III. B. 51 and IIV.C.2I above, there is strong evidence 
full helium reionization occurred at z ~ 3.5. The most 
likely He ii reionization sources are QSOs. The epoch 
of hydrogen reionization is much less certain. The topic 
attracted considerable attention following measurements 
of the IGM Lya optical depth in high redshift QSOs and 
CMB measurements of intergalactic Thomson scattering. 
Unlike full helium reionization, the sources of hydrogen 
reionization may well be the first sources of radiation to 
turn on in the Universe, ending the Dark Ages that fol- 
lowed the recombination of the primordial hydrogen and 
helium plasma. Because of the numerous recent reioniza- 
tion reviews (see § HI , I concentrate here only on some of 
the most salient aspects. 

The Lya optical depth of the IGM measured in the 
spectra of several z > 5.8 QSOs discovered by the 
Sloan Digital Sky Surve y has show n that the IGM was 
reionized by z ~ 5.7 (|Fan et all l2001ai l2000l) . but 



the Ly g optical depth rises rap i dly towards higher red- 
shifts (iBecker a/.l. l2007ll200lUDiorgovski a/] .l200ll: 



iFan etall |2004 l2006bl : iPentericci et all I2002D . By z 

6, the spectra become so black shortward of Lya , usually 
only lower limits to the optical depth may be obtained, 
even using the higher order Lyman resonance transitions. 
The rapid rise is suggestive of an approach to the EoR. 
This interpretation is not unique, as only a small neutral 
fraction is adequate for providing a large Lya optical 
depth by Eq. (fT51) . A QSO at z c:^ 6.4 moreove r shows 
transmitted flux at z ~ 6.1 (jWhite et all l2003l 12005). 
Gaps, however, are expected as H ii re gions develop and 
evolve prior to their complete overlap ( Furlanetto et all . 
12004 IWvithe and Loebl . l2006^ . 

The interpretation of the IGM transmission measure- 
ments at z > 6 as evidence for the approach to the EoR 
is further confounded by WMAP measurements of the 
intergalactic Thomson optical depth 



dzHeCaT 



(112) 



0.004 



i^m{l + Zuf + 



1/2 



n^{i + zif + n^ 



1/2 



where ctt is the Thomson cross section, ne(0) is the co- 
moving electron density (assumed constant and without 
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full helium reionization) , and a flat universe is assumed. 
The 5-year WMAP polarization data yi eld a marginal- 
ized optical depth of = 0.087 ± 0.017 ( Dunklev et all 
I2OO8I ) up to the EoR, corresponding to an EoR redshift 
of Zu = 11. 0± 1.4, assuming reionization was a single, in- 
stantaneous event throughout the Universe. The 2a and 
3g lower limits are Zu > 8.2 and z„ > 6.7, respectively 
(iDunklev et a/.l . |2008| ) . These values are lowered slightly 
to Te = 0.084 ± 0.016 and z„ = 10.8 ± 1.4, allowing for 
additional astrophysical constraints on the preferred cos- 
mological model from Type la supernovae a nd Baryonic 
Acoustic Oscillations ( Hinshaw et aZ.I . [20081 ). 



IShapirol . Il986l : IShapiro and Girouxl . Il987^ 



It is worth noting that a non-negligible contribution to 
the total Thomson optical depth will arise from the IGM 
prior to the EoR because of the residual electron den- 
sity frozen out after the Recombination Era. A residual 
fraction of as little as (2 — 4) x 10~'^ over the redshift 
range 7 < z < 500 ( Page et a/] . l2007t) will produce an 
additional optical depth of Are = 0.005 — 0.01, using 
Eq. (112). Careful modeling of the contribution to the 
polarization signal from residual electrons between the 
Recombination Era and the Epoch of Reionization is re- 
quired to estimate the contribution due to reionization 
sources alone. 



B. The growth of cosmological H n regions 

The reionization of the Universe occurs as a conse- 
quence of the overlap of H 11 regions driven into the IGM 
by discrete radiation sources. Unlike the H 11 region 
around a star in a galaxy, the growth of cosmological 
H II regions will generally not terminate in a Stromgren 
sphere, defined as the surface within which the rate of 
radiative recombinations balances the rate of production 
of ionizing photons by the source. As the ionization front 
(I-front) of the H 11 region advances, the expansion of the 
Universe reduces the density of the gas into which the I- 
front moves. Unless the I-front encounters a large density 
inhomogeneity, it will travel unimpeded, lagging furt he r 
and further behind the Stromgren radius l|Shapird . [l98^ . 
until either the source dies, or the I-front overlaps with 
another I-front produced by another source or sources. 

The advance of the I-front of an isotropic cosmologi- 
cal H II region into a m ulticompon ent IGM i s described 
by (jPonahue and Shulli , ,1987, : Meiksin and Madau . ,1993; 



Excepted are scenarios invoking r eionization by decaying or an- 
nihilating dark matter partic les l[Ma^ell^^taZ ]. 120061; ISciamal . 



119821. 1199a IScott et all ll99lD . iMapelli et al\ bo06[) find that 
while some popular dark matter candidates may cause partial 
reionization, none are able to fully reionize the IGM prior to 
2 = 6. 
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'S'(0)-^7rrfnH ^ CifeH,iaB{Ti), 



(113) 



where r/ is the proper radius of the I-front, nn = 
J2i fi''^H,i is the mean hydrogen density of the IGM av- 
eraged over the various components with volume filling 
factors fi, H is the Hubble parameter, S'(O) is the emis- 
sion rate of ionizing photons from the central sources, 
Ci — fiTi^ i/'^H is partial recombination clumping fac- 
tor of the gas from component i, fe¥i,i is the number of 
electrons per hydrogen atom in component i (which will 
contain helium and possibly metals), and asiTi) is the 
Case B radiative recombination rate to the n = 2 level 
of hydrogen, allowing for ionizations from the ground 
state by the diffuse radiation emitted following recom- 
binations to the ground state (see § IIV.A.II above) . If 
the diffuse radiation is treated separately, then Case 
A recombination should be used, and a second set of 
partial ionization clumping factors must be introduced 
to represent the volume-averaged product of the local 
io nization rate and H i density in each subcomponent 
i ( Kohler et all l2007t) . If the electron fractions feHA 
and temperatures of the various components are approxi- 
mately the same, then fen,iOtB{Ti) may be approximated 
as feKCtsiT) and it is convenient to introduce a total 
recombination clumping factor over all subcomponents, 
C(z) = J2i Ci{z). The recombination term in Eq. (|113p is 
then {ATr/3)r^nf^CfeHaB{T). (If the diffuse radiation is 
treated separately, then a similar total ionization clump- 
ing factor may be similarly defined.) 

Eq. (|113p is an approximation for which it is assumed 
the I-front, within which the gas transforms from neutral 
to highly ionized, is narrow compared with the radius of 
the H II region. Is also assumes the I-front propagates 
rapidly compared with the sound speed, so that the ef- 
fects of gas motion induced by the pressure within the 
H II region may be neglected. It assumes as well that 
the I-front propagates slowly compared with the speed 
of light, which formally breaks down in the early stages. 
Near luminal expansion may be approximately described 
by replacing Airrju-a on the LHS of Eq. (I113|) by 47rrrnH+ 
S'/c where S' is t he RHS of Eg. (fTlB (IShapiro et al\ . 
I2OO6I : I White et al. I. I2OO3I : IWvithe and Loebl . l2004^ . This 
will ensure drj/dt < c. 

A key uncertainty in the rate at which an H 11 region 
expands is the recombination clumping factor. Because 
the clumping factor may be dominated by rare very dense 
structures, establishing its value from numerical simula- 
tions is not straightforward. In particular, systems that 
produce ionizing photons should be excluded. The dy- 
namical response of clumped structures to an impinging 
I-front must also be computed, which requires radiative 
transfer and high spatial resolution. The degree to which 
a dense cold clump bound by a minihalo impedes the 
propagation of an I-front in general depends on the mass 
of the clump and the intensity and spectral shape of the 



50 



incident radiation field ( Iliev et oil . l2005t IShapiro et al\ . 
I2004D . A wide range of values for the clumping factor 
are available in the literature. On th e basis of high res- 
olution simulations, llliev et ^ (I2007D find the recombi- 
nation clumping factor is well-fit over the redshift range 
6 < z < 30 by 



C{z) ~ 26.2917exp[-0.1822z + 0.003505z2 



(114) 



for models usin g the preferred WMA P 3-year cosmologi- 
cal parameters ( Spergel et al\ . l2007l ) . 

The evolution in the filling factor of ionized gas as the 
H II regions grow may be quantified by the porosity pa- 
rameter Qhii, which is the product of the volume ion- 
ized per source of a given ionizing photon luminosity, 
i ntegrated over the l umin o sity function of the source s 
([Donahue and Shulll . Il987t iMeiksin and Madaul . Il993l ). 
For Qhii > 1, the IGM becomes conipletely ionized . 
Its gr owth is governed by ( Madau et al\ . Il999l : iMeiksinl . 
I2005D 
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where nH(0) is the comoving number density of hydro- 
gen, ns{z) is the production rate of ionizing photons per 
comoving volume, and ti-cc{z) = l/[nH(z)/eHaB(T)] is 
the recombination time at the mean electron density at 
redshift z. If the comoving emissivity of the sources may 
be approximated as (z) = (i^/i^^)""® (1 + z)'''^ , the 
ionizing photon production rate is 
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To solve Eq. (|115p . it is convenient to define rH,Hi(z) — 
C(z)tH(z)/trcc ^ 0.02C(z)(l-|-z)3/2, which is the Hubble 
time, tjj = 2/3H{z), measured in units of the recombina- 
tion time, normalized here at T = 10^ K and including 
the recombination clumping factor. It is equivalent to the 
mean number of recombinations per hydrogen atom over 
a Hubble time. Similar definitions may be introduced for 
the reionization of He i, rH,Hoi — 0.03C(z)(l -I- z)'^/^ (as- 
suming all the ionized He i is in the form of He ii), and of 
He II, th.HcII — 0.2C(z)(H-z)^/2. Adopting the clumping 
factor from Eq. pi4p gives nearly constant values for the 
various th over 6 < z < 8 of th,hi — 4.7, rH,Hoi — 5.3, 
and th,hoII — 32. For constant th,hi, Eq. (jllSp may be 
directly integrated, yielding 
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where Zj is the turn-on redshift of the ionizing sources, 
and the assumption z^ 3> z has been made to obtain the 



final expression. Similar expressions apply for He i and 
He II reionization, on replacing th.hi by the appropriate 
form for He I or He 11, using the production rate of ioniz- 
ing photons for He I or He 11 for ng, and replacing nn^O) 
by nHc(O). In principle, a recombination clumping factor 
appropriate to each species must by used since the depth 
of self-shielding in clumps will differ for different species. 
To date, simulations have focused on the clumping of 
ionized hydrogen gas only. 

The (comoving) production rate of ionizing photons 
from the UV background emissivity estimate of Eq. (|106p 
follows from Eq. (|116p . 



As["-±^yi + zr hphs-'Mpc-^ 



where As = 1.3+g;^;tg x lO^^, 7 
0.6 ln(A5/1.3 X lO'^^), and as is t he spectral inde x at fre- 
quencies above the Lyman edge ( Meiksinl . |2005[ ) . (This 
neglects any contribution to the required emissivity from 
the diffuse recombination radiation from the IGM. If the 
IGM contributes 30% of the emissivity. As should be re- 
duced by the factor 1.3.) Substitution into Eq. (|117p and 
extrapolating Eq. I|118[) to high redshifts shows that the 
Universe would be reionized by Zri — 5.5 for umg = ols = 
0.5, but as late as Zri ~ 3.8 for q;mg — ols — 1.8. Reion- 
ization by Zri 6, as suggested by observations, requires 
either a hard spectrum for the sources or an increase over 
the evolution rate of Eg. (llOSp in the comoving emissiv- 
ity at z > 6 (iBolton and Haehneltl . 120071 : iMeiksinl . l2005t 
iMiralda-Escudel l2003f) . 



(118) 
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C. Sources of reionization 



1. Galaxies 



High redshift star-forming galaxies are generally con- 
sidered the most likely candidates for the reionization 
of the Universe. The discovery of high redshift galax- 
ies in deep surveys has sparked much lively debate over 
whether the population of galaxi es that reioniz e d the 
Universe has n ow b e en ide n tified ( Bouwens et al\ . 
20061 . 12004 



Stiavelli et al 



2006 . 



2004 



2004; Bunker et al 

lYan and Windhorst . 12004) . The uncertainties involved 
preclude any definite decision. 

Consider th e (comoving) luminosi ty density at 5.5 < 
z < 6.5 from iBouwens et al\ ( 2006t ). corrected upward 
by 50% to allow for internal extinction, of piaso — 
IQi^ "' WHz~^ Mpc""^. The conversion rates between 
Li,(1350 A) (in WHz~^), and the production rate of ion- 
izing photons for a Salpeter Initial Mass Function with 
-^'^lowor — 1 Mq and A/uppor = 100 Mq and solar mctal- 
licity are 

ATg^Hi = 1022"phs~iL^(1350A); HI 
iVs.Hoi = 10^^-^phs"iL^(1350A); Hcl 
A^s,Hoii = 10^^%hs-iL^(1350A); HeH (119) 
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( Leitherer et aZj . Il999l ). For a metallicity of 1/20 solar, 
the rates are about 0.1 dex higher for H i, 0.2 dex higher 
for He I, and 0.4 dex smaller for He ii. A s i milar con- 
version for H i is in ferred from iMadau et all ( 19981 ) and 
iMadau et all ( 1999) for an alternative Salpeter model 



with Mio 



0.1 M0 and Afuppor = 125 Mq. Assuming 



that all the ionizing photons escape, applying these rates 
to the luminosity density at 5.5 < z < 6.5, Eq. (|117p (and 
its variants for He i and He ii reionization, assuming the 
same recombination clumping factors, which may in fact 
differ) shows that the galaxies are able to reionize the H i 
by Zri — 12 and the He i by Zri = 21, while they are un- 
able to reionize the He ii even by the present epoch. The 
results for H i and He i, however, are very sensitive to 
the uncertainty in the counts and the escape fractions. 
If only 10% of the ionizing radiation escapes, then the 
reionization epochs for H i and He i are reduced dramat- 
ically to Zri = 1-8 and 3.7, respectively. Reionization of 
H I by z = 6 requires /esc > 0.6, while for He i, f^sc > 0.3 
is needed. Allowing for the upper limits on the counts 
eases these constraints, as would lower metallicities for 3. Other sources 



ulation of very low luminosity Active Galactic Nu- 
clei (AGN) (iBolton and Haehneltl . 120071: iMeiksinl . 120051: 
ISrbinovskv and Wvithd . 2007^ . The estimated produc- 
tion rate of ionizing photons from accretion onto massive 
black holes in the Univer se is comparable to the rate from 
nucleosynthesis in stars (jMeiksinl . l2006ct ) : low luminosity 
AGN may still await discovery at high redshifts. By con- 
trast, QSOs likely fully ionized helium by z ~ 3, and 
possibly as early as z ~ 5 if they have hard sp ectra ex- 
tendi ng down to the He ii photoelectric edge ( Meiksinl . 
I2005D . 

An interesting alternative class of reionization sources 
is a possible population of miniquasars, systems pro- 
duced by intermediate mass black holes with masses 
200—1000 Mq growing by accretion onto seed black holes 
produ ced by the co l lapse of very massive Population III 
stars (|Madau et all . \2004 . Such sources could plausibly 
reionize the Universe by z ~ 15. 



the stars ( Tumlinson and Shulll . [2OO0I ) or an Initial Mass 
Function heavi l y tilte d to very massive (~ 300 Mq) stars 
( Bromm et a/.l . BoOll ). On the other hand, a decline in 
the comoving luminosity density of the sources at z > 6 
would make the case that these sources reionized the Uni- 
verse less secure. While the population of sources that 
reionized the Universe may have been discovered in high 
redshift surveys, a definite conclusion is still premature. 
There is some evidence for a population of low- luminosity 
star-forming galaxies at redshifts as high as z ~ 8 — 10 
that may contribut e substantia l ly to th e total budget of 
reionizing photons ( Stark et all . l2007"al ) . 

Several semi-analytic estimates for the expected con- 
tribution of galaxies to the reionization of both hy drogen 
and helium have been made ([Haiman and BrvanJ, 12006 ; 



Various alternative sources have been proposed 
for the reionizati on of t he IG M, including globu- 
lar star clusters (|Ricottil . I2004D . massive Population 
HI stars more generally, with their formation pos- 
sibly triggered early by pr imordial magnetic fields 
( Tashiro and Sugivama . 20061) . and an early generation 
of black hol es feicotti and Ostrikeil 120041 : iRicotti et all . 
120051: IVenka t csan et al, 2001^ While these mechanisms 
are intriguing and have received a variety of theoretical 
investigations, observational evidence for any of them is 
as yet wanting. The possibility that sources other than 
galaxies reionized the IGM, however, should not be dis- 
counted: the Universe can surprise. 



and neiium nave be e n ma.ae ([iia iman and aryan. 
Haiman and H older. 2 0031 : IShuU and Venkatesaij . 



20081: 



Venkatesan et al.. . .200^) . The possible halo masses of the ^- 21cm signature of reionization 



objects that dominate reionization cover a broad range, 
including systems as small as minihalos with masses 
-M^haio < 10^ M© and virial temperatures Tvir < 1000 K 
which would have so far evaded detection. Definitive pre- 
dictions, however, are hampered by the need to incorpo- 
rate several poorly determined quantities, including the 
fraction of baryons in a collapsed halo that are converted 
into stars, the number of ionizing photons produced by 
the stars, and the escape fraction of the ionizing radiation 
from the galaxies. Recent estimates suggest that the ef- 
ficiency of ionization per galaxy cannot be too high, or a 
Thomson optical depth would be produced that is too 
large compared with CMB polarization me asurements 
( Haiman and Brvanll2006l : IShull and Venkatesan. . ,2008i) . 



An altogether novel method for discovering the Epoch 
of Reionization is through the detection of the 21cm sig- 
nature from the neutral IGM before and during the reion- 
ization process. The emission or absorption of 21cm ra- 
diation from a neutral IGM is governed by the spin tem- 
perature Ts of the hydrogen, defined by 



no 



3 exp 



(120) 



2. QSOs 

The emissivity of QSOs lies well below the require- 
ments for reionization, unless there is a large pop- 



where no and ni are the singlet and triplet n = 1 
hyperfine levels and = h-pUio/k-Q ~ 0.07 K, where 
1^10 is the frequency of the 21cni transition. In the 
presence of only the CMB radiation, the spin temper- 
ature will come into equilibrium with the CMB on a 
timescale of T,/(Tcmb^io) =i 3 x 10^(1 + z^^ yr, 
where Aiq — 2.85 x lO"^'^ s^^ is the spontaneous de- 
cay rate of the transition. In this case, the spin tem- 
perature will be the same as the CMB temperature, 
and no emission or absorption relative to the CMB 
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will be detectable. A mechanism is needed to decou- 
ple the two temperatures. Two mechanisms may cou- 
ple the spin temperature to the kinetic temperature of 
the gas. One is the co llisions between hydrogen atoms 
( Purcell and Fieldl . [1956) . The collision-induced cou- 
pling between the spin and kinetic temperatures is dom- 
inated by spin-exchange between the colliding hydro- 
gen atoms, with possible significant contributions from 
electron-hydrogen collisions and proton-hydrogen colli- 
sions, depending on the ionizat ion f raction and tem- 
perature ( Allison and DalgarncJ Il969l: [b algarng. 1967 : 
Furlanetto a.nd Furlanettol. 2007: H irata and Sigurdson . 
20071: ISmithl . Il966t IZygel^an. . ■2005t ). The secon d, less 
obvious mechanisni , is sca ttering by Lya photons ( Fieldl . 
119581 : IWouthuvseiH . Il952^. The Wouthuysen-Field ef- 
fect mixes the hyperfine levels of neutral hydrogen in 
its ground state via an intermediate transition to the 
2p state. An atom initially in the n = 1 singlet state 
may absorb a Lya photon that puts it in an n = 2 
state, allowing it to return to the triplet n = 1 state 
by a spontaneous decay. The combined atomic collision 
rates are generally too small at the average IGM densi- 
ties at the redshifts of interest, although collisions may be 
i mportant in o yerde nse regions, (5p/yO^30[(l -|- 2:)/10]~^ 
( Madau et a/.l . [l997[ l. Instead the dominant mechanism 
at detectable redshifts for the near future is likely to be 
Lya scattering. 

A patchwork of either 21cm emission, or absorption 
against the Cosmic Mi crowave Background, will result 
( Hogan and Reesl . ll979l ). Large radio telescopes like LO- 
FAR, the MWA, PAPER, or the SKA offer the prospect 
of measuring this signature, and so detecting the transi- 
tional epoch from a dark universe to one with light. The 
spin temperature of the neutral hydrogen is 



Ts 



TcMB + UaTa + VcTr 



where 



Va = 



Pio 



and 
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are the Lya and coUisional pumping efficiencies, respec- 
tively. Here Pio is the indirect de-excitation rate of the 
triplet state via absorption of a Lya photon to the n = 2 
level, and Ts ^ T* is assumed. Consideration of the net 
transition rates between the various hyperfine n = 1 and 
n = 2 levels shows that the 1 — > transition rate via 
Lya scattering is related to the total rate 



Pa = c 



dv a I, 



(123) 



where a^, is the resonance line cross section given by 
Eq. and U i, is the ra diation energy density, by 

Pio = 4P„/27 (|Fieldl . Il958l ). The quantity is the 
color temperature of the radiation field, given by the har- 
monic mean = l/{T~^{v)) of the frequency-specific 
temperature Tu{v) — —{hp /kB){d\nu„ /dv)~^ , weighted 



by u^ipvia,u) (jMeiksinl . [2"006bl ). The factor do is the 
rate of collisional de-excitation of the triplet level. In the 
absence of collisions, there exists a critical value of Pa 
which, if greatly excee ded, would dr ive Ts Tk- This 
thermalization rate is ( Madau et a/.ljl997 ) 
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It corresponds to a Lya intensity of ~ O.&S^^ x 
lO^^"* Wm~2 Hz~-i sr~i, where Sa is a suppression fac- 
tor of order unity that takes into account the ef- 
fect of recoils in an expanding medium on the shape 
of t he intensity near the c e nter of the line pro- 
file (IChen and Miralda-Escudel |2004 iFurlanetto et all 
I2OO6 ). The value of the intensity corresponds to a com- 
parable number of Lya photons as hydrogen atoms at 
z = 8, and so is comparable to the inte nsity expected 
from sources able to reionize the Universe ( Madau et al\ . 
I1997D . It is this coincidence that promises that the 
Wouthuysen-Field mechanism may reveal the IGM prior 
to its complete reionization. 

To illustrate the basic principle of the observations, 
consider a patch of neutral material with spin temper- 
ature Ts 7^ JcMB, having an angular size on the sky 
which is large compared with a beamwidth, and radial 
velocity extent due to the Hubble expansion which is 
larger than the bandwidth. Following an analysis sim- 
ilar to the derivation of the Gunn-Peterson optical depth 
(Eq. [15]), and taking into account the reduction of the 
absorption coefficient due to stimulated emission by the 
factor 1 - e-'^*/^s ~ T^/Ts for T^<s:Ts, the intergalac- 
tic optical depth at 21(1 -I- z) cm along the line of sight 
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which will typically be much less than unity. Here, cchi 
is the neutral fraction, 6 accounts for fluctuations about 
the mean hydrogen density, nni = ''^HI(1 + S), and a 
factor of 1 /4 is introduced since only the hydrogen atoms 
in the singlet state will take part in the absorption. A 
flat universe is assu med f o r the s econd form in Eq. (|125[) 



and ^hh'^ = 0.022. iFieldl ^dEM) applied Eq. (fT^ to a 



search for intergalactic emission or absorption against a 
background radio galaxy. No clear emission or absorption 
signature was detected, with an upper limit on the 21cm 
optical depth of r < 0.0075. For = 0.3, h = 0.7, and 
Ts = TcMB, this corresponds to a;Hi(l + S) < 1.5 along 
the line of sight. Adopting = 1 and S — would 
correspond to an upper limit on the baryon density of 
flh < 0.2 assuming the hydrogen was all neutral, showing 
that early measurements could already have ruled out a 
baryon-dominated universe closed by a neutral IGM. 
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In the context of radio measurements of the IGM 
planned today, an alternative observational strategy con- 
sists of cross-correlating measurements of patches sepa- 
rated in either angle or frequency. Possible detections 
include large isolated H i i regions, as would be pro- 
duced by the first QSOS (|Madau et all Il997h or, pos- 
sibly more straightforwardly, a global statistically aver- 
aged fluctuation signal from the IGM (,Toz zi et aL, 2000). 
Foreground contamination, particularly fluctuations in 
the ionosphere and discrete sources in the Galaxy, pose 
formidable challenges to detecting the 21cm signature 
from the Epoch of Reionization. But the development 
of novel detection strategies suggests the challenges may 
not be insurmountable. Discussio ns of the feasibility o f 
the measurements are provided by Wvithe et al\ ( 2005( ). 
iBowman et~ai\ (l2006l) iDatta et al\ (l2007l). and in th e 
reviews bv iFan et M (l2006al) . iFurlanetto et al\ (|2006t) . 
and lBarkana and Loebl ( 2007 ). which also discuss specific 
reionization scenarios. 



VIII. THE ABSORBER-GALAXY CONNECTION 

Ever since their discovery, intervening absorption sys- 
tems have held the promise of opening a new window on 
the evolution of galaxies by probing the gas in galaxies 
as it was converted into stars, and tracing the feedback 
of the stars on the gas as supernovae and stellar winds 
pollute their environment with metals. Despite long and 
hard searches, direct evidence relating intervening ab- 
sorbers to galaxies has been scant until only the past 
few years. Today, the absorber-galaxy connection is fast 
becoming a boom area in IGM studies. 



A. Galaxy-associated IGM absorption 

1. Association of DLAs with galaxies 

The motivation for the original systematic Damped 
Lya Absorber survey was the sea rch for the gaseou s 
disks of galaxies at high redshifts ( Wolfe et al\ . Il986f ). 
A large number of surveys have been cond ucted since, 
summ arized in the extensive review by IWolfe et al\ 
(I2005D . There are several properties of DLAs sugges- 
tive of galactic disk material: 1. They are the dom- 
inant re servoir of neutral hydrogen at high r e dshift s 
[z > 3) (jJorgenson et all l2006t iProchaska et all l2005t) . 



2. with the cosmic gas mass density decreasing and 



metallicity increasing with cosmic time (|Prochaska et all 
I2003D . 3. There is evidence for depletion of Fe and Si 
onto dust grains simil ar in strength to that in cold gas 
in the Galaxy (Pctitiean et a/.l . |2002( ). 4. There is photo- 
metric evidence for reddening by dust, with the amount 
of extinction increasing with iron column density simi- 
lar to th e trend found for in terstellar clouds within the 
Galaxy ( Vladilo et all I2006D . 5. Molecul ar hydrogen is 
present in about 15-20% of the systems ( Ledoux et all 



l2006bl ). with the metallicity and depletion factor increas- 
ing with the amount of H2 ( Ledoux et aU [20031 ) . 6. The 
metallicity of DLAs increases with the velocity width 
of the line profiles of the low-ionization metals, which 
presumably trace the gravitational potentials of the sys- 
tems, aiid_soJs_suggestiw of a mass-metallicity correla- 
tion (^ Ledoux "et 7. In addition to a cold low 
ionization phase, a warm phase is indicated by C iv and 
Si IV absorption systems, and a hot phase is indicated 
by O VI absorption detections ( Fox et aH . l2007l ). demon- 
strating a complex interstellar medium in DLAs much as 
is found in disk galaxies. 8. Using a method that deduces 
the star formation rate from the strength of C ii*A 1335.7 
as a measure of the amount of gas heating by the pho- 
toelectric dust mec hanism due to UV heating by mas- 
sive stars, IWolfe et ai. (,2003b. ) infer star formation rates 
in a presumed cold neutral medium in DLAs similar to 
that in the interstellar medium of the Galaxy. 9. The 
inferred global cosmic star formation rate in DLAs at 
z ~ 3 is comparable to that of Lyman B reak Galaxies 
at the same redshift (|Wolfe et a/.l . i2003al) . 10. A search 



for Lyman Break Galaxies at z > 1 associated with 
DLAs shows evidence for a significant spatial correlation 
between the two, with a cross-correlation length com- 
parable t o the autoco r relati on length of Lyman Break 
Galaxies (|Cooke et all l2006l) . 11. Optical counterparts 
to DLAs at z < 1 suggest the DLAs arise as part of large 
(radii of ~ 30 kpc), rotating H i disks (jChen etaU . 
I2005al) . 12. The incidence rate and column density dis- 
tribution of DLAs at z < 1 is comparable t o those of low 
redshift galaxies directly measured in H i ( Zwaan et all 
I2OO5D . 

Despite the circumstantial evidence favoring a galactic 
disk connection for DLAs, there is strong contrary evi- 
dence as well: 1. While they are the dominant reservoir 
of neutral gas, DLAs contain less than half the mass den- 
sity of present day stars. 2. The amount of neutral gas 
i n DLAs exceeds t he stellar mass in present day disks 

dmHIilai], iooi), which suggests they trace the disk 

component of galaxies, but their metal depletion pat- 
terns most often resemble that of warm disk clouds in the 
Galaxy or halo clouds (Ledoux e t al,, ,20Q2,) . The differ- 
ence in the UV metagalactic backgrounds at low and high 
redshifts may partially account for the apparent discrep- 
ancy. 3. The C II* inferred star formation rates suggest 
that a mass in metals should have b een produced a facto r 
30 larger than measured in DLAs ( Wolfe et aLl . l2003al ). 
4. The C II* inferred star formation rates so mewhat ex- 
ceed t he predict i on of the Kennicutt-Schmidt ( Kenniciittl . 
Il998l : ISchmidtl . Il959i) relation between the star for- 
matio n rate and gas mass in galaxies ( Hopkins et all 
I2OO5I ). Possibly the principal source of heat powering 
the measured C 11* cooling rates is star formation in cen- 
trally located, compact Lyman Break Galaxies rather 
than in the bulk of a spatially extended DLA system 
(IWolfe and Chenl . |2006[ ) . 5. DLAs are sufficiently abun- 
dant that if they obeyed the Kennicutt-Schmidt rela- 
tion, then 3% of the sky should be lit up with star- 
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forming DLAs at high redshifts, contrary to observations 
which find a comoving star for mation density a facto r 
30-100 smaller than predicted (| Wolfe and Chenl . |2006^ . 
6. Most puzzling is the rarity of successful follow-up op- 
tical identifications. While there is mounting circum- 
stantial evidence that DLAs and Lyman Break Galaxies 
are associated with the same systems at z > 1.6, only 
a han dful of DLAs show Lya in emission ( Moller et all . 
\2004) . Galactic counterparts have been confirmed in the 
optical a nd near-i nfrared for about a dozen DL As at 
z < 1.6 (|Chen et a l. 2005a; Chen and Lanzetta^, 120031: 
' l2003t ISchulte-Ladbeck et all 120041) . The 



iRao et ai. 



numbers are too few, however, to establish whether or 
not the galaxies repres ent a random samp ling of the main 
population of galaxies (IWolfe et al ]. l2005D . The nature of 
DLAs and their relation to present day galaxies remain 
an enigma. 

Simulations set in a ACDM cosmology reasonably well 
reproduce the number of DLAs and their evolution at 
high redshifts when dissipation due to cooling, star for- 
mation and subsequent supern ovae feedback and galac- 
tic winds arc included (Nag amine et all . l2004al ). The 
DLAs are found to arise in halos no smaller than those 
with total masse s of abo ut 10^h~^ Mq at high redshifts 
( Nagamine et all . l2004a[) . and almost all in dark mat - 
ter halos less than lO^^/i'^M© (jNagamine et aLl . [2007[) . 
Detailed comparisons between the simulated DLAs and 
observations, however, reveal some discrepancies. The 
predictions for the number of systems along a line of 
sight lie somewhat below the measured number. This 
may be due to overly compact DLAs in the simulations, 
which also show narrower velocity widths than measured. 
High resolution simulations with radiative transfer sug- 
gest that gas in tidal tai ls and quasi-filamentar y struc- 
tures also produce DLAs ( Razoumov et aLl . l2006[ ). which 
may help to reconcile the predicted and measured num- 
bers. The star formation rate in the simulations is con- 
sistent with the observational estimates in DLAs, how- 
ever an order of magnitude mo re metals are produced i n 
the simulations than measured ("Nagami ne aLl . l2004b[) . 
Simulations with radiative transfer suggest that the UV 
radiation from stars produced in DLAs and the ambient 
UV metagalactic backgroun d may regulate t he rate of 
star formation within DLAs ( Iliev et all . l2006l ). Detailed 
kinematic comparisons between the simulations and the 
observations have yet to be made. While the simulations 
have gone a long way toward accounting for the DLAs 
within the framework of the ACDM model for structure 
formation, much work still remains to demonstrate that 
this model correctly describes the origin of most DLAs. 



2. Association of metal absorbers with galaxies 

Optical and infrared searches for counterparts to 
metal absorption systems, such as Mg ii or C iv, 
have proven an effective means of discovering galax 



ISteidel et al. I. ll997lll994l . ll995l) . On the whole, the galax- 
ies appear to sample the general population of galaxies 
in the field, although there occur fewer dim blue galax- 
ies in the absorber-selected samples. The morphologi- 
cal and kinematic properties o f galaxies assoc i ated w ith 
Mg II systems are reviewed bv lChurchill et all (120051 ). 

The spectra of QSOs proxim ate to galaxies on t he sky 
probe galactic environments. lAdelberger et all (l2005f ) 



ies ( Bergeron and Boissel Il99lt iBouche et al 



s gaiax- 



find galaxies at 2 < z < 3 generally show very strong C iv 
absorption, A'^civ ^ lO^^'cm"^, within 40 kpc (proper) 
of the galaxy, and C iv absorbers with typical column 
densities of A^civ — lO^'^cm^^ are found as far out 
as 80 kpc. The key question is why the metals are 
near the galaxies. Galaxies are known to be associ- 
ated with large scale baryonic overdensities, at least as 
probed by the Lya forest. In the nearby Universe, the 
stronger absorbers (13.2 < log^g A'^hi < 15.4) are spa- 
tially correlated with galaxies, although with mixed views 
over the strength of the cross-correl ation compared with 
the galaxy autocorrelation function (iPenton et all . l2002t 
iRvan-W eber. '2006'; 'Wilman et^, |2007|). At higher red- 
shifts, lAdelberger et al.. (.20051 ) similarly find galaxies are 
associated with enhanced amounts of Lya absorption. 
If the IGM is polluted by metals, then the proximity 
of metal absorption systems with galaxies may merely 
reflect the correlation of H i with galaxies. Alterna- 
tively, the metals associated with galaxies may have been 
ejected by the galaxies in a wind following an episode of 
intense star formation. Either alternative is intriguing. 
The possibility that galaxies polluted their environment 
with metals suggests these metals may serve as tracers of 
the past star formation history of the stars in the galaxy, 
and so provide independent constraints on their forma- 
tion. On the other hand, if metal pollution is widespread 
throughout the IGM, it raises the question of whether the 
pollution was in situ, produced by some ubiquitous un- 
detected class of stellar systems like a population of high 
redshift globular star clusters, or are the collective ejecta 
of galaxies that underwent a superwind phase. In the 
latter case, it may mean that the evolution of the stellar 
population of galaxies was not isolated, but that galaxies 
"talked" with each other as their stars formed, entrain- 
ing metal enriched material produced by other galaxies. 
A galactic wind may also delay the onset of star for- 
mation in nearby low-mass halos by delivering t urbulent 
kinetic support to the gas ( Fuiita et all . |2004[ ). In ex- 
treme cases, a wind may eradicate the gaseous content 
of nearby low-mass halos through evaporation or shock- 
driven strippi ng, suppressing star forma tion in these sys- 
tems entirely ( Scannapieco et a/] . l2000t ). 

Evidence that hig h redshift galaxies do imp act on the 
IGM is provided bv lAdelberger et all (|2005l i. who find 
that about a third of the galaxies in their sample at 
2 ^; -2 ^; 3 show clearings of weak or absent H i absorp- 
tion out to distances of lh~^ Mpc (comoving). Simu- 
lations suggest, however, that the clearings are c onsis- 
tent with random fluctua tions ( Desiacaues et all . l2006t 
iKawata and Rauchl . |2007( ). The kinematic studies of 
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ISoneailal (|2006 h show that C iv systems at 2 < z < 3.5 
with A'^civ > 2 X lO^'^ cm"^ he in complexes having veloc- 
ity widths of up to 300 — 400kms~-^, suggestive of galactic 
outflows. Lower column density systems lie in complexes 
with substantially lower velocity widths of ~ 50 km s~^, 
suggesting they are not part of any galactic outflow, al- 
though they may be part of residual outflows that were 
initiated at much higher redshifts (Fujita et ai, 2004). 
The evidence for galactic outflows has prompted an in- 
creasing interest in the impact of galactic winds on the 
IGM, discussed next. 



B. Galactic winds and the IGIVI 

Galactic winds occur throughout the Universe, in a 
variety of environments, over a wide range of galac- 
tic masses, incl u ding starbursts in the local Universe 
( Heckman et al. . 1998*) and Lyman Break Galaxies at 
z ~ 3 (jPettini et ai. 200 2,) and Lya emitti ng systems 
associated with them (jWilman et all . |2005[ ). A com- 
prehensive r eview of galac t ic win ds in the Universe is 
provided bv IVeilleux oZ. (12005). Reviews specif ic to 
the IGM are provide d by lAguirre and Schavd (I2007D and 
iBertone et all (|2008l ). 

A natural mechanism for enriching the gas near galax- 
ies with metals is through winds. The possible enrich- 
ment by wind ejecta raises further questions: What sort 
of galaxies ejected the metals? When did they do it? 
How does the metal pollution relate to the star forma- 
tion history of the Universe? 

Geometric estimates of the impact of galactic winds 
may be made on the basis of the spatial corre l ations 
of the metals and galaxies. iPorciani and Madaul (I2005D 
find that metal-enriched bubbles of comoving radius ^ 
100 kpc expelled at z = 9 by a biased population of dwarf 
galaxies with a halo mass of lO^M©, polluting the region 
around Lyman Break Galaxies, would account both for 
the metal abundances measured at z = 3 and the mea- 
sured spatial cross-correlation between metals and Ly- 
man Break Galaxies. Using hig h resolut ion, high signal- 
to-noise ratio spectra, Scannapi eco et all (2006a, b) study 
the clustering properties of high (C iv and Si iv) and low 
(Mg II and Fe ii) ionization metal absorption systems. In 
conjunction with numerical simulations, they find that a 
uniform distribution of metals fails to reproduce the cor- 
relations measured. Modulating the metallicity by the 
density field similarly fails. They do find, however, that 
the C IV correlation function is well-reproduced by as- 
signing the metals to clumps confined to bubbles of co- 
moving radius 2 Mpc around systems with masses ex- 
ceeding IO^^Mq at z = 3. A similar model with a slightly 
larger bubble radius recovers the measured correlations 
of Mg II and Fe n at z — 0.5 — 2. Both models are de- 
generate with a biased high-redshift enrichment model 
with the metals confined to 2.5 Mpc (comoving) radius 
bubbles around 3 x 10^ Mq systems at z = 7.5. The re- 
sults are consistent with two scenarios, metals deposited 



by winds from massive galaxies at z ~ 3, or by winds 
from dwarf galaxies at high redshifts, z > 7. An early 
wind model is also favored by the near constant mass 
density associated wit h C iv systems over the broad red- 
shift range 2 < z < 6 ( Rvan- Weber et aTl . 120061 : fSongailal 
I2005D . Whether any of these scenarios are viable depends 
on whether wind models that support them may be phys- 
ically realized. 

Correlations between low redshift absorption systems 
and galaxies provide a potentially power ful means of in- 
ferring the origins of the metals. Tumlinson and Fans! 
(I2005D find that dwarf galaxies are able to reproduce the 
line-of-sight number density of O vi systems if they pro- 
duce winds extending out to about 200 kpc, while lu- 
minous galaxies would need to produce winds extending 
to 0.5 — 1 Mpc. Simil ar results are obtain ed from the 
O VI survey analysis of lStocke et al. Idlooii), who match 
O VI systems at z < 0.15 with galaxies in the field. They 
find the metal systems have median distances of 350- 
500 kpc (for h — 0.7) from galaxies, and 200-270 kpc 
around O.IL* galaxies. A contribution from O.IL* galax- 
ies is required to reproduce the line-of-sight number den- 
sity of O VI absorbers. 

While these results are consistent with the ejection of 
metals from the galaxies detected, this is not the only 
possible picture. The metals may have arisen at earlier 
times in smaller systems that eventually merged into the 
galaxies seen. Alternatively, they may have been ejected 
from much smaller systems, perhaps as small as globular 
star clusters or very faint dwarf galaxies, that are spa- 
tially correlated with the visible galaxies. The metals 
may also have been deposited very early in the history 
of the Universe, perhaps by blowout from the projenitors 
of present day galaxies in superwinds or possibly driven 
out by the impact of AGN that may have been active in 
the projenitors. 

As a galactic wind expands into the IGM, it will drive 
a shockfront that sweeps up the surrounding gas into 
a thin shell just behind the shock. In this limit, the 
thin-shell approximation may be made to describe the 
advance of the shockfront. Denoting the distance of the 
shockfront from the center of the galaxy by Rg , the mass 
of the expanding gas interior to Rs by Mg, the mean 
interior pressure by Pi, and enthalpy density by hi,, the 
external pressure and density by Pq and po, respectively, 
the velocity of the external gas relative to the galaxy by 
vh, the gravitational potential energy of the wind gas 
by W, and assuming mass injection from the center at 
the rate Min with velocity Vin, the equations describing 
the evolution of the shockfront radius and the interior 
thermal energy Ei, of the bubble powered by energy in- 
jected at the rate L[^{t) and coohng at the rate Lr„n\{t) 
are, in the thin-shell appro ximation (|Madau et a/.l , [200]] : 
lOstriker and McKed . fl988l) . 



dt 



(MgRs) = A7TR^^il-U)[iPb-Po)+PoVHiRs-VH)] 



-MgR^JX,l + W/Rs + MinV-n 



(126) 
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at 



(127) 



where isothermal jump conditions across the shock front 
have been assumed for the energy equation, presuming 
the wind expands into a photoionized medium. A clumpy 
interior is allowed for due to the presence of clouds with 
a volume filling factor /d and spatial density and cross 
section corresponding to a mean free path Ad, which en- 
ters as a drag term on the advance of the shockfront in 
Eq. (|126p . The interior net cooling Lcooi will be domi- 
nated by Compton cooling at high redshifts and two body 
atomic processes, which in turn depend on the gas clump- 
ing factor, and radiative heating from the galaxy or UV 
metagalactic background at lower redshifts. The second 
and fourth terms on the RHS of Eq. (|127p account for en- 
ergy loss due to work done by the shock on the surround- 
ing gas and on crushing interior clouds, respectively. 

The factors in Eqs. (fT26)) - (fT27|) highhght the various 
uncertainties in any galactic wind calculation. The rate 
of mass and energy injection by the stars and supernovae 
depend on the rate of star formation and the Initial Mass 
Function of the stars, which may have a lower upper 
mass cutoff in starbursts comp ared with normal galaxies 
( Leitherer and Heckmad . [T995h . or m ay be dominated by 
mass ive stars if the metallicity is low ( Bromm and Loebl . 
I2OO3I ). The cooling rate will depend on the metallicity 
of the bubble interior, its dumpiness, and on how well 
mixed the metals are. It will also depend on the temper- 
ature of the interior gas, which requires a specification of 
the amount of interior mass even though it is negligible 
compared with the mass in the shell. Most of the inte- 
rior mass results from evaporative losses from the con- 
tact surface between the hot interior and the cool shell 
( Madau et al\ . l200l[ ). The rate of evaporation depends 



on the degree to which thermal heat conduction may be 
suppressed by interior magnetic fields, although mechan- 
ical heat input balancing adiabatic expansion losses dom- 
inates the interior thermodynamics, so this may not be 
critical. The strength of the shock front may be sub- 
stantially dissipated as the shock overtakes small clouds 
and crushes them. Lastly, the progress of the shockfront 
depends on the environment of the galaxy, which will 
differ in different directions, so that the winds will in 
general be anisotropic. The acceleration of dense gas by 
hot rarefied gas and the motion of surfaces past one an- 
other as the bubble penetrates the IGM will also gener- 
ally be subject to Ray leigh- Taylor and Kelvin-Helmholtz 
instabilities, giving rise to the mixing of gas components 
with vastly contrasting densities and temperatures and 
the formation of eddies both at the surface of the bub- 
ble and in its interior as the shell fragments. Thermal 
instabilities may also play a role. Magnetic fields will 
generally also affect the fragmentation of the expanding 
shell, and may affect the crushing of interior clouds. No 
numerical simulation is able to contend with all these fac- 
tors: the required spatial resolution, range in timescales. 



and number of physical processes involved are far too de- 
manding. Instead, simulators must resort to the art of 
approximation. 

Several estimates of the possible impact of galactic 
winds on the IGM have been made. An instantaneous 
stellar burst of IO^Mq, typical of a dwarf galaxy with 
total halo mass ~ IO^Mq and allowing for a conver- 
sion efhciency /, ~ 0.01 of gas mass into stars, produces 
a mechanical luminosity provided by supernov ae and stel- 
lar winds of 5 X 10^3 W (|Leitherer et aLl . ll999D . assuming 
a Salpeter Initial Mass Function between 1 — IOOMq. The 
contribution of supernovae dies off after about 4 x 10^ yr, 
and the mechanical luminosity plummets, having yielded 
a total mechanical energy input of 2 x 10^* J. If the wind 
is allowed to expand until it reaches pressure equilibrium 
with the surrounding IGM, then a maximum distance the 
wind travels may be estimated. For example, if a fraction 
Vm of the mechanical energy is available to power the 
wind (taking into account losses internal to the galaxy 
and as the wind expands), balancing the mean internal 
pressure Ph = {2/'i)Ei,/Vs = i?h/27ri?f (assuming 7 = 
5/3 for the ionized interior gas), where Vg is the bubble 
volume, with the pressure of the external IGM Pjgm — 
5 X 10~^°T4(1 -I- z)'^Nm~^, gives a maximum comoving 

bubble radius of Rs,c ^ 300{iym/0-lV^^T~^^^ kpc. 

The maximum radius is comparable to the geomet- 
ric estimates discussed above. Taking account of radia- 
tive losses, however, stalls the expansion b efore it reaches 
the maximum value ( Madau et all . l200lt) . More massive 
galaxies driving stronger winds might be appealed to, but 
this is problematic on two grounds: they are much rarer 
than dwarf galaxies at high redshifts and the ram pres- 
sure of the surroundin g infalling IGM ca n smother the 
expansion of the wind ( Fuiita et a/.l . l2003 ). 

An estimate of the maximum mass halo able to pro- 
duce an escaping wind may be made as follows. Galaxies 
have nearly flat rotation curves, corresponding to en- 
velope d ensity profiles slightly steeper than p ^ l/r^ 
(jPeebled . |i993). The infall velocity of the gas is then 
VH ^ i2GMv/rvy^^, where My = (47r/3)/9yr^ is the 
mass of the virialized halo of mean density pv = 2[3(37r-|- 
2)/4]2pcrit(0)(H-2:y)3 and radius ry f^ lV.Bp . Assuming 
the gaseous component follows the profile of the envelope 
of dark matter, the ram pressure at the surface of the 
bubble is Pram = pvjj = (3/27r)(r!fc/r!„0GM2/(4p2) 
for Rs ^ ry. If the ram pressure dominates the IGM 
pressure, then it will determine the stalling radius where 
it matches Pb, according to Eq. (|126p . The criterion 
Pram < Pfc < i?h/(27rry P^) (for Rs > ry) imposes the 
upper limit on the virial mass of the galaxies able to pro- 
duce escaping winds of 



Afy < (5 X WH-'^Mq) 



^ X 3/2 /, , X -3/2 



0.001 



10 



(128) 

i n good agre e ment with hydrodynaniical computations 
(|Fuiita et all \2Q04 . For a more complete treatment, 
hydrodynamical simulations are required. Exam- 
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ales of recent wind simulations 



2002bD. lAguirre et aLl (1200 



20051). 



in clude iTheuns et al. 
Cen et al\ (^2005^. 



ICen and Ostrikei (20M}^^ lOppenheimer and Dav€ 

( 2006 r iKobavashi et al. n2003), and iKawata and RauchI 
20071) . 



IX. PROSPECTS FOR THE FUTURE 

The future of IGM studies is entering a new phase in 
both observations and computations that promises to set- 
tle a number of outstanding questions. Prominent among 
these are the circumstances that led to the reionization 
of the Universe. At the same time, the growing recogni- 
tion of the relationship between the IGM and galaxies is 
likely to shift the subject increasingly toward a merger 
with the more embracing subject of galaxy formation. 
Perhaps most spectacularly, imminent and planned de- 
tectors and space missions may open up a new era in 
IGM science: direct imaging of the structure of the IGM 
in the infrared and radio. 

Recent years have witnessed the detection of high red- 
shift Lya emitting galaxies, long-sought as candidates for 
the first forming galaxies (,Partridgc and P eebles, 19671 ) . 
Measurements of the Lya emitter luminosity function 
and its evolution would in principle provide evidence for 
the reionization epoch, as the Lya radiative damping 
edge of neutral patches would increasingly erode the lines 
and produce a sharp decrease in the luminosity func- 
tion. Evidence for such rapid evolution over the red- 
shift range 5.7 < z < 7 has been o btained from the 
Suba ru Deep Field ( Ive et all . l2006l : iKashikawa et all . 
I2OO6I ). The interpretation of any decline, however, is 
complicated by systematics, such as the unknown evo- 
lution in the dust content, the structure of the local 
IGM that may produce a Lya radiative damping wing, 
and the mass of the systems producing star forma- 
tion, as well as the possible infall of the IGM towards 
the source and the impact of winds produced by the 
source on the surrounding IGM, all of which could af- 
fect the intensity and shape of the Lya emission line 



(iDave et a/.l, 12006 ;lHaiman and Cen 


.l2005HNilsson et all 


12007; Santod. 2004). McQuinn et al. 


( 20071) suggest these 



uncertainties may be circumvented by using the angular 
correlation function of the Lya emitters. The angular 
correlations will be strongly modulated by H 11 regions 
during the EoR on large scales where the dark matter 
fluctuations are weak, providing a distinctive signature 
of reionization. Much improved Lya emitter statistics 
are anticipated with current or upcoming narrow band 



(iGardner et al. 




20061 


Horton et al\ 


I2OO8I: Stark et al 


.. 2007b). 



Anticipated CMB polarization measurements are ex- 
pected to further narrow the redshift range during which 
reionization occurred. Polarization anisotropy measure- 



ments on arcminute scales by the Planck satellite, to 
be launched in 2008, may reveal the expected patchiness 
of reionization, which may be used to discriminate be- 
tween competing re ionizatio n scenario s uinn et all 
l2005t ISantos et all i2003; Zahn et all \2QQ$ The Ata- 
cama Cosmology Telescope (ACT)^'* and the South Pole 
Telescope (SPT)^^ have already embarked on such mea- 
surements. 

The relationship between the IGM and galaxies will be 
further clarified by upcoming galaxy surveys, especially 
at low redshifts. UV spectroscopy with the Cosmic Ori- 
gins Spectrograph (COS),^^ to be installed on HST in 
2008, will refine the measurements of low to moderate 
redshift metal absorption systems, and increase sample 
sizes by an order of magnitude, greatly improving the 
statistics on the correlations with galaxies. Currently de- 
tections of the WHIM have been confined almost entirely 
to O VI measurements. Detections of further high ioniza- 
tion species would be possible using COS as well as the 
spectroscopic capabilities of the proposed x-ray satellite 
Constellation- X or a possible smaller scale sat ellite with 
high resolution x-ray spectroscopic capability ( Nicastrol . 
|2006( ). providing important diagnostics of the physical 
properties of the WHIM. LOFAR is expected to detect 
the synchrotron emission from lO*" starburst galaxies, 
which with follow-up optical, UV and sub-millimeter ob- 
servations, will provide a wealth of data that may estab- 
lish the impact of winds on the IGM over a wide range of 
redshifts, including the spreading of metals. Even greater 
numbers may be detected by a Square Kilometer Array, 
with an expected detection rate of several thousand at 
2; > 5 per square degree surveyed. Optical spectroscopy 
using a Giant Segmented Mirror Telescope (GSMT),^'' 
such as the Thirty Meter Telescope (TMT)^^ and the 
Giant Magellan Telescope (GMT),^^ and an Extremely 
Large Telescope (ELT)'^° could measure the temperature 
of the moderate and even low density components of 
the IGM by measuring the velocity widths of multiple 
metal line absorption systems. Deep photometric sur- 
veys with such very large telescopes could detect very low 
luminosity galaxies that may be associated with low red- 
shift metal absorption systems. The Ultra Deep Survey 
of the James Webb Space Telescope, to be launched 
in 2013, will detect dwarf galaxies out to redshifts of 
6 < 2; < 10, and so may positively ident ify the popula- 
tion of galaxies which reionize d the IGM ( Gardner et all 
I2OO6I : [Windhorst et all |2006|) . These multiband mea- 
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surements, extending from the radio through x-rays, will 
help to establish the metal pollution patterns throughout 
the IGM and their relationship with star-forming galax- 
ies. 

These advances will enable simulations to provide a 
complete description of the evolution of the IGM. Such 
a prospect would open up new opportunities for other 
areas of astrophysics and cosmology to exploit the IGM 
as a tool for testing models of interest. For example, the 
abundances and distribution of metals as well as the line 
widths of metal and hydrogen absorption features around 
galaxies could place stringent constraints on galaxy for- 
mation models and the winds young galaxies may pro- 
duce, including the possible impact on th e host galaxy 
and i ts surroundings of an embedded AGN ( Siiacki et all 
I2007D . The imaging of the IGM around QSO sources 
could directly reveal the opening angles of their emission 
cones and relative orientation to the observer, placing 
strong constraints on unification models of AGN. 

Precise predictions for the absorption statistics of 
the IGM may provide valuable constraints on modifi- 
cations of standard cosmological models such as those 
allowing for a "dark energy" equation of state. Dark 
energy will affect the growth rate of structures in 
the I GM, possib l y at a detectable level in future sur- 
veys (IViel et al\ . l2003bD . The density fluctuations in 
the IGM along the line of sight and transvere direc- 
tions will differ as a consequence of the difference be- 
tween radial and ang u lar d istances at high redshifts 
( Alcock and Paczvnskil Il979l ). With sufficient statis- 
tics, the difference may be detectable in the correlation 
function or power spectrum measurem ents of the trans- 
mitted Lya flux th ro ugh the IGM dHui et all. 11998 : 
Marble et all . l2008al lbl; iMcDonald and Miralda-Escude . 
19991) . Even more ambitious would be the direct de- 



tection of the expansion of the Universe through the 
monitoring of absorption lines along fixed lines of sight 
over a period of several years through spectroscopy 
on extremely lar ge telescopes (iCorasaniti et all . l2007t 
iLiske et all i2008l ). Since realizing such a measurement 
would require the development of an extraordinarily sta- 
ble detection system and adequate taming of systemat- 
ics associated with variations in the QSO brightness and 
spectral shape and their effects on the noise statistics, 
such a project may spur efforts to investigate the techni- 
cal characteristics of large telescopes and the behavior of 
QSOs at unprecedented levels of accuracy. 

The next few years may produce the first images of 
the IGM. Ionizing radiation reprocessed by the neutral 
gas in Lyman Limit Systems into Lya photons may 
soon be detectable by upcorning Lya narrow-band detec- 



tors fCantalupo et aV., 2007 



l2005HGould and Weinberg 



[1996; Hogan and WcvmaniJ. 1l987i) . It my even be pos- 



sible to detect the Lya photons emitted by a central 
source like a QSO and scattere d by the surrounding neu - 
tral IGM prior to reionization ( Loeb and RvbickH . Il999l ) 
or Lya photons produced by coUisional excitation within 
the rapidly expanding ionization fronts during the EoR 



( Cantalupo et all . l2008l ). providing Lya tomographic 
imaging of the neutral IGM. Such detections would com- 
plement the planned IGM imaging by the LOFAR, MWA 
and SKA radio facilities using 21cm tomography. Cross- 
correlations between the Lya and radio images may help 
to establish statistically significant detections. 

At even higher redshifts, 21cm tomography prior to 
reionization may become possible through the production 
of a 21cm signature via Lya and coUisional decoupling of 
the spin temperature from the CMB. At z > 30, coUi- 
sional decoupling alone is adequate, offering the oppor- 
tunity to explore the redshift desert 30 < z < 200 and ex- 
amine the primordial density fluctuation power spectrum 
on scales smaller than may be reached by CMB measure- 
ments. Such measurements could probe the structure of 
an inflati onary potential and induc e d non-gaussian flue 
tuatiqns ( Barkana and Loea . 2005 ^ Bharadwai and Ali 



2004': 'Loeb and Zaldarriagal ^2004!: iMcOu inn e^ aZ.1 . 12006 : 
Pillepich et ai, 200_7.). Astrophysical and instrumental 
systematics and sources of radio frequency interference 
(RFI), both man made and natural, such as arises from 
the ionosphere, likely preclude su ch measurements usin g 
currently planned radio facilities (jBowman et aLl . l2007f ). 
These limitations, however, may be circumvented by 
placing a radio telescope in space. Even more adven- 
turous would be to build a much larger facility on the 
far side of the Moon, where it would be shielded from 



terrestrial RFI (jCarilh et all . l2007l ). The exploration of 
long wavelengths unavailable to terrestrial radio astron- 
omy would open up a new panorama on the Universe, 
introducing a new era of "Astronomy on the Far Side." 
While not all of these prospects may materialize, others 
not currently envisaged will surely arise, ensuring that 
IGM science will continue to flourish for many years to 
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